Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 25,028 G→C 7.3% A880P (GCT→CCT) ‡ ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,0280GC 7.3% 412.7 / 18.0 124A880P (GCT→CCT) ‡ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (62/53);  new base C (6/3);  total (68/56)
Fisher's exact test for biased strand distribution p-value = 5.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.26e-01

Only 100 of 127 total aligned reads displayed.

AGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCC  >  NC_000913/24879‑25177
                                                                                                                                                     |                                                                                                                                                     
aGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                                                                       >  2:767179/1‑150 (MQ=255)
aGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTAGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                                                                       >  2:948019/1‑150 (MQ=255)
   caGTGGGTGGCTCGCTGGAAGCGTCAGTAACCTTGGATGCAGAACCGGAACTGTCGGCGAACCTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTc                                                                                                                                                    <  2:1643069/149‑1 (MQ=255)
      tgggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagc                                                                                                                                                 >  2:729305/1‑150 (MQ=255)
      ggggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagc                                                                                                                                                 <  1:316269/149‑1 (MQ=255)
       gggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagca                                                                                                                                                >  1:1032408/1‑150 (MQ=255)
       gggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagca                                                                                                                                                <  1:1337923/150‑1 (MQ=255)
        ggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagcag                                                                                                                                               >  1:1282490/1‑150 (MQ=255)
        ggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagcag                                                                                                                                               >  2:1307616/1‑150 (MQ=255)
         gtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAACCTGACCGCGCTGGGCGATGAATTACGAGTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCaga                                                                                                                                              >  2:540026/1‑150 (MQ=255)
         gtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagcagc                                                                                                                                              >  2:1517613/1‑149 (MQ=255)
         gtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCaga                                                                                                                                              >  2:14061/1‑150 (MQ=255)
          tggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGCCCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCagag                                                                                                                                             >  1:25990/1‑150 (MQ=255)
          gggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCagag                                                                                                                                             <  2:1287451/149‑1 (MQ=255)
           ggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtg                                                                                                                                            <  1:713426/150‑4 (MQ=255)
            gCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCg                                                                                                                                           >  1:81738/1‑150 (MQ=255)
            gCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCg                                                                                                                                           >  1:68318/1‑150 (MQ=255)
             cgcGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGa                                                                                                                                          <  2:294190/148‑1 (MQ=255)
             cTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGa                                                                                                                                          <  2:975015/150‑1 (MQ=255)
             cTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGa                                                                                                                                          <  1:402435/150‑1 (MQ=255)
              tCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtgcgt                                                                                                                                         <  2:45643/150‑7 (MQ=255)
               cGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAg                                                                                                                                        <  1:1555925/150‑1 (MQ=255)
                 cTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTa                                                                                                                                      >  2:262017/1‑150 (MQ=255)
                 cTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTa                                                                                                                                      >  2:227417/1‑150 (MQ=255)
                  tGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGtgc                                                                                                                                     <  2:1370247/150‑3 (MQ=255)
                  tGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTAc                                                                                                                                     >  1:1551945/1‑150 (MQ=255)
                    gAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTc                                                                                                                                   >  1:1517966/1‑150 (MQ=255)
                      aGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCaa                                                                                                                                 <  2:1104553/150‑1 (MQ=255)
                        cGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAg                                                                                                                               >  1:1718128/1‑150 (MQ=255)
                         ggCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCGTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                              >  2:656357/1‑150 (MQ=255)
                         ggCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                              <  1:556622/150‑1 (MQ=255)
                         ggCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                              <  1:563419/150‑1 (MQ=255)
                          gCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGAGGAATTACGATTTGTCCTGTTGACCTCCGGCTCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAggg                                                                                                                             <  1:1543494/150‑1 (MQ=255)
                               aaCCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                        >  2:320785/1‑150 (MQ=255)
                                 ccTTGTATGCAGAACCGGAATTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAg                                                                                                                      >  2:755603/1‑150 (MQ=255)
                                  cTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGt                                                                                                                     >  2:1392704/1‑150 (MQ=255)
                                   ttGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTc                                                                                                                    <  1:390700/150‑1 (MQ=255)
                                      tATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTcgcg                                                                                                                 <  1:1232705/150‑1 (MQ=255)
                                        tGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCCCCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGtt                                                                                                               >  2:1614704/1‑150 (MQ=255)
                                         gCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTg                                                                                                              >  1:1496755/1‑150 (MQ=255)
                                         gCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTg                                                                                                              >  1:660448/1‑150 (MQ=255)
                                           aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAg                                                                                                            <  2:1250378/150‑1 (MQ=255)
                                            gAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGATGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGt                                                                                                           >  1:1055925/1‑150 (MQ=255)
                                                cGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAg                                                                                                       <  2:480364/150‑1 (MQ=255)
                                                   aaCTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCg                                                                                                    <  1:1517613/150‑1 (MQ=255)
                                                    aCTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGa                                                                                                   >  2:799313/1‑150 (MQ=255)
                                                      tGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAg                                                                                                 <  1:997211/150‑1 (MQ=255)
                                                       gTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAgg                                                                                                >  1:1622375/1‑150 (MQ=255)
                                                         cGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTg                                                                                              <  2:670763/150‑1 (MQ=255)
                                                           gCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTgag                                                                                            <  1:948019/150‑1 (MQ=255)
                                                                 cTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGc                                                                                      >  2:1146604/1‑150 (MQ=255)
                                                                   gACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGccc                                                                                    <  2:660448/150‑1 (MQ=255)
                                                                   gACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGccc                                                                                    >  2:1440439/1‑150 (MQ=255)
                                                                    aCCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCa                                                                                   <  1:262017/150‑1 (MQ=255)
                                                                      cgcgCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACg                                                                                 <  2:1465812/150‑1 (MQ=255)
                                                                       gcgcTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgc                                                                                <  1:755603/150‑1 (MQ=255)
                                                                         gcTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctg                                                                              <  2:719381/150‑1 (MQ=255)
                                                                          cTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgc                                                                             >  2:1492146/1‑150 (MQ=255)
                                                                          cTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgc                                                                             <  1:1705014/150‑1 (MQ=255)
                                                                           tGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgct                                                                            <  1:415700/150‑1 (MQ=255)
                                                                           tGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgct                                                                            <  1:1664762/150‑1 (MQ=255)
                                                                               cGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCa                                                                        <  2:1550213/150‑1 (MQ=255)
                                                                                    aaTTACGATATGACCTGTTGACCTCCGGCGCTACCGTGGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTaca                                                                   <  2:1282490/150‑1 (MQ=255)
                                                                                        acgaggagTCCTGTTCCACTCCGGCGCAACCGTTGCAGAATATACCTACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAACGGGCTCAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGGGCCCACGCGGCTGGCACTACACCCa                                                               <  1:1061515/143‑1 (MQ=255)
                                                                                         cGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGACGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCCCGCTGCTGGCACTACCCCCac                                                              >  2:607328/1‑149 (MQ=255)
                                                                                          gagttGTCCTGTGGACCTCCGGCGCTACCGGTGCAGACTATAACGACGCACCTGCTGATGCACAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGATGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAgg                                                             <  1:110954/147‑1 (MQ=255)
                                                                                               gTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTc                                                        >  1:1589346/1‑150 (MQ=255)
                                                                                                 ccTGTTGACCTCCGGCGCTACCGTTGCATCCTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCgg                                                      <  2:1055925/150‑1 (MQ=255)
                                                                                                 ccTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCgg                                                      >  2:585529/1‑150 (MQ=255)
                                                                                                 ccTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCgg                                                      <  2:81738/150‑1 (MQ=255)
                                                                                                    gTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCaa                                                   <  2:1718128/150‑1 (MQ=255)
                                                                                                     ttGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAg                                                  >  2:42186/1‑150 (MQ=255)
                                                                                                       gACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGt                                                >  1:1091085/1‑150 (MQ=255)
                                                                                                        aCCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTg                                               >  1:1290563/1‑150 (MQ=255)
                                                                                                         ccTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTCAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGGCGGCAAGGTgg                                              >  2:1519587/1‑150 (MQ=255)
                                                                                                          cTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGc                                             >  2:1536779/1‑150 (MQ=255)
                                                                                                            ccGGAGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCgg                                           <  2:1064005/150‑1 (MQ=255)
                                                                                                              ggCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGaa                                         >  2:428933/1‑150 (MQ=255)
                                                                                                                 gcTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAacac                                      >  2:1019339/1‑150 (MQ=255)
                                                                                                                    aCCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCa                                   <  1:1146604/150‑1 (MQ=255)
                                                                                                                    aCCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCa                                   >  1:339283/1‑150 (MQ=255)
                                                                                                                    aCCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCa                                   <  2:575547/150‑1 (MQ=255)
                                                                                                                      cGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGa                                 <  1:1374767/150‑1 (MQ=255)
                                                                                                                         tGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAAt                              >  1:245442/1‑150 (MQ=255)
                                                                                                                           cAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTTCTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCt                            >  1:732979/1‑150 (MQ=255)
                                                                                                                               cTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCgg                        <  1:1022626/150‑1 (MQ=255)
                                                                                                                                tataACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTAGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCAAACGCTCCTGACCCTACACCCAGGATGTCGGCAAGGTGGCGTAACACGCAGAAATCTGCGGc                       >  2:1059742/1‑150 (MQ=255)
                                                                                                                                 ataACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGcc                      >  1:732243/1‑150 (MQ=255)
                                                                                                                                   aaCGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGc                    <  1:773060/150‑1 (MQ=255)
                                                                                                                                   aaCGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGc                    >  2:1123177/1‑150 (MQ=255)
                                                                                                                                         gCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTc              <  1:1307616/150‑1 (MQ=255)
                                                                                                                                          gctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtg                  <  1:751290/142‑1 (MQ=255)
                                                                                                                                           ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGACCCCGCAGAAATCTGCGGCCGCTGTGTCAg            >  2:877241/3‑150 (MQ=255)
                                                                                                                                           ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg            >  2:702338/3‑150 (MQ=255)
                                                                                                                                           ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg            >  1:1201220/3‑150 (MQ=255)
                                                                                                                                           aCCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg            >  2:1522737/1‑150 (MQ=255)
                                                                                                                                            ccTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGc           >  1:984697/1‑150 (MQ=255)
                                                                                                                                             ctgctgATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCa          >  2:1028779/1‑150 (MQ=255)
                                                                                                                                                    tGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGc   >  2:696788/1‑150 (MQ=255)
                                                                                                                                                     gCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGcc  >  1:1214773/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCC  >  NC_000913/24879‑25177

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: