Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,882,476 A→G 5.7% T418T (ACT→ACC)  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,4760AG 5.7% 327.9 / 7.1 106T418T (ACT→ACC) dnaAchromosomal replication initiator protein DnaA, DNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base A (46/53);  new base G (3/3);  total (49/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

TGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTC  >  NC_000913/3882332‑3882619
                                                                                                                                                |                                                                                                                                               
tGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGt                                                                                                                                               >  2:32155/1‑147 (MQ=255)
tGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTc                                                                                                                                              <  5:402241/148‑1 (MQ=255)
 gACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGt                                                                                                                                               >  2:319534/1‑146 (MQ=255)
        ttCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGACAATCTCCGGCAGACTGTGGTTAGTAAGCTCttt                                                                                                                                      <  6:330157/148‑1 (MQ=255)
         tCTGATTAAATTTGAAAAATCTTCATTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTc                                                                                                                                     >  1:323405/1‑148 (MQ=255)
           tGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGc                                                                                                                                   >  7:316201/1‑148 (MQ=255)
            gATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACTCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCga                                                                                                                                   <  6:126237/147‑2 (MQ=255)
             aTTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCg                                                                                                                                    <  6:397697/145‑1 (MQ=255)
             aTTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCg                                                                                                                                    <  7:141559/145‑1 (MQ=255)
             aTTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGc                                                                                                                                   >  6:1756/1‑146 (MQ=255)
              ttAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAg                                                                                                                                >  3:416099/1‑148 (MQ=255)
               tAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTGGTCAGCTCTTTCGCCAgc                                                                                                                               >  4:325590/1‑148 (MQ=255)
               tAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCAGGCAGACTGTGGTTGGTCAGCTCTTTCGCCAgc                                                                                                                               >  4:319715/1‑148 (MQ=255)
               tAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCAAATCTCCGGCAGACTGTGGTTGGTCAGCTCTTTCGCCAgc                                                                                                                               >  4:319710/1‑148 (MQ=255)
                   tttGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCCCCGTCGTGTGGTCACGGCCGCCAAACGCATCGCCAGTCTCCGGCAGACGGGGGTGAGTCTGCTCTTTCGCCGGcgccg                                                                                                                           >  2:422185/1‑147 (MQ=255)
                     tGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccat                                                                                                                          <  5:216448/147‑1 (MQ=255)
                        aaaaTCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCATACGGCAGGCATGAAGCACCGTAGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGccccatcgc                                                                                                                       >  1:424107/1‑147 (MQ=255)
                         aaaTCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccatcgcca                                                                                                                     <  8:411177/148‑1 (MQ=255)
                          aaTCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGcgccatcgccat                                                                                                                    >  3:396589/1‑148 (MQ=255)
                                cttTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCt                                                                                                                  <  8:133585/144‑1 (MQ=255)
                                   tGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGgcg                                                                                                              <  2:467760/145‑1 (MQ=255)
                                   tGATATCGTGGCTCTCTCCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGGCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGcgcg                                                                                                            <  8:414557/147‑1 (MQ=255)
                                        tCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGAcgcg                                                                                                      >  8:153097/1‑148 (MQ=255)
                                         cGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCgg                                                                                                           <  5:266881/142‑1 (MQ=255)
                                         cGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCgg                                                                                                           >  6:266881/1‑142 (MQ=255)
                                           tGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCa                                                                                                   >  4:287145/1‑148 (MQ=255)
                                            ggCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCCTCTGGCGCGGACGCGCCAc                                                                                                  >  6:195810/1‑148 (MQ=255)
                                              ctctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGTCCGCGCCagat                                                                                                <  5:1756/148‑4 (MQ=255)
                                                ctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGac                                                                                              >  6:80937/1‑147 (MQ=255)
                                                ctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAg                                                                                              >  3:45269/1‑148 (MQ=255)
                                                ctctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACAGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTGCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCAGGCGCGGACGCGCCACCGAg                                                                                              >  8:111435/1‑148 (MQ=255)
                                                  ctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCg                                                                                            >  1:45735/1‑148 (MQ=255)
                                                  ctTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCg                                                                                            >  8:180243/1‑148 (MQ=255)
                                                       cGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGa                                                                                                          >  8:453891/1‑129 (MQ=255)
                                                       cGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGa                                                                                                          <  7:453891/129‑1 (MQ=255)
                                                        gCAACTGCGCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCg                                                                                      <  1:422185/148‑1 (MQ=255)
                                                         cAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGa                                                                                     >  3:333852/1‑148 (MQ=255)
                                                               cTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCt                                                                                 <  4:416099/146‑1 (MQ=255)
                                                                 cGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTgg                                                                              >  6:60677/1‑147 (MQ=255)
                                                                   aTCTTACGGCAGGCATGAAGCACCGTCGTGTCGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATAGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGa                                                                             >  1:123268/1‑146 (MQ=255)
                                                                        aCGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaa                                                                            <  1:90473/142‑1 (MQ=255)
                                                                        aCGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaa                                                                            <  1:261715/142‑1 (MQ=255)
                                                                        aCGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAg                                                                          <  1:268267/144‑1 (MQ=255)
                                                                        aCGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGga                                                                        <  1:32155/146‑1 (MQ=255)
                                                                         cGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGAt                                                                     <  8:316091/148‑1 (MQ=255)
                                                                           gcaggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                           >  8:388949/1‑140 (MQ=255)
                                                                           gcaggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGaaa                                                                           <  7:388949/140‑1 (MQ=255)
                                                                            caggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCcg                                                                  <  5:60677/148‑1 (MQ=255)
                                                                             aggcaTGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCcgc                                                                 >  3:365099/1‑148 (MQ=255)
                                                                                  tGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGAGCGACGCTTGGAAAGGAGATCCGCGACtt                                                            <  6:404351/148‑1 (MQ=255)
                                                                                   gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACttt                                                           >  4:128097/1‑148 (MQ=255)
                                                                                   gAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACttt                                                           >  1:262988/1‑148 (MQ=255)
                                                                                       cACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                       <  2:262988/148‑1 (MQ=255)
                                                                                         ccGTCTTGTGTTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCtt                                                     <  4:330075/148‑1 (MQ=255)
                                                                                         ccGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGTCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGCTCCGCGACTTTGATCtt                                                     >  3:149836/1‑148 (MQ=255)
                                                                                         ccGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCtt                                                     <  2:323405/148‑1 (MQ=255)
                                                                                         ccGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCtt                                                     <  4:333852/148‑1 (MQ=255)
                                                                                          cgtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTg                                                    >  2:459684/1‑148 (MQ=255)
                                                                                           gtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                   <  2:45735/148‑1 (MQ=255)
                                                                                           gtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                   <  4:396589/148‑1 (MQ=255)
                                                                                           gtcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTg                                                    >  5:169311/1‑147 (MQ=255)
                                                                                           ggcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                   <  2:87182/146‑1 (MQ=255)
                                                                                           ggcgggtgGTCACGGCCACCAACCGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                   <  1:90661/143‑1 (MQ=255)
                                                                                            tcgtGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATc                                                       >  1:214222/1‑143 (MQ=255)
                                                                                              gtgtGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtagt                                                <  7:48757/148‑1 (MQ=255)
                                                                                              gtgtGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtagt                                                <  7:180243/148‑1 (MQ=255)
                                                                                                 tgGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTAc                                              <  7:153097/147‑1 (MQ=255)
                                                                                                  ggTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtagt                                                <  4:169681/144‑1 (MQ=255)
                                                                                                  ggTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgtagt                                                >  3:169681/1‑144 (MQ=255)
                                                                                                         gCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGcc                                        >  5:97072/1‑145 (MQ=255)
                                                                                                         gCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGcc                                        <  6:97072/145‑1 (MQ=255)
                                                                                                           caccaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGt                                   <  1:76598/148‑1 (MQ=255)
                                                                                                            accaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGCGCGGGATCGAGGCTTGGAAAGGAGATCCGCGCCTTTGATCTTGTAGTACTCCGCCGCCGtc                                  >  4:384390/1‑148 (MQ=255)
                                                                                                            accaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGga                                                                        >  4:154069/1‑110 (MQ=255)
                                                                                                            accaAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGtc                                  <  7:21351/148‑1 (MQ=255)
                                                                                                                      tCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATAGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAatat                        <  2:123268/148‑1 (MQ=255)
                                                                                                                       cGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATAtt                       >  8:237764/1‑148 (MQ=255)
                                                                                                                         ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCTTCTGGCGCGGACGCGCCACCGAGCGGGAGCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                     <  3:411490/148‑1 (MQ=255)
                                                                                                                         ccAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGt                     <  3:156098/148‑1 (MQ=255)
                                                                                                                           aaTCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTc                    >  3:132464/1‑147 (MQ=255)
                                                                                                                            aTCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGa                  <  1:459684/148‑1 (MQ=255)
                                                                                                                             tctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                 >  3:356414/1‑148 (MQ=255)
                                                                                                                              ctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGAt                                                                     >  2:381084/1‑95 (MQ=255)
                                                                                                                              ctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGAt                                                                     <  1:381084/95‑1 (MQ=255)
                                                                                                                              ctcCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATg                >  2:341870/1‑148 (MQ=255)
                                                                                                                                 gggCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCTCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGt              <  3:326292/146‑1 (MQ=255)
                                                                                                                                 cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                 <  3:287145/144‑1 (MQ=255)
                                                                                                                                 cGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGAt                 <  3:128097/144‑1 (MQ=255)
                                                                                                                                  ggCAGTCTGTGGTTGGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGTCGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa            <  3:319710/148‑1 (MQ=255)
                                                                                                                                  ggCAGACTGTGGTTGGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa            <  3:325590/148‑1 (MQ=255)
                                                                                                                                  ggCAGACTCTGGTAGGTCAGCTCTTGCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGa            <  3:319715/148‑1 (MQ=255)
                                                                                                                                   gCAGACTGTGGTTAGTCAGCTATTTCGCCACCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACCCTTGGAAAGCAGATCCCAGACTTTGATATTGTAGTACTACGCCACCGTCTTCTGAATATTGTCGAAGGTGa            >  6:234706/1‑147 (MQ=255)
                                                                                                                                    cAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                   >  2:322391/1‑107 (MQ=255)
                                                                                                                                    cAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTgt                                                   <  1:322391/107‑1 (MQ=255)
                                                                                                                                       aCTGTGGATAGTCCGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGCCGCTTGGAAAGGAGATCCGCAACTTTGATCTTGTAGTACTCCGCCACCTTCTTCTGAATATTGTCGATGGTGCCCatt       >  2:381905/1‑146 (MQ=255)
                                                                                                                                           tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGttttt   >  2:232000/1‑148 (MQ=255)
                                                                                                                                           tgGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtttat   <  4:356414/148‑3 (MQ=255)
                                                                                                                                             gTTTGTCTGCTCTTTGGCCAGCGCCATCGCCATTTGGCGAGGACGCGCCACCGAGGGGGTTCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATTGTGACCAGTTTTTc  <  5:234706/147‑1 (MQ=255)
                                                                                                                                              ttAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtt      >  8:31365/1‑142 (MQ=255)
                                                                                                                                              ttAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGtt      <  7:31365/142‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TGACAATGTTCTGATTAAATTTGAAAAATCTTCTTTGATATCGTGGCTCTCTTCACGCAACTGCTCGATCTTACGGCAGGCATGAAGCACCGTCGTGTGGTCACGGCCACCAAACGCATCGCCAATCTCCGGCAGACTGTGGTTAGTCAGCTCTTTCGCCAGCGCCATCGCCATCTGGCGCGGACGCGCCACCGAGCGGGATCGACGCTTGGAAAGGAGATCCGCGACTTTGATCTTGTAGTACTCCGCCACCGTCTTCTGAATATTGTCGATGGTGACCAGTTTTTC  >  NC_000913/3882332‑3882619

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