Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,882,514 G→T 5.7% P406T (CCG→ACG)  dnaA ← chromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,882,5140GT 5.7% 547.1 / 17.4 157P406T (CCG→ACG) dnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (75/73);  new base T (4/5);  total (79/78)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

Only 100 of 157 total aligned reads displayed.

ATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAG  >  NC_000913/3882369‑3882659
                                                                                                                                                 |                                                                                                                                                 
 tatCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGAcgc                                                                                                                                               <  3:140652/149‑1 (MQ=255)
  atCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGAcgcg                                                                                                                                              >  1:855286/1‑149 (MQ=255)
   tCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACgcgc                                                                                                                                             >  3:868474/1‑149 (MQ=255)
       ggCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCaga                                                                                                                                         <  3:1123445/149‑3 (MQ=255)
         ctctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagatg                                                                                                                                       <  1:1090486/149‑5 (MQ=255)
         ctctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGa                                                                                                                                       >  1:1187683/1‑149 (MQ=255)
            tctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCg                                                                                                                                    >  3:396764/1‑149 (MQ=255)
              ttCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagatgg                                                                                                                                      >  2:454858/1‑140 (MQ=255)
              ttCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCg                                                                                                                                    >  3:536877/1‑147 (MQ=255)
                cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGAt                                                                                                                                <  2:1202890/149‑1 (MQ=255)
                cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGAt                                                                                                                                <  3:154716/149‑1 (MQ=255)
                cACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGAt                                                                                                                                <  4:868474/149‑1 (MQ=255)
                 aCGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                               <  1:24299/149‑1 (MQ=255)
                 aCGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATc                                                                                                                               >  4:605811/1‑149 (MQ=255)
                        tGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCtt                                                                                                                        <  3:779428/149‑1 (MQ=255)
                        tGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCtt                                                                                                                        <  1:762916/149‑1 (MQ=255)
                           tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCgg                                                                                                                                   >  2:708999/1‑135 (MQ=255)
                           tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                     >  4:814181/1‑149 (MQ=255)
                           tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                     >  1:16079/1‑149 (MQ=255)
                           tCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                                                                     >  4:1453587/1‑149 (MQ=255)
                             gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   >  2:768820/1‑149 (MQ=255)
                             gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   >  1:525338/1‑149 (MQ=255)
                             gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   >  3:1380948/1‑149 (MQ=255)
                             gATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                                                                   <  4:459717/149‑1 (MQ=255)
                              aTCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAg                                                                                                                  <  2:38544/149‑1 (MQ=255)
                              aTCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAg                                                                                                                  >  2:224383/1‑149 (MQ=255)
                               tCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAgg                                                                                                                 >  4:1172068/1‑149 (MQ=255)
                               tCTTACGGCAGGCATGAAGCACCGTCGTGTGGTAACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAgg                                                                                                                 >  3:924384/1‑149 (MQ=255)
                                cTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGga                                                                                                                <  2:49307/149‑1 (MQ=255)
                                cTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGga                                                                                                                <  1:926453/149‑1 (MQ=255)
                                 ttACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGgag                                                                                                               >  4:491023/1‑149 (MQ=255)
                                  tACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGgaga                                                                                                              <  4:886739/149‑1 (MQ=255)
                                         ggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGa                                                                                                       <  4:599323/149‑1 (MQ=255)
                                              gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagat                                                                                                                                        <  3:926780/111‑4 (MQ=255)
                                              gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTg                                                                                                  >  1:594846/1‑149 (MQ=255)
                                              gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTg                                                                                                  >  1:1236916/1‑149 (MQ=255)
                                                aGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGAt                                                                                                >  1:250047/1‑149 (MQ=255)
                                                aGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGAt                                                                                                >  2:132994/1‑149 (MQ=255)
                                                 gCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                                                               <  2:1104769/149‑1 (MQ=255)
                                                 gCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                                                               >  3:797441/1‑149 (MQ=255)
                                                      gtcgtGTGGTCACGGCCACCAATCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGAACGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                                          <  4:1202741/149‑1 (MQ=255)
                                                      gtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                                          <  4:313733/149‑1 (MQ=255)
                                                      gtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                                          <  4:452612/149‑1 (MQ=255)
                                                       tcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtag                                                                                         >  2:1126151/1‑149 (MQ=255)
                                                         gtgtGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtagta                                                                                       >  4:7200/1‑149 (MQ=255)
                                                          tgtgGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTAc                                                                                      >  1:50483/1‑149 (MQ=255)
                                                             ggTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTcc                                                                                   <  4:185754/149‑1 (MQ=255)
                                                               tCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                                 >  3:629571/1‑149 (MQ=255)
                                                               tCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGc                                                                                 >  1:102960/1‑149 (MQ=255)
                                                                cACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGcc                                                                                <  3:491023/149‑1 (MQ=255)
                                                                   ggCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCAcc                                                                             >  1:698145/1‑149 (MQ=255)
                                                                         caAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                                       <  3:171417/149‑1 (MQ=255)
                                                                          aaaCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGTGGGATCGACGCTTGGAAAGGAGATCCGCGACCTTGATCTTGTAGTACTCCGCCACCGtcttct                                                                      <  4:396764/149‑1 (MQ=255)
                                                                          aaaCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttct                                                                      <  4:1029363/149‑1 (MQ=255)
                                                                               cATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAata                                                                 <  3:686069/149‑1 (MQ=255)
                                                                               cATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAata                                                                 <  3:814181/149‑1 (MQ=255)
                                                                                aTCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                                                                >  2:708992/1‑149 (MQ=255)
                                                                                aTCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAAATCCGCGGCTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                                                                >  2:1262868/1‑149 (MQ=255)
                                                                                   gCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                                                             >  4:535741/1‑149 (MQ=255)
                                                                                   gCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                                                             <  3:743999/149‑1 (MQ=255)
                                                                                    ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                                                            <  2:250047/149‑1 (MQ=255)
                                                                                    ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                                                            <  2:1236916/149‑1 (MQ=255)
                                                                                      aaTCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgta                                                                                          >  2:272475/1‑117 (MQ=255)
                                                                                          tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt                                                      >  3:1064956/1‑149 (MQ=255)
                                                                                          tcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt                                                      <  1:946937/149‑1 (MQ=255)
                                                                                            cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa                                                    >  3:523099/1‑149 (MQ=255)
                                                                                             ggCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGAc                                                   >  2:1396651/1‑149 (MQ=255)
                                                                                                    tgtgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttt                                            <  2:168398/149‑1 (MQ=255)
                                                                                                    tgtgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttt                                            <  4:1064956/149‑1 (MQ=255)
                                                                                                      tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTc                                          <  1:15187/149‑1 (MQ=255)
                                                                                                      tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCGTCTGAATATTGTCGATGGTGACCAGTTTTTc                                          >  4:156026/1‑149 (MQ=255)
                                                                                                        gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCg                                                                                                                                    >  4:771992/1‑57 (MQ=255)
                                                                                                        gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCg                                                                                                                                    <  3:771992/57‑1 (MQ=255)
                                                                                                        gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                                                                <  4:1400251/125‑1 (MQ=255)
                                                                                                        gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCt                                        >  3:370040/1‑149 (MQ=255)
                                                                                                        gTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCt                                        >  4:1046961/1‑149 (MQ=255)
                                                                                                             tCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAAt                                   <  2:16079/149‑1 (MQ=255)
                                                                                                               aGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                                       <  2:852524/111‑1 (MQ=255)
                                                                                                               aGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtcttc                                                                       >  1:852524/1‑111 (MQ=255)
                                                                                                                gCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGcc                                >  1:1035167/1‑149 (MQ=255)
                                                                                                                 ctctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCa                               >  1:940015/1‑149 (MQ=255)
                                                                                                                  tctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAg                              >  2:771937/1‑149 (MQ=255)
                                                                                                                  tctTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAg                              >  2:1308762/1‑149 (MQ=255)
                                                                                                                    tttCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCa                            >  2:1476200/1‑149 (MQ=255)
                                                                                                                          cAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTcgc                      >  1:362111/1‑149 (MQ=255)
                                                                                                                                 atcgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGCCTTTGATCTTGTAGTACTcc                                                                                   <  2:953084/81‑1 (MQ=255)
                                                                                                                                  tcgccatcTGGCGCGGACGCGCCCCCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGcc              >  2:948925/1‑149 (MQ=255)
                                                                                                                                   cgccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCt             <  2:855286/149‑1 (MQ=255)
                                                                                                                                    gccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTc            <  1:1053775/149‑1 (MQ=255)
                                                                                                                                     ccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCa                               >  2:189857/1‑129 (MQ=255)
                                                                                                                                     ccatcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCa                               <  1:189857/129‑1 (MQ=255)
                                                                                                                                      tcggtGGCGCGTCCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcac          >  1:287099/5‑149 (MQ=255)
                                                                                                                                       atcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                                       <  3:1331947/103‑1 (MQ=255)
                                                                                                                                       atcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATgg                                                       >  4:1331947/1‑103 (MQ=255)
                                                                                                                                       atcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCGATGCCAGCAAGTCGCGCAGCGCCTcacg         >  2:41038/1‑149 (MQ=255)
                                                                                                                                       atcTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcacg         <  4:452730/149‑1 (MQ=255)
                                                                                                                                         cTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTcacgca       <  1:69843/149‑1 (MQ=255)
                                                                                                                                             cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGaa   <  3:459227/149‑1 (MQ=255)
                                                                                                                                             cgcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGaa   <  2:931765/149‑1 (MQ=255)
                                                                                                                                              gcgGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAg  >  1:714280/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
ATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTCCTGCAATGCCAGCAAGTCGCGCAGCGCCTCACGCACGAAG  >  NC_000913/3882369‑3882659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: