Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,918,296 C→T 5.2% intergenic (‑31/+20) atpG ← / ← atpA F1 sector of membrane‑bound ATP synthase, gamma subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,918,2960CT 5.2% 472.1 / 9.0 152intergenic (‑31/+20)atpG/atpAF1 sector of membrane‑bound ATP synthase, gamma subunit/F1 sector of membrane‑bound ATP synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (76/68);  new base T (2/6);  total (78/74)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

Only 100 of 155 total aligned reads displayed.

ATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCCAG  >  NC_000913/3918150‑3918444
                                                                                                                                                  |                                                                                                                                                    
aTGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCGTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTaa                                                                                                                                                   >  2:95102/1‑150 (MQ=255)
  gCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAgg                                                                                                                                                 >  2:1089854/1‑150 (MQ=255)
   cGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGc                                                                                                                                                >  1:919827/1‑150 (MQ=255)
   cGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGc                                                                                                                                                >  1:1739409/1‑150 (MQ=255)
    gATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCa                                                                                                                                               <  1:761542/150‑1 (MQ=255)
        cTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATTAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgcc                                                                                                                                           <  1:362640/150‑1 (MQ=255)
         tGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCTTGCCTTAGGGCAGgccg                                                                                                                                          <  2:711726/150‑1 (MQ=255)
         tGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgccg                                                                                                                                          <  1:1999831/150‑1 (MQ=255)
          gCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGTACGCTTGCGCTCTTTCTACGTATCTCTTTTGCGACGGCCATGAGCTTCTCCTCACTGCCTTGCGGCCTGCCCTAAGGCAAgccgc                                                                                                                                         >  1:1233755/1‑150 (MQ=255)
          gCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAgccgc                                                                                                                                         >  1:214920/1‑150 (MQ=255)
            gATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCa                                                                                                                                       <  2:444349/150‑1 (MQ=255)
              cttACGCATTTTGGAAGCGGCGACCATCTCCCTCGCTTGAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGa                                                                                                                                     <  1:81739/149‑1 (MQ=255)
               ttACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGAc                                                                                                                                    <  1:831763/150‑1 (MQ=255)
                 aCGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGt                                                                                                                                  <  2:1190374/150‑1 (MQ=255)
                  cGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGtt                                                                                                                                 >  2:449417/1‑150 (MQ=255)
                  cGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGtt                                                                                                                                 <  2:1851848/150‑1 (MQ=255)
                   gCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTGGCGGCTTGCCTTAGGGCAGGCCGCAAggcattg                                                                                                                                <  2:1130865/150‑7 (MQ=38)
                    cATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGCGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCACGCCGCCCGACGTTAc                                                                                                                               >  2:1277330/1‑150 (MQ=255)
                       tttGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAg                                                                                                                            >  1:693504/1‑150 (MQ=255)
                            aaGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTg                                                                                                                       <  2:874710/150‑1 (MQ=255)
                             aGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCACGCCGCCAGAAGTTACCAGGATTgg                                                                                                                      >  1:160252/1‑150 (MQ=255)
                               cggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGt                                                                                                                    <  2:859272/150‑1 (MQ=255)
                                ggcgACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGtt                                                                                                                   >  1:1522949/1‑150 (MQ=255)
                                     ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCttt                                                                                                              <  2:1387962/150‑1 (MQ=255)
                                     ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCttt                                                                                                              >  1:2051945/1‑150 (MQ=255)
                                     ccATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGACCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCttt                                                                                                              <  1:1005124/150‑1 (MQ=255)
                                       aTCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGa                                                                                                            <  1:919150/150‑1 (MQ=255)
                                        tctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGCCGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGaa                                                                                                           >  2:630537/1‑150 (MQ=255)
                                         ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTTCCAGGATTGGGTTGCTTTGAAg                                                                                                          >  2:1002334/1‑150 (MQ=255)
                                         ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                          <  1:639883/150‑1 (MQ=255)
                                         ctcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAg                                                                                                          >  2:2785/1‑150 (MQ=255)
                                          tcCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAgg                                                                                                         >  2:495851/1‑150 (MQ=255)
                                                gCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGATTCTCCTCAATGCCTTGCGGACTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCgc                                                                                                   >  1:1273288/1‑149 (MQ=255)
                                                 cTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAg                                                                                                  <  1:301141/150‑1 (MQ=255)
                                                  tttAGTGATCTTTTGCGTGTGCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAgg                                                                                                 <  2:694367/150‑1 (MQ=255)
                                                   ttAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGa                                                                                                >  2:1211605/1‑150 (MQ=255)
                                                     aGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCAAAGGAAAGCCGCCAGACGTTAACAGGATTGGGTTGCTTTGAAGGAATAGAGGATg                                                                                              >  1:1721294/1‑150 (MQ=255)
                                                      gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGc                                                                                             >  1:652475/1‑150 (MQ=255)
                                                      gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCAAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGc                                                                                             >  1:409577/1‑150 (MQ=255)
                                                      gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCAATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGc                                                                                             >  1:1146934/1‑150 (MQ=255)
                                                      gTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGAGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCCACCCGCCCGACGGGACCAGGATTGGGTTGCTTTTAAGGAATCGAGGCTGc                                                                                             >  1:1785047/1‑150 (MQ=255)
                                                       tGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGcc                                                                                            <  1:1706563/150‑1 (MQ=255)
                                                        gATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCt                                                                                           >  2:1766605/1‑150 (MQ=255)
                                                         aTCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCtt                                                                                          >  2:1377449/1‑150 (MQ=255)
                                                          tCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCttt                                                                                         >  1:190636/1‑150 (MQ=255)
                                                           cTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTc                                                                                        <  1:1406430/150‑1 (MQ=255)
                                                              ttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                                                                                     <  1:495851/150‑1 (MQ=255)
                                                              ttGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGc                                                                                     >  1:365367/1‑150 (MQ=255)
                                                                gCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGCCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGcgt                                                                                   >  2:1350687/1‑148 (MQ=255)
                                                                  gtgtTCTGTTCGCTTGCGATATTACTACGTATATATTTTGCGGCGGCCATGATCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGc                                                                                 <  1:2014399/150‑1 (MQ=255)
                                                                      tCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTc                                                                             <  1:1598642/150‑1 (MQ=255)
                                                                       cTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCg                                                                            >  2:65810/1‑150 (MQ=255)
                                                                       cTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCg                                                                            <  2:746589/150‑1 (MQ=255)
                                                                        tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGa                                                                           <  2:1192983/150‑1 (MQ=255)
                                                                        tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGa                                                                           <  2:1543949/150‑1 (MQ=255)
                                                                        tGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGa                                                                           <  2:286278/150‑1 (MQ=255)
                                                                         ggACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGAt                                                                          <  1:1654400/150‑1 (MQ=255)
                                                                          gACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGAGTGGGTTGCTTTGAAGGAATCGAGGATTCCTTTCAGCTTGCCTTCCAtt                                                                         >  2:778510/1‑150 (MQ=255)
                                                                              cTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgt                                                                     >  2:553388/1‑150 (MQ=255)
                                                                                tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg                                                                   >  1:1070739/1‑150 (MQ=255)
                                                                                tGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgtcg                                                                   >  2:677593/1‑150 (MQ=255)
                                                                                    aTCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGt                                                               <  2:698491/150‑1 (MQ=255)
                                                                                     tCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                                                              <  1:931956/150‑1 (MQ=255)
                                                                                      cTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAg                                                             <  2:3435/150‑1 (MQ=255)
                                                                                        tactacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc                                                           <  2:655808/150‑1 (MQ=255)
                                                                                        tactacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc                                                           >  2:21838/1‑150 (MQ=255)
                                                                                          ctacGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                                                         >  1:93300/1‑150 (MQ=255)
                                                                                             cGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCa                                                                                       >  2:1385812/1‑117 (MQ=255)
                                                                                             cGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCgg                                                      <  1:107096/150‑1 (MQ=255)
                                                                                                 tctcGTATGCGACGGCTATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGTAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGGCTg                                                  <  1:321843/150‑1 (MQ=255)
                                                                                                  ctctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTgg                                                 <  1:1170356/150‑1 (MQ=255)
                                                                                                  ctctTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAACGCAAGCCGCCAGACGTTACCATGTTTGGGTTGCTTTGAAGGAAGCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGAACCCACCGGTCTgg                                                 >  2:814918/1‑150 (MQ=255)
                                                                                                           gcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCCAGCCGCCAGACGTTACCAGGGTTGGGTTGCTTTGAAGGAATCGAGGTTGCCTTTCAGCGTGCCTTCGATTTCGTCGTTGTGGCCACCGGTCTGGTTGTtctcg                                        >  1:951230/1‑149 (MQ=255)
                                                                                                           gcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATctct                                        <  2:2012011/150‑1 (MQ=255)
                                                                                                           gcCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATctct                                        <  1:1310199/150‑1 (MQ=255)
                                                                                                             cGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTg                                      <  2:1413212/150‑1 (MQ=255)
                                                                                                             cGGCCATGAGCTTCTCCTCAAGGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTg                                      <  1:240185/150‑1 (MQ=255)
                                                                                                              ggCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGc                                     <  1:7302/150‑1 (MQ=255)
                                                                                                               gCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCTTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCa                                    <  1:1960265/150‑1 (MQ=255)
                                                                                                                 cATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATc                                  >  2:897721/1‑150 (MQ=255)
                                                                                                                     aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGTATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTTGTTGATCTCTTGCATCAACg                              <  1:470839/150‑1 (MQ=255)
                                                                                                                     aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCCGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACg                              >  2:163611/1‑150 (MQ=255)
                                                                                                                        ttCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAg                           <  1:1241521/150‑1 (MQ=255)
                                                                                                                         tctcCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGc                          >  1:522191/1‑150 (MQ=255)
                                                                                                                         gctcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCg                         <  2:1711994/149‑1 (MQ=255)
                                                                                                                          ctcctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCg                         >  2:1106056/1‑150 (MQ=255)
                                                                                                                            cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTAGATTTCGTCGTTGTAGCCACCGGGCTGGTTGTTCTAGTGCCTCAACGGAGCGTg                       >  2:132549/1‑150 (MQ=255)
                                                                                                                             ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGa                      <  1:65810/150‑1 (MQ=255)
                                                                                                                                aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa                   <  2:1070738/150‑1 (MQ=255)
                                                                                                                                    ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGt               >  1:1894647/1‑150 (MQ=255)
                                                                                                                                         gggCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg          <  2:1739408/149‑1 (MQ=255)
                                                                                                                                         cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGCTTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg          <  2:833045/150‑1 (MQ=255)
                                                                                                                                         cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg          >  1:1954051/1‑150 (MQ=255)
                                                                                                                                         cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg          <  2:825860/150‑1 (MQ=255)
                                                                                                                                         cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGGTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg          <  1:1006778/150‑1 (MQ=255)
                                                                                                                                           gCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa        >  2:1980542/1‑150 (MQ=255)
                                                                                                                                               gCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGcc    <  1:1338534/150‑1 (MQ=255)
                                                                                                                                                cccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCca   >  1:669950/1‑150 (MQ=255)
                                                                                                                                                 ccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcag  >  2:1423625/1‑150 (MQ=255)
                                                                                                                                                 ccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTACGGCGATGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCcag  >  1:98196/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATGCGATCCTGCGATTTACGCATTTTGGAAGCGGCGACCATCTCCATCGCTTTAGTGATCTTTTGCGTGTTCTGGACGCTTGCGATCTTACTACGTATCTCTTTTGCGCCGGCCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGCCAG  >  NC_000913/3918150‑3918444

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