Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,659,407 T→G 6.7% S21R (AGC→CGC)  hscB ← HscA co‑chaperone, J domain‑containing protein Hsc56; IscU‑specific chaperone HscAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,659,4070TG 6.7% 246.7 / 5.4 83S21R (AGC→CGC) hscBHscA co‑chaperone, J domain‑containing protein Hsc56; IscU‑specific chaperone HscAB
Reads supporting (aligned to +/- strand):  ref base T (40/37);  new base G (3/3);  total (43/40)
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

TTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTT  >  NC_000913/2659261‑2659549
                                                                                                                                                  |                                                                                                                                              
ttCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCa                                                                                                                                              <  2:36828/149‑1 (MQ=255)
   cgcgcgCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGGTAATGGTTGCAGATTGCGGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAgcg                                                                                                                                           <  2:421801/149‑1 (MQ=255)
     cgcgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGcc                                                                                                                                         >  2:837515/1‑149 (MQ=255)
      gcgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCt                                                                                                                                        <  1:1280158/149‑1 (MQ=255)
       cgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTg                                                                                                                                       <  2:1017019/149‑1 (MQ=255)
       cgcATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTg                                                                                                                                       <  2:299579/149‑1 (MQ=255)
       cgcATTAACGGAGGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTg                                                                                                                                       <  2:304057/149‑1 (MQ=255)
          aTTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGt                                                                                                                                    >  1:608885/1‑149 (MQ=255)
             aaCGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATc                                                                                                                                 >  2:1325062/1‑149 (MQ=255)
             aaCGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATc                                                                                                                                 >  2:118713/1‑149 (MQ=255)
             aaCGGATGACGCAGCGTTTGCCAGGCATGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATc                                                                                                                                 <  1:255165/149‑1 (MQ=255)
              aCGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCg                                                                                                                                >  2:1590886/1‑149 (MQ=255)
              aCGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGACTGCTGTACGGCGGCGAGTTGGTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCg                                                                                                                                <  2:116684/149‑1 (MQ=255)
               cGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGa                                                                                                                               <  2:1029313/149‑1 (MQ=255)
               cGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGa                                                                                                                               <  2:911067/149‑1 (MQ=255)
                ggATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAg                                                                                                                              <  2:725552/149‑1 (MQ=255)
                ggATGACGCAGCGTTAGCCAGTCCTCGTTAATGGCTGCAGCTTGCGGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTTGCGAATTTATCAGCATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAg                                                                                                                              <  2:233934/149‑1 (MQ=255)
                 gATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGt                                                                                                                             >  1:511819/1‑149 (MQ=255)
                   tGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTg                                                                                                                           <  1:199061/149‑1 (MQ=255)
                    gACGCAGCGTTTGCCAGGCCTGGTTAATGGTGGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGa                                                                                                                          <  2:316411/149‑1 (MQ=255)
                     aCGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGAt                                                                                                                         <  2:575288/149‑1 (MQ=255)
                     aCGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGAt                                                                                                                         >  2:525005/1‑149 (MQ=255)
                          gCGTTTGCCTGGCCTGTTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGTCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGATAACGCAGGCGCAGCGCCTGTGTATCGAGTTGATAGCgt                                                                                                                    <  2:574148/149‑2 (MQ=255)
                           cGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCggg                                                                                                                   <  1:656673/149‑1 (MQ=255)
                               tgccctgCCTGGTTAATGGTTGCAGAATGCTGGACGGCGGCGAGTTGTGCCGCCTGGCTTCCGCTGGGGAATTTATCAGGAAGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcagg                                                                                                               <  1:259653/143‑1 (MQ=255)
                                gCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGggcaggc                                                                                                              >  1:388711/1‑149 (MQ=255)
                                  cAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCaa                                                                                                            <  1:447743/149‑1 (MQ=255)
                                  cAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCACTGGCGAATTTATCAGGATGAAACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTAACGAGTTGATAGCGGGCAGGCaa                                                                                                            >  2:685449/1‑149 (MQ=255)
                                   aGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTAGATAGCGGGCAGGCAAg                                                                                                           >  2:173849/1‑149 (MQ=255)
                                    ggCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGc                                                                                                          <  2:1140872/149‑1 (MQ=255)
                                    ggCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGc                                                                                                          <  2:1505901/149‑1 (MQ=255)
                                           ttAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAACCAAagag                                                                                                   >  2:594890/1‑149 (MQ=255)
                                           ttAATGGGAGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAagag                                                                                                   <  2:1417073/149‑1 (MQ=255)
                                            tAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAgg                                                                                                  >  2:383415/1‑149 (MQ=255)
                                            tAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAATTGATAGCGGGCAGGCAAGCCAAAGAgg                                                                                                  >  2:646661/1‑149 (MQ=255)
                                                       gATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATc                                                                                       <  2:160249/149‑1 (MQ=255)
                                                          tGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCAt                                                                                    <  1:1075761/149‑1 (MQ=255)
                                                               tACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACAt                                                                               >  2:307608/1‑149 (MQ=255)
                                                                             gttCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGa                                                                 >  1:825629/1‑149 (MQ=255)
                                                                             gttCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGa                                                                 >  1:824362/1‑149 (MQ=255)
                                                                             gttCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGa                                                                 <  2:608885/149‑1 (MQ=255)
                                                                             gttCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGa                                                                 <  1:1500572/149‑1 (MQ=255)
                                                                                ccGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATcc                                                              >  1:1390232/1‑149 (MQ=255)
                                                                                   ccTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTTAAGTAATCCATAACATTCTCAGGGCTGCGATccgcc                                                           >  2:944599/1‑149 (MQ=255)
                                                                                    cTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATccgccg                                                          >  1:1104204/1‑149 (MQ=255)
                                                                                             ggtGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACgg                                                 <  2:757445/147‑1 (MQ=255)
                                                                                             gCTGGCGAATTTTTCATGATGATACTGACGTTGTAGATCCTGATAACGCAGGATAAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGTGTGAAGTAATACATAACATTCTCAGGTCTGCGATCCGCCGCGATAACgg                                                 <  2:1555309/149‑1 (MQ=255)
                                                                                               tGGCGAATTTATCAGTATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTGTCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCg                                               <  2:1012612/149‑1 (MQ=255)
                                                                                                ggCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGa                                              >  2:539822/1‑149 (MQ=255)
                                                                                                   gAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATc                                           >  2:483635/1‑149 (MQ=255)
                                                                                                      tttATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAgg                                        <  2:1383011/149‑1 (MQ=255)
                                                                                                         aTCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCATGTaa                                     >  2:1013321/1‑149 (MQ=255)
                                                                                                         aTCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTaa                                     >  1:296927/1‑149 (MQ=255)
                                                                                                         aTCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTaa                                     <  2:1561588/149‑1 (MQ=255)
                                                                                                          tCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAACCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTaaa                                    <  1:594890/149‑1 (MQ=255)
                                                                                                            aGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATATCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAaaac                                  <  2:1265993/149‑1 (MQ=255)
                                                                                                              gatgatACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAaaacaa                                >  2:100795/1‑149 (MQ=255)
                                                                                                              gatgatACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAaaacaa                                >  1:1307406/1‑149 (MQ=255)
                                                                                                              gatgatACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAaaacaa                                >  2:1594839/1‑149 (MQ=255)
                                                                                                               atgatACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAaaacaaa                               <  1:413873/149‑1 (MQ=255)
                                                                                                                  atACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAAccc                            <  2:825629/149‑1 (MQ=255)
                                                                                                                      tGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTCTCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCAGAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCAcc                        >  1:64995/1‑148 (MQ=255)
                                                                                                                             gTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGc                 >  2:1501426/1‑149 (MQ=255)
                                                                                                                             gTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCACCCGCGCAGGc                 >  2:873062/1‑149 (MQ=255)
                                                                                                                               aGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAACGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGa               <  1:525005/149‑1 (MQ=255)
                                                                                                                               aGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGa               <  2:254860/149‑1 (MQ=255)
                                                                                                                                gATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGTGATCCGCCGCGATACCGGCGAATCAGGTAAAACATACCCCACGCGCCGGCGAc              >  1:123065/1‑149 (MQ=255)
                                                                                                                                 aTCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGAcc             >  1:906140/1‑149 (MQ=255)
                                                                                                                                  tCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGAccc            >  1:328690/1‑148 (MQ=255)
                                                                                                                                  tCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCa            >  2:1543986/1‑149 (MQ=255)
                                                                                                                                    cTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGATGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGGCCACg          >  2:290594/1‑149 (MQ=255)
                                                                                                                                     tGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAATTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACgg         >  2:902623/1‑149 (MQ=255)
                                                                                                                                     tGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATACCGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACgg         >  1:1058618/1‑149 (MQ=255)
                                                                                                                                     tGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACgg         >  1:565375/1‑149 (MQ=255)
                                                                                                                                      gAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCCCggg        >  2:171455/1‑148 (MQ=255)
                                                                                                                                        atacccAGGCGCTGAGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTgg      >  1:261408/6‑149 (MQ=255)
                                                                                                                                        atacccAGGCGCTGAGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTgg      <  1:1560412/144‑1 (MQ=255)
                                                                                                                                        atacccAGGCGCTGAGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTgg      <  1:29902/144‑1 (MQ=255)
                                                                                                                                        atacccAGGCGCTGAGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTgg      >  2:272660/6‑149 (MQ=255)
                                                                                                                                          aCGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTgggg    <  2:511819/149‑1 (MQ=255)
                                                                                                                                           cccAGGCGCTGAGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGt   >  2:687738/3‑149 (MQ=255)
                                                                                                                                           cGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCAGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGt   >  1:241571/1‑149 (MQ=255)
                                                                                                                                            gCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGtt  >  1:1317250/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TTCCGCGCGCATTAACGGATGACGCAGCGTTTGCCAGGCCTGGTTAATGGTTGCAGATTGCTGTACGGCGGCGAGTTGTTCCGCCTGGCTTCCGCTGGCGAATTTATCAGGATGATACTGACGTTGTAGATCCTGAAAACGCAGGCTCAGCGCCTGGGTATCGAGTTGATAGCGGGCAGGCAAGCCAAAGAGGGTGAAGTAATCCATAACATTCTCAGGGCTGCGATCCGCCGCGATAACGGCGAATCAGGTAAAACAAACCCCACGCGCAGGCGACCACGGTGGGGTT  >  NC_000913/2659261‑2659549

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