Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 934,195 T→A 7.2% A324A (GCT→GCA ftsK → cell division DNA translocase FtsK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913934,1950TA 7.2% 250.4 / 8.1 83A324A (GCT→GCAftsKcell division DNA translocase FtsK
Reads supporting (aligned to +/- strand):  ref base T (43/34);  new base A (3/3);  total (46/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

CGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTG  >  NC_000913/934048‑934343
                                                                                                                                                   |                                                                                                                                                    
cGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAgctgc                                                                                                                                                    >  1:287226/1‑150 (MQ=255)
  aaGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAgctgctg                                                                                                                                                  >  2:85827/1‑150 (MQ=255)
   agagATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAgctgctgc                                                                                                                                                 >  1:1538715/1‑150 (MQ=255)
   agagATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAgctgctgc                                                                                                                                                 <  1:1808577/150‑1 (MQ=255)
    gagaTTACCTACACTGCACGCGTTGTGGCTGCCGACCCTGACGACGTACTATTTTCGGTCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACTGTGCGCCAATTACCGAACCTGTCGCTGTAGCAgctgctgct                                                                                                                                                <  2:773357/150‑1 (MQ=255)
    gagaTTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAgctgctgct                                                                                                                                                <  1:1428040/150‑1 (MQ=255)
      gaTTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTAc                                                                                                                                              >  1:306640/1‑150 (MQ=255)
        ttACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCa                                                                                                                                            >  2:1215893/1‑150 (MQ=255)
         tACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCAc                                                                                                                                           >  1:967956/1‑150 (MQ=255)
         tACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCAGCTGCTAcagc                                                                                                                                           <  1:1818992/150‑4 (MQ=255)
         tACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCAGCTGCTAcagc                                                                                                                                           <  1:59401/150‑4 (MQ=255)
             tACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCg                                                                                                                                       >  1:612392/1‑150 (MQ=255)
                 cTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTTCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGacac                                                                                                                                   >  1:1418241/1‑150 (MQ=255)
                 cTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGacac                                                                                                                                   >  1:713132/1‑150 (MQ=255)
                   gCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACaa                                                                                                                                 >  2:1353025/1‑150 (MQ=255)
                    cACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCCGTGGTAGCAGCAGCTGCTACAGCGacaggtt                                                                                                                                <  1:713286/150‑5 (MQ=38)
                           tgtgGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGc                                                                                                                         <  1:1215893/150‑1 (MQ=255)
                           tgtgGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGc                                                                                                                         <  2:967956/150‑1 (MQ=255)
                              ggCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTgc                                                                                                                      >  2:978646/1‑150 (MQ=255)
                              ggCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTgc                                                                                                                      >  2:1120617/1‑150 (MQ=255)
                                   ccGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCgg                                                                                                                 <  2:306640/150‑1 (MQ=255)
                                    cGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGt                                                                                                                <  1:217423/150‑1 (MQ=255)
                                        ccGGACGCCGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAAGTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGaa                                                                                                            >  2:761061/1‑150 (MQ=255)
                                             cgacgTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCtgt                                                                                                       <  1:865123/150‑1 (MQ=255)
                                              gacgTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCtgtg                                                                                                      <  1:417546/150‑1 (MQ=255)
                                                  tCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTc                                                                                                  >  2:1591265/1‑150 (MQ=255)
                                                   ccTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCa                                                                                                 >  1:903630/1‑150 (MQ=255)
                                                     tATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGa                                                                                               <  1:1646749/150‑1 (MQ=255)
                                                      aTTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGAc                                                                                              >  1:1602759/1‑150 (MQ=255)
                                                          tCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGAcgccg                                                                                          >  1:872517/1‑150 (MQ=255)
                                                            gggCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACgccgcc                                                                                        <  2:872517/150‑1 (MQ=255)
                                                            gagCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACgccgcc                                                                                        >  1:1966437/3‑150 (MQ=255)
                                                              gCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTg                                                                                      <  2:222261/150‑1 (MQ=255)
                                                               cAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGt                                                                                     >  2:425307/1‑150 (MQ=255)
                                                                 aTCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTg                                                                                   <  1:85827/150‑1 (MQ=255)
                                                                 aTCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTg                                                                                   >  1:1053418/1‑150 (MQ=255)
                                                                 aTCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTg                                                                                   <  2:713132/150‑1 (MQ=255)
                                                                   cGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGcc                                                                                 <  1:1320903/150‑1 (MQ=255)
                                                                      gCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCctct                                                                              >  1:196382/1‑150 (MQ=255)
                                                                      gCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCctct                                                                              >  1:1953182/1‑150 (MQ=255)
                                                                       cAACGCAGCCAGAATATGACTAATACCATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTg                                                                             >  2:1840622/1‑150 (MQ=255)
                                                                        aaCGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGt                                                                            <  2:1538715/150‑1 (MQ=255)
                                                                        aaCGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGt                                                                            >  2:813525/1‑150 (MQ=255)
                                                                         aCGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGtt                                                                           >  1:315955/1‑150 (MQ=255)
                                                                         aCGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGtt                                                                           <  1:1945065/150‑1 (MQ=255)
                                                                          cGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAAc                                                                                                           >  1:1423343/1‑117 (MQ=255)
                                                                          cGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAAc                                                                                                           <  2:1423343/117‑1 (MQ=255)
                                                                           gCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGa                                                                         >  1:1424067/1‑150 (MQ=255)
                                                                              gCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGt                                                                      >  2:1368963/1‑150 (MQ=255)
                                                                                  gAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCa                                                                  <  2:287226/150‑1 (MQ=255)
                                                                                  gAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCa                                                                  <  2:159412/150‑1 (MQ=255)
                                                                                   aaTATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCAc                                                                 <  2:1966437/150‑1 (MQ=255)
                                                                                      atGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTg                                                              >  2:1783578/1‑150 (MQ=255)
                                                                                       tGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTgc                                                             >  1:1293535/1‑150 (MQ=255)
                                                                                         aCGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTgcgc                                                           >  2:360481/1‑150 (MQ=255)
                                                                                          cGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCa                                                          >  2:1548092/1‑150 (MQ=255)
                                                                                              tACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCt                                                      <  1:1851695/150‑1 (MQ=255)
                                                                                                cGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTAc                                                    >  2:723964/1‑150 (MQ=255)
                                                                                                 gATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACa                                                   >  1:16521/1‑150 (MQ=255)
                                                                                                   tCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGt                                                 >  1:1377451/1‑150 (MQ=255)
                                                                                                   tCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGt                                                 <  1:425307/150‑1 (MQ=255)
                                                                                                     cATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAg                                               >  1:937042/1‑150 (MQ=255)
                                                                                                       ttattaAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGcc                                             <  2:315955/150‑1 (MQ=255)
                                                                                                           taAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGtt                                                                     >  1:1542996/1‑122 (MQ=255)
                                                                                                           taAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGtt                                                                     <  2:1542996/122‑1 (MQ=255)
                                                                                                             aaCGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAg                                       <  2:612392/150‑1 (MQ=255)
                                                                                                               cGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGcc                                     <  1:1548092/150‑1 (MQ=255)
                                                                                                               cGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGcc                                     <  1:723964/150‑1 (MQ=255)
                                                                                                               cGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGcc                                     <  1:546452/150‑1 (MQ=255)
                                                                                                                         aTTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGt                           <  1:1840622/150‑1 (MQ=255)
                                                                                                                           tACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTcc                         <  1:978646/150‑1 (MQ=255)
                                                                                                                           tACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTcc                         >  2:341816/1‑150 (MQ=255)
                                                                                                                            aCCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCca                        >  2:688483/1‑150 (MQ=255)
                                                                                                                            aCCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCca                        >  1:749141/1‑150 (MQ=255)
                                                                                                                            aCCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCAGGTTGAGCCTGTGACTCAGACGCCGCCTGTTGCCTCTGGTGATCTTCCACCTGCGCAACCTACAGTAGCCTGCCAGCCTGAACCGGGTCca                        >  2:996358/1‑150 (MQ=255)
                                                                                                                              gctttgtGTCGCCGTGGTAGCAGCAGCTGCTACAGCGACACAAAGCTGGGCTGCGCCGGGTGAACCTGGGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTACACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCcaca                      >  2:1049455/7‑150 (MQ=37)
                                                                                                                                 aCCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAAc                   <  2:967147/150‑1 (MQ=255)
                                                                                                                                    cgccgtggTAGCAGCAGCTGCTACA‑‑GCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGAATCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTAACGGGTACACAAACGGgaga             >  2:846459/8‑150 (MQ=255)
                                                                                                                                      tCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGgag              <  1:1909323/150‑1 (MQ=255)
                                                                                                                                          tGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGcc          >  2:1548630/1‑150 (MQ=255)
                                                                                                                                          ggTAGCAGCAGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGcc          >  2:456672/2‑150 (MQ=255)
                                                                                                                                             agcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGt       <  1:1783578/150‑1 (MQ=255)
                                                                                                                                             agcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGt       <  2:903630/150‑1 (MQ=255)
                                                                                                                                             agcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGt       <  2:1293535/150‑1 (MQ=255)
                                                                                                                                             agcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGt       >  1:1973119/1‑150 (MQ=255)
                                                                                                                                             agcagcTGCTGCTACCACGGAGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGt       >  2:743117/1‑150 (MQ=255)
                                                                                                                                              gcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGtt      >  1:148184/1‑150 (MQ=255)
                                                                                                                                              gcagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGtt      >  2:44900/1‑150 (MQ=255)
                                                                                                                                               cagcTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTa     <  1:44900/150‑1 (MQ=255)
                                                                                                                                                 gctgctgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTAtt   <  2:1053418/150‑1 (MQ=255)
                                                                                                                                                 cctgctgctACCACGGCGACACAAAGCTGGGCGGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGGGGCCTCTGTGGAGGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTAtt   <  1:1748035/149‑1 (MQ=255)
                                                                                                                                                  ctgctgctACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTg  <  2:16521/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGAAGAGATTACCTACACTGCACGCGGTGTGGCTGCCGACCCGGACGACGTCCTATTTTCGGGCAATCGTGCAACGCAGCCAGAATATGACGAATACGATCCATTATTAAACGGTGCGCCAATTACCGAACCTGTCGCTGTAGCAGCTGCTGCTACCACGGCGACACAAAGCTGGGCTGCGCCGGTTGAACCTGTGACTCAGACGCCGCCTGTTGCCTCTGTTGATGTTCCACCTGCGCAACCTACAGTAGCCTGGCAGCCTGTACCGGGTCCACAAACGGGAGAGCCGGTTATTG  >  NC_000913/934048‑934343

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: