Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 25,027 T→G 9.0% D879E (GAT→GAG)  ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091325,0270TG 9.0% 260.5 / 17.0 89D879E (GAT→GAG) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base T (32/48);  new base G (5/3);  total (38/51)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

AAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAA  >  NC_000913/24878‑25170
                                                                                                                                                     |                                                                                                                                               
aaGAATGTGGGTGGCTCGCTGGAATCGGCTGAAACCTTGTATGCAGAACCGGATCTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAt                                                                                                                                                 <  2:1253707/150‑1 (MQ=255)
aaGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAt                                                                                                                                                 <  1:1950937/150‑1 (MQ=255)
 aGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                                                                <  1:1668508/150‑1 (MQ=255)
 aGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                                                                >  2:138761/1‑150 (MQ=255)
 aGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATg                                                                                                                                                >  1:1581575/1‑150 (MQ=255)
  gAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGc                                                                                                                                               >  1:46362/1‑150 (MQ=255)
     aGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGATCTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGCCTATAACGACGCACCTGCTGATGCTca                                                                                                                                            >  1:623424/1‑150 (MQ=255)
       tgggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagc                                                                                                                                          <  2:1825960/150‑1 (MQ=255)
        gggtggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTcagca                                                                                                                                         >  2:1879683/1‑150 (MQ=255)
            ggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtg                                                                                                                                     <  2:1049310/150‑4 (MQ=255)
            ggCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtg                                                                                                                                     <  1:1901447/150‑4 (MQ=255)
              cTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGa                                                                                                                                   <  2:664631/150‑1 (MQ=255)
               tCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtg                                                                                                                                     >  2:1568762/1‑144 (MQ=255)
               tCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGAGCATcagcaggtg                                                                                                                                     <  1:1568762/147‑4 (MQ=255)
                 gCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGt                                                                                                                                >  1:1976775/1‑150 (MQ=255)
                   gggAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTAc                                                                                                                              <  2:520210/149‑1 (MQ=255)
                    ggAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACt                                                                                                                             <  1:559878/150‑1 (MQ=255)
                      aaGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCa                                                                                                                           <  2:200945/150‑1 (MQ=255)
                      aaGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCa                                                                                                                           <  2:1193678/150‑1 (MQ=255)
                       aGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCaa                                                                                                                          >  1:1910586/1‑150 (MQ=255)
                         cGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACAACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAg                                                                                                                        <  1:966116/150‑1 (MQ=255)
                          ggCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAgg                                                                                                                       <  1:1061526/150‑1 (MQ=255)
                                aaCCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                 >  2:1683706/1‑150 (MQ=255)
                                aaCCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                 <  1:182334/150‑1 (MQ=255)
                                aaCCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGaa                                                                                                                 <  2:561834/150‑1 (MQ=255)
                                   cTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGt                                                                                                              <  2:1581575/150‑1 (MQ=255)
                                      gTATGCAGAACAGGAACTGTCGGCGAAACTGAACGCGCTGGCCGATGAATTACGATTTGTACTGTTCACCTCCGGCGCTACCGTTGAAGACTAAAACGACGCACCTGCTGAAGCTCACCAGAGCGAAGTACTCAAAGGCCTGAAAGTcgc                                                                                                           >  2:345222/1‑150 (MQ=255)
                                        aTGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGt                                                                                                         >  1:1592464/1‑150 (MQ=255)
                                        aTGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGt                                                                                                         >  1:763549/1‑150 (MQ=255)
                                          gCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTg                                                                                                       <  1:204101/150‑1 (MQ=255)
                                          gCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTg                                                                                                       <  1:945197/150‑1 (MQ=255)
                                           cAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                                                      >  2:1148255/1‑150 (MQ=255)
                                           cAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                                                      >  2:244725/1‑150 (MQ=255)
                                           cAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGa                                                                                                      >  2:80946/1‑150 (MQ=255)
                                            aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAg                                                                                                     <  2:1965437/150‑1 (MQ=255)
                                            aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAg                                                                                                     <  1:1439130/150‑1 (MQ=255)
                                            aGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCGGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAg                                                                                                     <  2:1664635/150‑1 (MQ=255)
                                             gAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGt                                                                                                    <  1:80946/150‑1 (MQ=255)
                                               aCCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTaa                                                                                                  >  2:1633458/1‑150 (MQ=255)
                                                 cGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAg                                                                                                <  1:649753/150‑1 (MQ=255)
                                                    aaCTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCg                                                                                             <  1:1879683/150‑1 (MQ=255)
                                                     aCTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGa                                                                                            <  2:1910586/150‑1 (MQ=255)
                                                       tGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAg                                                                                          <  2:736820/150‑1 (MQ=255)
                                                        gTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAgg                                                                                         <  1:1345718/150‑1 (MQ=255)
                                                          cGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTg                                                                                       <  2:1832781/150‑1 (MQ=255)
                                                               aaaCTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAg                                                                                  <  1:160255/150‑1 (MQ=255)
                                                                 aCTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTg                                                                                <  1:1427366/150‑1 (MQ=255)
                                                                   tGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTAGCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAATGTGAGAAGTGcc                                                                              <  2:1516860/150‑1 (MQ=255)
                                                                    gACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGccc                                                                             <  2:627899/150‑1 (MQ=255)
                                                                     aCCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCa                                                                            >  2:1643767/1‑150 (MQ=255)
                                                                        gcgcTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACgc                                                                         <  2:1592464/150‑1 (MQ=255)
                                                                              ggCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCTTTTCAGACTATAACGACGCACCTGCGGATGCTAAGCAGAGCGATGTACTCAAAGGGATGATAGTCTAGTTGATTAAAGCCTATGGTGATAAGTGCCCAAgctgcttt                                                                   <  2:269454/150‑3 (MQ=255)
                                                                               gCGATGAATTACGATTTGTCCTGTTGCCCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGc                                                                  >  2:115173/1‑150 (MQ=255)
                                                                               gCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGc                                                                  <  1:812380/150‑1 (MQ=255)
                                                                                     aaTTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGtt                                                                                                        <  2:875108/106‑1 (MQ=255)
                                                                                     aaTTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGtt                                                                                                        >  1:875108/1‑106 (MQ=255)
                                                                                        tACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACAccc                                                         <  2:1597227/150‑1 (MQ=255)
                                                                                          ggATTTGTGCTGTTGACCTCCGGCGCTACCGTTGGAGTCTATATCGACGCACGTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAg                                                       <  1:345222/149‑1 (MQ=255)
                                                                                           gATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCATCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAgg                                                      <  2:46362/150‑1 (MQ=255)
                                                                                           gATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACt                                                               <  1:271960/141‑1 (MQ=255)
                                                                                           gATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACt                                                               >  2:271960/1‑141 (MQ=255)
                                                                                               tGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGt                                                  >  2:1300139/1‑150 (MQ=255)
                                                                                                gTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTc                                                 <  1:1633458/150‑1 (MQ=255)
                                                                                                gTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTc                                                 >  1:807252/1‑150 (MQ=255)
                                                                                                  ccTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCCCTACACCCAGGATGTCgg                                               >  1:1917753/1‑150 (MQ=255)
                                                                                                    tGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCa                                             >  2:1546515/1‑150 (MQ=255)
                                                                                                            tCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCg                                     >  2:686894/1‑150 (MQ=255)
                                                                                                             ccGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCgg                                    <  2:744826/150‑1 (MQ=255)
                                                                                                             ccGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCgg                                    <  1:1507912/150‑1 (MQ=255)
                                                                                                                gcgcTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAac                                 <  2:1976775/150‑1 (MQ=255)
                                                                                                                     aCCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCa                            >  1:502145/1‑150 (MQ=255)
                                                                                                                      ccGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAg                           <  1:1218502/150‑1 (MQ=255)
                                                                                                                             aGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTg                    >  2:1571938/1‑150 (MQ=255)
                                                                                                                               aCTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCg                  >  1:156272/1‑150 (MQ=255)
                                                                                                                                   taACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGc                >  2:1127164/1‑148 (MQ=255)
                                                                                                                                   taACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGc                <  1:1127164/148‑1 (MQ=255)
                                                                                                                                   taACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCg              >  2:353694/1‑150 (MQ=255)
                                                                                                                                      cgacgCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtg           <  1:1457073/150‑1 (MQ=255)
                                                                                                                                        acgCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtgtg         <  1:1300139/150‑1 (MQ=255)
                                                                                                                                        acgCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCtgtg         <  1:686894/150‑1 (MQ=255)
                                                                                                                                         cgCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTgtgt        >  2:177094/1‑150 (MQ=255)
                                                                                                                                         cgCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTgtgt        >  2:1349706/1‑150 (MQ=255)
                                                                                                                                          gCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTc       <  2:763549/150‑1 (MQ=255)
                                                                                                                                            ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGACCACGCAGAAATCTGCGGCCGCTGTGTCAg     >  1:778942/3‑150 (MQ=255)
                                                                                                                                            ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg     >  1:665399/3‑150 (MQ=255)
                                                                                                                                            ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg     >  1:1138271/3‑150 (MQ=255)
                                                                                                                                            ctctgctgAGCATCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAg     >  2:1727726/3‑150 (MQ=255)
                                                                                                                                              ctgctgATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCa   <  1:981275/150‑1 (MQ=255)
                                                                                                                                               tgctgATGCTCAGAACAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAACCTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGACAGCaa  >  2:427255/1‑150 (MQ=255)
                                                                                                                                                gctgATGCTCAGCAGAGCGAAGTACTCATAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGa                                   <  2:108629/116‑1 (MQ=255)
                                                                                                                                                gctgATGCTCAGCAGAGCGAAGTACTCAGAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGa                                   >  1:108629/1‑116 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AAGAAAGTGGGTGGCTCGCTGGAAGCGGCAGTAACCTTGTATGCAGAACCGGAACTGTCGGCGAAACTGACCGCGCTGGGCGATGAATTACGATTTGTCCTGTTGACCTCCGGCGCTACCGTTGCAGACTATAACGACGCACCTGCTGATGCTCAGCAGAGCGAAGTACTCAAAGGGCTGAAAGTCGCGTTGAGTAAAGCCGAAGGTGAGAAGTGCCCACGCTGCTGGCACTACACCCAGGATGTCGGCAAGGTGGCGGAACACGCAGAAATCTGCGGCCGCTGTGTCAGCAA  >  NC_000913/24878‑25170

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