Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,400,313 A→C 100% D9A (GAT→GCT)  hfq → global sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,400,3130AC99.3% 588.6 / ‑4.2 170D9A (GAT→GCT) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base C (82/84);  total (84/86)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

GACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATC  >  NC_000913/4400217‑4400408
                                                                                                |                                                                                               
  cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGt                                                                                            >  1:1183849/1‑100 (MQ=255)
    tATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTc                                                                                          >  2:479281/1‑100 (MQ=255)
       cGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTg                                                                                       <  1:2010898/100‑1 (MQ=255)
       cGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTg                                                                                       <  1:34734/100‑1 (MQ=255)
        gTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGa                                                                                      >  2:1273108/1‑100 (MQ=255)
          gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAAc                                                                                    <  2:1385441/100‑1 (MQ=255)
          gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAAc                                                                                    >  1:3002737/1‑100 (MQ=255)
          gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAAc                                                                                    <  1:3072620/100‑1 (MQ=255)
           cgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACg                                                                                   >  2:3177070/1‑100 (MQ=255)
             cAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCa                                                                                 >  2:350562/1‑100 (MQ=255)
               aTTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACt                                                                               >  2:2677051/1‑100 (MQ=255)
                ttttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGTAATCTTTACAAGCTCCGTTCCTGAACGCACTg                                                                              <  1:3631372/100‑1 (MQ=255)
                ttttttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTg                                                                              >  2:1990190/1‑100 (MQ=255)
                     tCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCg                                                                         <  1:3326183/100‑1 (MQ=255)
                       aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCggg                                                                       <  2:1039878/100‑1 (MQ=255)
                       aGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCgg                                                                        >  1:3139727/1‑99 (MQ=255)
                        gAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGa                                                                      >  2:181823/1‑100 (MQ=255)
                        gAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGa                                                                      >  1:3099629/1‑100 (MQ=255)
                          aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCCGACGCAGTGCGttcaggaac                                                                   <  1:1707048/100‑9 (MQ=37)
                          aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCCGACGCAGTGCGttcaggaac                                                                   <  2:2339013/100‑9 (MQ=37)
                          aTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCCGACGCAGTGCGttcaggaac                                                                   <  2:3294291/100‑9 (MQ=37)
                          aTCGAAAGGCTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTCAGTGGCAAACTTTACAAGCTCCGTTCCTGAACGCCCTACGTCGGGAAc                                                                    >  1:1022307/1‑100 (MQ=255)
                           tCGACAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGTCAATCTATACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACg                                                                   >  1:517964/1‑100 (MQ=255)
                            cGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgt                                                                  >  2:1591310/1‑100 (MQ=255)
                             gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 >  1:2142618/1‑100 (MQ=255)
                             gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 >  2:1512329/1‑100 (MQ=255)
                             gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 >  2:172136/1‑100 (MQ=255)
                             gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 >  2:185211/1‑100 (MQ=255)
                             gAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                 <  2:3471566/100‑1 (MQ=255)
                                 tcttCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             <  1:983944/98‑1 (MQ=255)
                                 ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             <  1:37990/100‑1 (MQ=255)
                                 ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             <  2:1660817/100‑1 (MQ=255)
                                 ggTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                             <  2:2025808/100‑1 (MQ=255)
                                  gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCa                                                            >  1:671158/1‑100 (MQ=255)
                                  gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCa                                                            >  1:2203432/1‑100 (MQ=255)
                                    tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAga                                                          >  2:983952/1‑99 (MQ=255)
                                      aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGttt                                                        >  2:1261987/1‑100 (MQ=255)
                                       aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                       <  1:1837864/100‑1 (MQ=255)
                                       aaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                       <  1:185209/100‑1 (MQ=255)
                                       aaGTACAAATAAGCATATAAGGAAAAGAGAGAAGGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAATGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                       >  1:2327124/1‑100 (MQ=255)
                                        aGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCt                                                      >  2:3433655/1‑100 (MQ=255)
                                         gTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta                                                     <  1:3604158/100‑1 (MQ=255)
                                         gTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta                                                     <  1:2888985/100‑1 (MQ=255)
                                           aCAAATAAGCATATAAGGAAAAGATAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtat                                                    <  2:646901/99‑1 (MQ=255)
                                           aCAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtat                                                    <  1:1018989/99‑1 (MQ=255)
                                           aCAAATAAGCATACAAGGAAAAGAGAGCATCGCTAACGGGCAATCATTCCAAGCTCCGTTCCTGAATGCACTGCGTCGGGAACGTGTGCCAGTTTCtatt                                                   <  1:1506455/100‑1 (MQ=255)
                                            cAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtattt                                                  >  2:3362212/1‑100 (MQ=255)
                                               aTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtatttatt                                               <  1:1105160/100‑1 (MQ=255)
                                                 aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTg                                             <  2:1803609/100‑1 (MQ=255)
                                                 aaGAATATAAGGAAAAGAGAGAATGACTAAGGGGCAATCGTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTg                                             <  1:753325/100‑1 (MQ=255)
                                                  aGCATATAAGGAAAAGAGAGAATGGCTAAGGTGCAATCTTTACAAGCTCCTTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTgg                                            <  2:1787060/100‑1 (MQ=255)
                                                   gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                           >  1:621356/1‑100 (MQ=255)
                                                   gCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                           <  1:1949547/100‑1 (MQ=255)
                                                    cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTg                                          <  2:2485498/100‑1 (MQ=255)
                                                    cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCGATTTAGTTGGTg                                          >  1:1203705/1‑100 (MQ=255)
                                                     atatAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGa                                         <  2:2970465/100‑1 (MQ=255)
                                                      tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGaa                                        >  2:1368568/1‑100 (MQ=255)
                                                       ataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                       <  2:624658/100‑1 (MQ=255)
                                                        taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                      >  1:1717328/1‑100 (MQ=255)
                                                        taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                      >  2:3660627/1‑100 (MQ=255)
                                                        taAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                      >  2:2145004/1‑100 (MQ=255)
                                                         aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                      <  1:202779/99‑1 (MQ=255)
                                                          aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGt                                    >  1:818623/1‑100 (MQ=255)
                                                           ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                   >  1:1570078/1‑100 (MQ=255)
                                                           ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                   >  1:338461/1‑100 (MQ=255)
                                                           ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                   <  2:1212176/100‑1 (MQ=255)
                                                            gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGACCGCACTGCGTCGGGAACGGGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                  >  2:1306007/1‑100 (MQ=255)
                                                            gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAt                                  <  1:2870/100‑1 (MQ=255)
                                                            gAAAAGAGAGAATGGCTAAGGGACAATCTTTACAAGCTACGTTCCTGAACGCACTGCGACGGAAACGTGTTCCCGTTTCTATTTATTTGGTTACTGGTAt                                  >  1:2579749/1‑100 (MQ=255)
                                                             aaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAtt                                 <  1:582359/100‑1 (MQ=255)
                                                               aaGAGAGAATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTan                               <  2:1076174/100‑2 (MQ=255)
                                                                agagagAATGCCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAg                              <  2:3541755/100‑1 (MQ=255)
                                                                 gagagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGc                             <  2:3650199/100‑1 (MQ=255)
                                                                  agagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCt                            <  2:2743569/100‑1 (MQ=255)
                                                                  agagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCt                            <  1:1779502/100‑1 (MQ=255)
                                                                   gagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                           >  1:687357/1‑100 (MQ=255)
                                                                   gagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                           <  1:86499/100‑1 (MQ=255)
                                                                   gagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                           <  2:3301316/100‑1 (MQ=255)
                                                                   gagaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTg                           >  2:3038932/1‑100 (MQ=255)
                                                                     gaATGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                         >  2:926232/1‑100 (MQ=255)
                                                                        tGGCTAAGGGGCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAgg                      >  2:1077390/1‑100 (MQ=255)
                                                                               ggggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc               <  2:2064620/100‑1 (MQ=255)
                                                                                gggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCg              <  1:1743610/100‑1 (MQ=255)
                                                                                gggCAATCTTTACAAGCTCCGCTCCTGAACGCACTGCGTCGGGAACGCGTTCAAGTGTCTATGTATTTGGTGAGTGGTAGTAATCTGCAAGGGCTAATc               >  1:1285637/1‑99 (MQ=255)
                                                                                 ggaaaTTTTTATAGGTTCCGTTCCTGAACGTATTGTGTCGGGAACGTGTTCCAGTTTTTATTTATTTTGTGAATGGTATTAAGCTGTAAGGGCAAATc               <  1:466428/95‑1 (MQ=255)
                                                                                 ggCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGa             >  1:2921002/1‑100 (MQ=255)
                                                                                 cgCAATCTTTACAAGCTCAGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGACTGGTATTAAGCTGCAAGGGCAAATCGa             >  1:383504/2‑100 (MQ=255)
                                                                                  gCAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAg            >  2:1470983/1‑100 (MQ=255)
                                                                                   cAATCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGt           >  1:1124570/1‑100 (MQ=255)
                                                                                    aagcTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGGTTCTAGTTCTTTGGTATATGGTATTAAGCGGCAAGTGCAGATCGATTc          >  1:2717352/4‑100 (MQ=255)
                                                                                    aaTCTTTACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc          >  1:619804/1‑100 (MQ=255)
                                                                                        tttACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCtttt      >  2:2152143/1‑100 (MQ=255)
                                                                                         ttACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTg     >  2:2847213/1‑100 (MQ=255)
                                                                                          tACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTg     <  1:520444/99‑1 (MQ=255)
                                                                                          tACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa    <  1:2587422/100‑1 (MQ=255)
                                                                                          tACAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa    >  1:1326099/1‑100 (MQ=255)
                                                                                           aCAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa    >  1:3723236/1‑99 (MQ=255)
                                                                                           aCAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt   <  1:737433/100‑1 (MQ=255)
                                                                                            cAAGCTCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc  <  1:3715340/100‑1 (MQ=255)
                                                                                            cAAGCTCCGTTCCTGAACGCACTGCGTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc  <  1:1990175/100‑1 (MQ=255)
                                                                                                |                                                                                               
GACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATC  >  NC_000913/4400217‑4400408

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: