Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 380,105 A→C 7.3% G246G (GGT→GGG)  yaiO ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913380,1050AC 7.3% 334.4 / 10.9 100G246G (GGT→GGG) yaiOouter membrane protein
Reads supporting (aligned to +/- strand):  ref base A (53/39);  new base C (6/2);  total (59/41)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTAAG  >  NC_000913/379962‑380253
                                                                                                                                               |                                                                                                                                                    
acTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTa                                                                                                                                                >  1:2704324/1‑150 (MQ=255)
  tGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTaaa                                                                                                                                              >  2:2296734/1‑150 (MQ=255)
                      gTGAATGTGGTCGGAAAACAAAGAGGAAA‑GGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTa                                                                                                                         <  2:541776/150‑1 (MQ=255)
                        gAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTAc                                                                                                                        >  2:1411023/1‑150 (MQ=255)
                                    acaacaacGAGGCACGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt                                                                                                            <  2:1654440/148‑1 (MQ=255)
                                    aaaaCACCGACGACAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCt                                                                                                            <  2:1342523/150‑1 (MQ=255)
                                         aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt                                                                                                       >  1:2677151/1‑150 (MQ=255)
                                         aaaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGt                                                                                                       >  1:2384270/1‑150 (MQ=255)
                                          aaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc                                                                                                      >  2:2576927/1‑150 (MQ=255)
                                          aaGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTc                                                                                                      >  2:1180182/1‑150 (MQ=255)
                                             agGAAAGGGGGGGGGGCTAATCGGCAGGTAAGGCCGCCCCGGATAGCGGGCGGCAGTAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCATACCGTTAAAATCATCGGTTGGGGTTTCGTACCACACTTTACCTGCCGTCAGc                                                                                                   >  2:534470/1‑150 (MQ=255)
                                             agGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGc                                                                                                   >  1:2313959/1‑150 (MQ=255)
                                                 aaagggggggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCga                                                                                               <  2:1864963/147‑1 (MQ=255)
                                                 aaGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCga                                                                                               <  2:1079538/150‑1 (MQ=255)
                                                      gggggggCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTa                                                                                          >  1:1913403/1‑150 (MQ=255)
                                                                      cgggACGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGGCAGACGGGATGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCg                                                                          <  1:1248548/149‑1 (MQ=255)
                                                                      aGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCg                                                                          <  1:1248533/150‑1 (MQ=255)
                                                                      aGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCg                                                                          >  2:1539909/1‑150 (MQ=255)
                                                                       gggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCgg                                                                         <  2:415206/150‑1 (MQ=255)
                                                                       gggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCgg                                                                         <  1:539702/150‑1 (MQ=255)
                                                                       gggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCgg                                                                         <  1:1338161/150‑1 (MQ=255)
                                                                       gggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCgg                                                                         <  1:1338155/150‑1 (MQ=255)
                                                                        ggAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGt                                                                        >  2:949884/1‑150 (MQ=255)
                                                                              ccgccCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATa                                                                  >  2:1457715/1‑150 (MQ=255)
                                                                                  cccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGtt                                                              >  1:793835/1‑150 (MQ=255)
                                                                                   ccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTCCCTGCCGTCAGCCCGAGATTAAGGTGCTCAGTCAGCGGTTGAATACGTTg                                                             <  2:1234562/150‑1 (MQ=255)
                                                                                   ccGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTg                                                             <  1:1709589/150‑1 (MQ=255)
                                                                                   ccGGAGAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTACCATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTg                                                             <  2:1830190/150‑1 (MQ=255)
                                                                                    cGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGc                                                            >  2:1681037/1‑150 (MQ=255)
                                                                                    cGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGc                                                            >  2:1680856/1‑150 (MQ=255)
                                                                                       aTAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGa                                                         <  2:2059653/150‑1 (MQ=255)
                                                                                         aGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACt                                                       >  1:1641338/1‑150 (MQ=255)
                                                                                         aGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACt                                                       >  1:1267710/1‑150 (MQ=255)
                                                                                         aGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACt                                                       <  1:2448124/150‑1 (MQ=255)
                                                                                         aGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACt                                                       >  1:327066/1‑150 (MQ=255)
                                                                                          gcgggcggCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTg                                                      <  2:2020196/150‑1 (MQ=255)
                                                                                          gcgggcggCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTg                                                      <  2:489411/150‑1 (MQ=255)
                                                                                          gcgggcggCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTg                                                      <  2:489404/150‑1 (MQ=255)
                                                                                          gcgggcggCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTg                                                      <  2:458213/150‑1 (MQ=255)
                                                                                          gcgggcggCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTg                                                      <  2:1391441/150‑1 (MQ=255)
                                                                                                 gcAGAAGGAATCAGAATGTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGctct                                               <  1:1669569/150‑1 (MQ=255)
                                                                                                  cAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCtt                                              >  1:2580825/1‑150 (MQ=255)
                                                                                                  cAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCtt                                              >  1:1754251/1‑150 (MQ=255)
                                                                                                   aGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTc                                             <  2:1516700/150‑1 (MQ=255)
                                                                                                   aGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTc                                             <  2:1516698/150‑1 (MQ=255)
                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                                           >  2:662463/1‑150 (MQ=255)
                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                                           >  2:1768953/1‑150 (MQ=255)
                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                                           <  2:2143064/150‑1 (MQ=255)
                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                                           >  2:2406054/1‑150 (MQ=255)
                                                                                                     aaGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGGTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAt                                           >  2:854753/1‑150 (MQ=255)
                                                                                                        gactCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGCCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAtgct                                        <  2:518173/147‑1 (MQ=255)
                                                                                                         aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAtgctg                                       >  1:1599334/1‑150 (MQ=255)
                                                                                                         aaTCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCAtgctg                                       >  1:2239612/1‑150 (MQ=255)
                                                                                                           tCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGcc                                     >  1:304700/1‑150 (MQ=255)
                                                                                                           tCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGcc                                     >  1:304699/1‑150 (MQ=255)
                                                                                                            cAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCg                                    >  2:1074433/1‑150 (MQ=255)
                                                                                                              gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTa                                  >  1:1840781/1‑150 (MQ=255)
                                                                                                              gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTa                                  >  1:1751007/1‑150 (MQ=255)
                                                                                                              gAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTa                                  >  1:891439/1‑150 (MQ=255)
                                                                                                               aaTTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAg                                 >  1:1773506/1‑150 (MQ=255)
                                                                                                                  ttCCAGGTCAGTCGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCgt                              <  2:1854239/150‑1 (MQ=255)
                                                                                                                  ttCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCgt                              <  2:1854163/150‑1 (MQ=255)
                                                                                                                     cAGGTCAGTCGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGtt                           <  2:2245756/150‑1 (MQ=255)
                                                                                                                     cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGtt                           >  1:1913935/1‑150 (MQ=255)
                                                                                                                     cAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGtt                           <  2:2182906/150‑1 (MQ=255)
                                                                                                                      gattttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGttt                          >  2:2582121/8‑150 (MQ=255)
                                                                                                                        ttttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCt                        >  2:265934/6‑150 (MQ=255)
                                                                                                                        ttttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCt                        >  2:2464661/6‑150 (MQ=255)
                                                                                                                        ttttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCt                        >  2:2464644/6‑150 (MQ=255)
                                                                                                                        ttttaaCGGTCTGCAACTTGCAGCCCGTTAAAATCATCGGGTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCt                        >  2:1981092/6‑150 (MQ=255)
                                                                                                                           aGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGc                     <  2:864417/150‑1 (MQ=255)
                                                                                                                             aCGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGt                   >  2:1349592/1‑150 (MQ=255)
                                                                                                                             aCGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGt                   >  2:2457343/1‑150 (MQ=255)
                                                                                                                              cGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTc                  >  1:2467393/1‑150 (MQ=255)
                                                                                                                               gggCTGCAAGTTGCAGACCGTTAAATTCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATTCGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTcc                 >  1:1363472/1‑150 (MQ=255)
                                                                                                                               gggCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTcc                 >  1:2542963/1‑150 (MQ=255)
                                                                                                                                ggCTGCAGGTCGCTGCCCGCTACAATCATCGGCTGGGGTGGCGTACCACACTTTACCTGGCGTCTGCGCGTGATTTTTTTGCTCATTCTGCGGTTGCATACGTTGCAGACTGTCGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCa                <  2:1228504/150‑1 (MQ=255)
                                                                                                                                ggCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCa                <  2:1877251/150‑1 (MQ=255)
                                                                                                                                ggCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCa                <  2:1228500/150‑1 (MQ=255)
                                                                                                                                   tGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTc             >  1:937049/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  2:673239/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  1:785815/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  2:1897053/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  1:359015/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  1:785946/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  1:1700976/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  2:2142552/1‑150 (MQ=255)
                                                                                                                                     cAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  2:2142667/1‑150 (MQ=255)
                                                                                                                                     cAACTTGCAGCCCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTCGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCAt           >  1:2075409/1‑150 (MQ=255)
                                                                                                                                       aGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAg         <  2:1476561/150‑1 (MQ=255)
                                                                                                                                         ttGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGgt       <  2:2647197/150‑1 (MQ=255)
                                                                                                                                         ttGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGgt       <  2:221514/150‑1 (MQ=255)
                                                                                                                                         ttGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGgt       <  2:1284216/150‑1 (MQ=255)
                                                                                                                                         gtGCAGCCCGTTAACCTCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGgt       <  2:967409/149‑1 (MQ=255)
                                                                                                                                         gtGCAGACCGTTAAAATCCTCGGTTGGGGGGTCGTACCACACTTTACCCGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGgt       <  2:1814081/149‑1 (MQ=255)
                                                                                                                                             aGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTaa   <  2:2614950/150‑1 (MQ=255)
                                                                                                                                              gACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTAAg  >  1:1375191/1‑150 (MQ=255)
                                                                                                                                              gACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTAAg  >  2:899755/1‑150 (MQ=255)
                                                                                                                                              gACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTAAg  >  2:935341/1‑150 (MQ=255)
                                                                                                                                              gACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTAAg  >  2:1316004/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACTGGAGAATAAAATTTATCCGGTGAATGTGGTCGGAAAACAAAGAGGAAAGGGGGGGGGGCTAATCGGCAGGGAAGGCCGCCCCGGATAGCGGGCGGCAGAAGGAATCAGAATTTCCAGGTCAGACGGGCTGCAAGTTGCAGACCGTTAAAATCATCGGTTGGGGTGTCGTACCACACTTTACCTGCCGTCAGCCCGAGATTAAGTTGCTCAGTCAGCGGTTGAATACGTTGCAGACTGACGCTCTTCATGCTGCCGTAGCGTGTTTCTGGCGTCCAGTCATAGGTGTAAG  >  NC_000913/379962‑380253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑