Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,400,330 C→G 5.5% L15V (CTG→GTG)  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,3300CG 5.5% 409.2 / 10.7 123L15V (CTG→GTG) hfqglobal sRNA chaperone; HF‑I, host factor for RNA phage Q beta replication
Reads supporting (aligned to +/- strand):  ref base C (54/62);  new base G (5/2);  total (59/64)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 125 total aligned reads displayed.

GTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCC  >  NC_000913/4400182‑4400478
                                                                                                                                                    |                                                                                                                                                     
         aTCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCggg                                                                                                                                             >  1:1578727/1‑149 (MQ=255)
          tCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGa                                                                                                                                            >  2:291791/1‑149 (MQ=255)
          tCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGa                                                                                                                                            >  2:514086/1‑149 (MQ=255)
               ggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                                                                                       >  1:641109/1‑149 (MQ=255)
               ggCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACgtg                                                                                                                                       >  1:135904/1‑149 (MQ=255)
                gCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCCGACGCAGTGCGTTCAGGAACggat                                                                                                                                      <  1:725680/149‑4 (MQ=255)
                   gAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTcc                                                                                                                                   <  2:591908/149‑1 (MQ=255)
                       gtgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGtt                                                                                                                               >  1:1362316/1‑149 (MQ=255)
                       gtgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGtt                                                                                                                               <  2:539341/149‑1 (MQ=255)
                        tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGttt                                                                                                                              >  2:1398316/1‑149 (MQ=255)
                        tgtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGttt                                                                                                                              <  2:582025/149‑1 (MQ=255)
                         gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                                                                                             <  1:390628/149‑1 (MQ=255)
                         gtACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTc                                                                                                                             <  2:185565/149‑1 (MQ=255)
                          tACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCt                                                                                                                            >  1:886408/1‑149 (MQ=255)
                           aCAATTGAGACGTATCGTGCGCGATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta                                                                                                                           <  1:101927/149‑1 (MQ=255)
                           aCAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta                                                                                                                           <  1:909749/149‑1 (MQ=255)
                           aCAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGAAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCta                                                                                                                           >  2:1153738/1‑149 (MQ=255)
                            cAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtat                                                                                                                          <  2:1578727/149‑1 (MQ=255)
                            cAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCtat                                                                                                                          <  2:1578726/149‑1 (MQ=255)
                                  agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTatttattt                                                                                                                    >  1:1539829/1‑149 (MQ=255)
                                  agaCGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTatttattt                                                                                                                    <  1:1420048/149‑1 (MQ=255)
                                     cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                                                                                                 <  1:605225/149‑1 (MQ=255)
                                     cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                                                                                                 <  2:299578/149‑1 (MQ=255)
                                     cGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAACATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGt                                                                                                                 <  2:65680/149‑1 (MQ=255)
                                       tATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGa                                                                                                               >  2:548989/1‑149 (MQ=255)
                                         tCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAt                                                                                                             <  1:107744/149‑1 (MQ=255)
                                         tCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTACTGAACGCACTCAGTCGGGAACGTTTTCCAGTTTCTATTTATTTGGTTAAt                                                                                                             >  1:639565/1‑149 (MQ=255)
                                          cGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATg                                                                                                            <  1:522532/149‑1 (MQ=255)
                                           gTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATgg                                                                                                           <  2:800510/149‑1 (MQ=255)
                                             gcgcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTa                                                                                                         <  1:461759/149‑1 (MQ=255)
                                               gcAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTAtt                                                                                                       >  2:677060/1‑149 (MQ=255)
                                                      tttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGc                                                                                                <  2:1105129/149‑1 (MQ=255)
                                                       ttCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCa                                                                                               >  1:544353/1‑149 (MQ=255)
                                                         cAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAg                                                                                             <  1:472165/149‑1 (MQ=255)
                                                           gAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                           >  1:747398/1‑149 (MQ=255)
                                                           gAATCGAAAGGTTCAAAGTAAAAATAAGCATATAAGGAAAAGAGAGTATGGCTAAGGGTCAATCTTTACATGATCCGTTCCTGTACGCTCTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAggg                                                                                           <  2:630417/149‑1 (MQ=255)
                                                            aaTCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGc                                                                                          >  2:1365903/1‑149 (MQ=255)
                                                              tCGAAAGGTTCATAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCAGTGCGTCGGGAACGTGTTCCAGTGTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCaa                                                                                        <  1:112799/149‑1 (MQ=255)
                                                                 aaaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTACGCTGCAAGGGCAAATc                                                                                     >  2:245568/1‑149 (MQ=255)
                                                                 aaaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc                                                                                     >  1:1094207/1‑149 (MQ=255)
                                                                 aaaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATc                                                                                     >  1:638899/1‑149 (MQ=255)
                                                                 aaaGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTCAATGGTATTAAGCTGCAAGGGCaaagc                                                                                     >  2:62460/1‑147 (MQ=255)
                                                                     gTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCATTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAAAGGTATTAAGCTGCAAGGGCAAATCGAGt                                                                                 >  1:359549/1‑149 (MQ=255)
                                                                      ttCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc                                                                                >  2:1291877/1‑149 (MQ=255)
                                                                      ttCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTc                                                                                >  2:763910/1‑149 (MQ=255)
                                                                       tCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCt                                                                               <  1:266045/149‑1 (MQ=255)
                                                                         aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAGGGGCAAATCGAGTCttt                                                                             <  2:1141913/149‑1 (MQ=255)
                                                                         aaaGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCttt                                                                             >  2:351556/1‑149 (MQ=255)
                                                                             tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                         >  2:659122/1‑149 (MQ=255)
                                                                             tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                         <  1:291791/149‑1 (MQ=255)
                                                                             tACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGAAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGAt                                                                         <  1:1153738/149‑1 (MQ=255)
                                                                              aCAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                        <  2:361758/149‑1 (MQ=255)
                                                                              aCAAATAAGCATATAAGGAAAAGAGAGAATGGCGAAGGGGCAATCTTTACAAGATCCGGTCCTGAACGCACTGCGTCGGGAACGTGGTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATc                                                                        <  2:273551/149‑1 (MQ=255)
                                                                                aaaTAAGCATATTAGGAAAAGAGAGAATGGCTTAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCat                                                                      >  2:773743/1‑148 (MQ=255)
                                                                                aaaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAg                                                                      <  1:388424/149‑1 (MQ=255)
                                                                                aaaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAg                                                                      >  1:458435/1‑149 (MQ=255)
                                                                                aaaTAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATACGTTCCTGAACGCACTGAGTCGGGAACGTGTTCCAGATTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAg                                                                      >  1:799880/1‑149 (MQ=255)
                                                                                    aaGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCg                                                                  <  1:514086/149‑1 (MQ=255)
                                                                                       cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGa                                                               <  1:1035198/149‑1 (MQ=255)
                                                                                       cATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGa                                                               >  2:485439/1‑149 (MQ=255)
                                                                                         tataAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATc                                                             >  1:726505/1‑149 (MQ=255)
                                                                                            aaGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTg                                                          >  1:1391023/1‑149 (MQ=255)
                                                                                             aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGt                                                         <  1:597722/149‑1 (MQ=255)
                                                                                             aGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGt                                                         >  2:880109/1‑149 (MQ=255)
                                                                                              gtaaaaGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                        <  1:524655/147‑1 (MQ=255)
                                                                                              ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCTTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTACAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                        <  1:911069/149‑1 (MQ=255)
                                                                                              ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGGACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                        <  2:667439/149‑1 (MQ=255)
                                                                                              ggAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGTGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGtt                                                        <  1:227459/149‑1 (MQ=255)
                                                                                               gAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGttt                                                       >  2:1180567/1‑148 (MQ=255)
                                                                                                    gagagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGATGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGaaaaa                                                  >  1:536012/1‑149 (MQ=255)
                                                                                                     agagaATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAac                                                 >  2:1594422/1‑149 (MQ=255)
                                                                                                         aaTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACgg                                             <  2:50333/149‑1 (MQ=255)
                                                                                                          aTGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGt                                            >  2:631009/1‑149 (MQ=255)
                                                                                                           tGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGa                                                                          <  1:196557/118‑1 (MQ=255)
                                                                                                            ggCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCa                                          >  1:1338853/1‑149 (MQ=255)
                                                                                                             gCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAg                                         <  2:726505/149‑1 (MQ=255)
                                                                                                                    gtcAATCTTTACAAGATCCGTTAATGAACGCACTTCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTTAATGGTATTAAGCTGCAAGGTCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                  <  2:909100/147‑1 (MQ=255)
                                                                                                                    ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                  >  1:601810/1‑149 (MQ=255)
                                                                                                                    ggCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATg                                  <  1:208049/149‑1 (MQ=255)
                                                                                                                     gCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATgg                                 >  2:400658/1‑149 (MQ=255)
                                                                                                                         tCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGTGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTa                             <  2:1084114/149‑1 (MQ=255)
                                                                                                                           tttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGACTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACa                           >  2:1238764/1‑149 (MQ=255)
                                                                                                                            ttACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATCGTTTACaa                          >  1:16722/1‑149 (MQ=255)
                                                                                                                               tggaacacGTTCCCGACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa                       >  2:1334467/8‑149 (MQ=38)
                                                                                                                               tggaacacGTTCCCGACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa                       >  2:350759/8‑149 (MQ=38)
                                                                                                                               cAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCa                       >  1:1545417/1‑149 (MQ=255)
                                                                                                                                aaGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAc                      <  1:659122/149‑1 (MQ=255)
                                                                                                                                aaGATCCGTTCCTGAACGCACTCCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGATCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAc                      <  2:63414/149‑1 (MQ=255)
                                                                                                                                   acacGTTCCCGACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCAcgcg                   >  1:1487427/4‑149 (MQ=255)
                                                                                                                                    tCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                  <  1:390772/149‑1 (MQ=255)
                                                                                                                                    cacGTTCCCGACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                  >  1:1083199/3‑149 (MQ=255)
                                                                                                                                    cacGTTCCCGACGCAGTGCGTTCAGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGa                  >  2:504403/3‑149 (MQ=255)
                                                                                                                                        ttCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTc              <  1:631009/149‑1 (MQ=255)
                                                                                                                                           cTGAACGCACTGCGTCGGTAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTACGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc           <  2:799880/149‑1 (MQ=255)
                                                                                                                                           cTGAACGCACTGCGTCGGGAACTTGTTCCAGTTTCTATTTATTTGGTGAATGGGATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTAc           <  2:1050005/149‑1 (MQ=255)
                                                                                                                                            tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt          <  1:245568/149‑1 (MQ=255)
                                                                                                                                            tGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACt          <  1:1487642/149‑1 (MQ=255)
                                                                                                                                             gAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACtg         <  2:1362316/149‑1 (MQ=255)
                                                                                                                                                   aCTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTcc   >  2:705640/1‑149 (MQ=255)
                                                                                                                                                    cTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTccc  >  2:137167/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTTGGTGCTATCGCAGGCTGAATGTGTACAATTGAGACGTATCGTGCGCAATTTTTTCAGAATCGAAAGGTTCAAAGTACAAATAAGCATATAAGGAAAAGAGAGAATGGCTAAGGGGCAATCTTTACAAGATCCGTTCCTGAACGCACTGCGTCGGGAACGTGTTCCAGTTTCTATTTATTTGGTGAATGGTATTAAGCTGCAAGGGCAAATCGAGTCTTTTGATCAGTTCGTGATCCTGTTGAAAAACACGGTCAGCCAGATGGTTTACAAGCACGCGATTTCTACTGTTGTCCC  >  NC_000913/4400182‑4400478

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