Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 828,749 A→T 5.4% L78Q (CTG→CAG)  ybhF ← putative ABC exporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913828,7490AT 5.4% 303.2 / 5.2 95L78Q (CTG→CAG) ybhFputative ABC exporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (34/55);  new base T (3/2);  total (38/57)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

AGACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGT  >  NC_000913/828601‑828897
                                                                                                                                                    |                                                                                                                                                    
agaCGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAg                                                                                                                                                     >  2:830613/1‑150 (MQ=255)
     tAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAgcgcgc                                                                                                                                                <  1:574160/150‑1 (MQ=255)
        aCTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGCGCTGCACGCcgt                                                                                                                                             <  1:142671/150‑1 (MQ=255)
          tCCAGCAGGAGAGCAAAATTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGTTTCTACATCACCGTTAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCcgtcg                                                                                                                                           <  2:449332/150‑1 (MQ=255)
           ccAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCcgtcgt                                                                                                                                          <  2:882482/150‑1 (MQ=255)
           ccAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCcgtcgt                                                                                                                                          <  2:1047387/150‑1 (MQ=255)
             agcagGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGCGCTGCACGCcgtgctcg                                                                                                                                        >  1:217761/1‑145 (MQ=255)
                agGCCAGCAAAAGTGTGCTTACGTGCCTCGCCGGGGACGCTGCGCAATTCCGCGTACAGATTGAGGTTCTCCATCACCGGGAGTGGTTCTTAGTGACCTACTTTCTGCGGGATATAACCGAGCAGGGCGGGGTGCGCGCCGTGGTTTTTg                                                                                                                                     <  2:1469253/150‑1 (MQ=255)
                     aGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCgg                                                                                                                                <  1:299748/150‑1 (MQ=255)
                      gCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGa                                                                                                                               >  2:1318477/1‑150 (MQ=255)
                       cAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGAt                                                                                                                              <  2:1202433/150‑1 (MQ=255)
                            gTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAg                                                                                                                         <  1:807827/150‑1 (MQ=255)
                            gTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAg                                                                                                                         <  1:997420/150‑1 (MQ=255)
                             tttGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGtttt                                                                                                                                       >  2:1551805/1‑135 (MQ=255)
                             tttGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGtttt                                                                                                                                       <  1:1551805/135‑1 (MQ=255)
                             tttGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGc                                                                                                                        <  1:123527/150‑1 (MQ=255)
                                   tACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCa                                                                                                                  >  1:970351/1‑150 (MQ=255)
                                    aCGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCAc                                                                                                                 >  2:28390/1‑150 (MQ=255)
                                        gCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTg                                                                                                             <  1:901993/150‑1 (MQ=255)
                                          cTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGc                                                                                                           >  1:1629059/1‑150 (MQ=255)
                                          cTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGc                                                                                                           <  2:582312/150‑1 (MQ=255)
                                          cTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGc                                                                                                           <  1:672373/150‑1 (MQ=255)
                                          cTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGc                                                                                                           <  1:267584/150‑1 (MQ=255)
                                               cGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGc                                                                                                      <  1:1120783/150‑1 (MQ=255)
                                                 gTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCg                                                                                                    <  1:1095597/150‑1 (MQ=255)
                                                  tGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGc                                                                                                   <  2:904223/150‑1 (MQ=255)
                                                  tGACGCTGCGCAAATCCGCGTACAGATTGAGGTGCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGACGc                                                                                                   <  2:485204/150‑1 (MQ=255)
                                                   gACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCt                                                                                                  >  1:939454/1‑150 (MQ=255)
                                                    aCGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTg                                                                                                 <  2:1294883/150‑1 (MQ=255)
                                                    aCGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTg                                                                                                 >  1:1366346/1‑150 (MQ=255)
                                                        tGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCgg                                                                                             <  1:343303/150‑1 (MQ=255)
                                                          cgcAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGt                                                                                           <  2:60443/150‑1 (MQ=255)
                                                             aaaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTc                                                                                        <  2:676342/150‑1 (MQ=255)
                                                              aaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCa                                                                                       <  2:257429/150‑1 (MQ=255)
                                                              aaTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCa                                                                                       >  1:1670586/1‑150 (MQ=255)
                                                               aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGTATCAAAGCCAATCAACGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                      <  2:1467780/150‑1 (MQ=255)
                                                               aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                      <  2:348060/150‑1 (MQ=255)
                                                               aTCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAg                                                                                      <  2:1403973/150‑1 (MQ=255)
                                                                tCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGt                                                                                     <  2:1629059/150‑1 (MQ=255)
                                                                  cgcgTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATAAAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTaa                                                                                   >  2:563339/1‑150 (MQ=255)
                                                                    cgTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATc                                                                                 <  1:1358501/150‑1 (MQ=255)
                                                                            ttGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAAc                                                                         >  1:894399/1‑150 (MQ=255)
                                                                                 gTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCg                                                                    <  2:1608456/150‑1 (MQ=255)
                                                                                 gTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCg                                                                    <  1:292260/150‑1 (MQ=255)
                                                                                  ttCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGc                                                                   <  1:1588949/150‑1 (MQ=255)
                                                                                   tctcCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCa                                                                  <  1:842335/150‑1 (MQ=255)
                                                                                   tctcCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCa                                                                  <  1:167170/150‑1 (MQ=255)
                                                                                   tctcCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCa                                                                  >  1:1037360/1‑150 (MQ=255)
                                                                                      ccATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCa                                                               >  2:1267393/1‑150 (MQ=255)
                                                                                          cACCGTGAGATCTTCATACAGACCAAATTTCTGAGGCATATAACCGAGCAAGGCGTGCAGCGCGCCGTAGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGt                                                           >  2:748314/1‑150 (MQ=255)
                                                                                             cGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGt                                                        <  2:1366346/150‑1 (MQ=255)
                                                                                               tGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGttt                                                      <  2:876678/150‑1 (MQ=255)
                                                                                                  gaTCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGAGCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTAc                                                   <  1:677543/150‑1 (MQ=255)
                                                                                                   aTCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTAcc                                                  <  2:970351/150‑1 (MQ=255)
                                                                                                    tCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCt                                                 >  2:1477134/1‑150 (MQ=255)
                                                                                                     cTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTg                                                <  2:1180858/150‑1 (MQ=255)
                                                                                                     cTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTg                                                >  2:1549577/1‑150 (MQ=255)
                                                                                                     cTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTg                                                <  2:910483/150‑1 (MQ=255)
                                                                                                      ttCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGc                                               <  2:692627/150‑1 (MQ=255)
                                                                                                       tCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCa                                              >  1:1005456/1‑150 (MQ=255)
                                                                                                         aTACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCAcc                                            >  1:1515129/1‑150 (MQ=255)
                                                                                                             aGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTcc                                        <  1:606519/150‑1 (MQ=255)
                                                                                                             aGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTcc                                        <  1:1556734/150‑1 (MQ=255)
                                                                                                              gACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGGCGTTTTTGATCGGATCAAAGCCAATCCCCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCGAACTTCCGCATCAGCGTGGTTTTACCTGCACCGTCCg                                       >  2:1047085/1‑150 (MQ=255)
                                                                                                                ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                     >  1:265460/1‑150 (MQ=255)
                                                                                                                ccAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGc                                     <  1:1212488/150‑1 (MQ=255)
                                                                                                                 cAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGcc                                    <  2:507320/150‑1 (MQ=255)
                                                                                                                  aaaTTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGccc                                   >  2:1266129/1‑150 (MQ=255)
                                                                                                                  aaaTTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGccc                                   >  1:548078/1‑150 (MQ=255)
                                                                                                                    aTTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCcca                                 >  2:104474/1‑150 (MQ=255)
                                                                                                                      ttCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCccacc                               <  1:28390/150‑1 (MQ=255)
                                                                                                                      ttCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCccacc                               <  1:1549577/150‑1 (MQ=255)
                                                                                                                          gCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGcc                                                                            <  2:1252089/101‑1 (MQ=255)
                                                                                                                          gCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGcc                                                                            >  1:1252089/1‑101 (MQ=255)
                                                                                                                           cGGCATATAACCGAGCACGGCGTGCAGCGCGCCGCCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAAccc                          >  1:1033214/1‑150 (MQ=255)
                                                                                                                               atatAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTc                      <  2:214031/150‑1 (MQ=255)
                                                                                                                               atatAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTc                      <  2:548078/150‑1 (MQ=255)
                                                                                                                                tataACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTca                     <  1:98610/150‑1 (MQ=255)
                                                                                                                                  taACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTccca                   >  2:1611034/1‑150 (MQ=255)
                                                                                                                                  taACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTcaca                   >  1:1404011/1‑150 (MQ=255)
                                                                                                                                  taACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTcaca                   >  1:1409805/1‑150 (MQ=255)
                                                                                                                                   aaCCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACAt                  >  2:425493/1‑150 (MQ=255)
                                                                                                                                   aaCCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACAt                  >  1:1637879/1‑150 (MQ=255)
                                                                                                                                    aaaaacgaCGGCGCGCTGCACGCCGTCGTTTTTGATCGGATCTAAGCCAATCACCGTGGCCCTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATa                 >  1:676584/8‑150 (MQ=255)
                                                                                                                                        acgaCGGCGCGCCGCACGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCg             >  2:1591561/4‑150 (MQ=255)
                                                                                                                                         gCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCgg            >  1:1646035/1‑150 (MQ=255)
                                                                                                                                         cgaCGGCGCGCTGCACGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCgg            >  2:445914/3‑150 (MQ=255)
                                                                                                                                           aCGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCg          >  1:606716/1‑150 (MQ=255)
                                                                                                                                            cGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGt         >  1:587449/1‑150 (MQ=255)
                                                                                                                                             ggCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTg        <  2:959463/150‑1 (MQ=255)
                                                                                                                                             ggCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTg        <  1:1266129/150‑1 (MQ=255)
                                                                                                                                             ggCGTGCAGCGCGCCGTCGTTTTTGATAGGATCAAAGCCAATCACCGTGCCAATGCCGCTGTAGGTTTTCAGTAATCCCGCCAACATCCGAATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGACACATAACCGGCGTg        >  2:1194143/1‑150 (MQ=255)
                                                                                                                                                  gCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATgg   <  2:1005456/150‑1 (MQ=255)
                                                                                                                                                   cAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGt  <  2:1037360/150‑1 (MQ=255)
                                                                                                                                                   cAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGt  <  2:12667/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGACGTAAACTCCAGCAGGCGAGCAAAAGTTTGCTTACGTGCCTCGCCGGTGACGCTGCGCAAATCCGCGTACAGATTGAGGTTCTCCATCACCGTGAGATCTTCATACAGACCAAATTTCTGCGGCATATAACCGAGCACGGCGTGCAGCGCGCCGTCGTTTTTGATCGGATCAAAGCCAATCACCGTGGCACTGCCGCTGTCGGGTTTCAGTAATCCCGCCAACATCCGCATCAGCGTGGTTTTACCTGCACCGTCCGGCCCCACCAACCCCGTCACATAACCGGCGTGAATGGT  >  NC_000913/828601‑828897

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