Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 25,029 C→A 7.3% A880D (GCT→GAT) ‡ ileS → isoleucyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,0290CA 7.3% 383.5 / 15.0 124A880D (GCT→GAT) ‡ileSisoleucyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (62/53);  new base A (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 = 2.75e-01

Only 100 of 127 total aligned reads displayed.

AGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCG  >  NC_000913/24879‑25178
                                                                                                                                                      |                                                                                                                                                     
   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                                                                                                                                                >  2:1307616/1‑150 (MQ=255)
        ggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagcag                                                                                                                                                >  1:1282490/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                                                                                                                                           <  1:402435/150‑1 (MQ=255)
             cTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGa                                                                                                                                           <  2:975015/150‑1 (MQ=255)
              tCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtgcgt                                                                                                                                          <  2:45643/150‑7 (MQ=255)
               cGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAg                                                                                                                                         <  1:1555925/150‑1 (MQ=255)
                 cTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTa                                                                                                                                       >  2:227417/1‑150 (MQ=255)
                 cTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTa                                                                                                                                       >  2:262017/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:563419/150‑1 (MQ=255)
                         ggCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                               <  1:556622/150‑1 (MQ=255)
                          gCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGAGGAATTACGATTTGTCCTGTTGACCTCCGGCTCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAggg                                                                                                                              <  1:1543494/150‑1 (MQ=255)
                               aaCCTTGTATGCAGAACCGGAACTGTCGGGGAAACTGCCCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                         >  2:832347/1‑150 (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:660448/1‑150 (MQ=255)
                                         gCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTg                                                                                                               >  1:1496755/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:1440439/1‑150 (MQ=255)
                                                                   gACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGccc                                                                                     <  2:660448/150‑1 (MQ=255)
                                                                    aCCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCa                                                                                    <  1:262017/150‑1 (MQ=255)
                                                                      cgcgCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACg                                                                                  <  2:1465812/150‑1 (MQ=255)
                                                                       gcgcTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgc                                                                                 <  1:755603/150‑1 (MQ=255)
                                                                          cTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgc                                                                              <  1:1705014/150‑1 (MQ=255)
                                                                          cTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgctgc                                                                              >  2:1492146/1‑150 (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:68318/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                                    <  2:575547/150‑1 (MQ=255)
                                                                                                                    aCCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCa                                    >  1:339283/1‑150 (MQ=255)
                                                                                                                    aCCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCa                                    <  1:1146604/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                     >  2:1123177/1‑150 (MQ=255)
                                                                                                                                   aaCGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGc                     <  1:773060/150‑1 (MQ=255)
                                                                                                                                         gCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTc               <  1:1307616/150‑1 (MQ=255)
                                                                                                                                          gctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtg                   >  2:751290/4‑145 (MQ=255)
                                                                                                                                          gctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtg                   <  1:751290/142‑1 (MQ=255)
                                                                                                                                           ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGACCCCGCAGAAATCTGCGGCCGCTGTGTCAg             >  2:877241/3‑150 (MQ=255)
                                                                                                                                           ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg             >  1:1201220/3‑150 (MQ=255)
                                                                                                                                           ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg             >  2:702338/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)
                                                                                                                                                      cTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCcg  <  2:708084/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
AGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAACGTCGCCG  >  NC_000913/24879‑25178

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: