Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 313,726 A→G pseudogene (80/90 nt) ykgP ← pseudogene, oxidoreductase family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913313,7260AG99.4% 611.3 / NA 158pseudogene (80/90 nt)ykgPpseudogene, oxidoreductase family
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (76/81);  minor base T (0/1);  total (76/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

ACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGT  >  NC_000913/313579‑313873
                                                                                                                                                   |                                                                                                                                                   
acatacTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACg                                                                                                                                                  >  1:1227421/1‑151 (MQ=255)
 catacTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACgg                                                                                                                                                 >  1:201451/1‑151 (MQ=255)
       ggATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAAcca                                                                                                                                           <  4:1354539/151‑1 (MQ=255)
       ggATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAAcc                                                                                                                                            <  2:266383/150‑1 (MQ=255)
       ggATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAAcc                                                                                                                                            <  4:1803216/150‑1 (MQ=255)
       ggATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAAcc                                                                                                                                            <  4:343511/150‑1 (MQ=255)
             cacTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGt                                                                                                                                     >  1:187999/1‑151 (MQ=255)
                 gCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTc                                                                                                                                   <  4:233236/149‑1 (MQ=255)
                   aGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTTACGGAAACCAGCCATTTCCttt                                                                                                                               <  1:1337736/151‑1 (MQ=255)
                   aGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCttt                                                                                                                               >  1:315688/1‑151 (MQ=255)
                    ggCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTc                                                                                                                              <  2:282099/151‑1 (MQ=255)
                             tATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCt                                                                                                                       <  4:1505695/149‑1 (MQ=255)
                             tATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTg                                                                                                                      >  4:48224/1‑150 (MQ=255)
                              aTGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGa                                                                                                                     >  4:237187/1‑150 (MQ=255)
                                 agagaACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAAttt                                                                                                                 >  4:423072/1‑151 (MQ=255)
                                   agaACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGa                                                                                                               <  2:166903/151‑1 (MQ=255)
                                         tGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCAt                                                                                                         >  1:363767/1‑151 (MQ=255)
                                          gCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCAt                                                                                                         >  1:1813858/1‑150 (MQ=255)
                                               aGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGttt                                                                                                    >  1:33455/1‑150 (MQ=255)
                                               aGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTa                                                                                                   <  4:1351053/151‑1 (MQ=255)
                                                gTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTAt                                                                                                  <  1:1736760/151‑1 (MQ=255)
                                                 tGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATc                                                                                                 >  1:666917/1‑151 (MQ=255)
                                                 tGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATc                                                                                                 >  1:800705/1‑151 (MQ=255)
                                                   aCGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCg                                                                                               >  1:1495628/1‑151 (MQ=255)
                                                   aCGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCg                                                                                               <  2:57542/151‑1 (MQ=255)
                                                    cGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCGTTTCGATGCCTGAATTTGATCCCATAGTTTATCCg                                                                                               <  2:357868/150‑1 (MQ=255)
                                                    cGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCgg                                                                                              <  4:280406/151‑1 (MQ=255)
                                                     gggTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGa                                                                                             <  4:1834232/151‑1 (MQ=255)
                                                      ggTTACTTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCTTGAATTTGATCCCATAGTTTATCCgg                                                                                              >  1:841650/1‑149 (MQ=255)
                                                       gTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAAt                                                                                           >  1:1109247/1‑151 (MQ=255)
                                                        ttAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATg                                                                                          >  3:126315/1‑151 (MQ=255)
                                                        ttAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATg                                                                                          >  1:1519943/1‑151 (MQ=255)
                                                         tAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATg                                                                                          <  1:337142/150‑1 (MQ=255)
                                                           aTTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGtt                                                                                       >  4:378811/1‑151 (MQ=255)
                                                              aTCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTaa                                                                                     <  2:419964/150‑1 (MQ=255)
                                                                  gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGGATGGTTAATTCAt                                                                                <  1:1364215/151‑1 (MQ=255)
                                                                  gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAAtt                                                                                   <  4:981311/148‑1 (MQ=255)
                                                                  gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCa                                                                                 <  4:86660/150‑1 (MQ=255)
                                                                  gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCAt                                                                                >  3:106622/1‑151 (MQ=255)
                                                                  gAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCAt                                                                                >  1:1603444/1‑151 (MQ=255)
                                                                    aaGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATg                                                                               >  1:198035/1‑150 (MQ=255)
                                                                      gAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAAc                                                                            <  1:371543/151‑1 (MQ=255)
                                                                          ttGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAAGTCATGAACAAg                                                                         <  1:1121970/150‑1 (MQ=255)
                                                                           tgtgAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGt                                                                        <  1:1441770/150‑1 (MQ=255)
                                                                              gAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTtgtg                                                                    <  2:403391/151‑1 (MQ=255)
                                                                              gAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATCCATGAACAAGTtgtg                                                                    >  1:1038685/1‑151 (MQ=255)
                                                                               aGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTgtgt                                                                   >  1:219737/1‑151 (MQ=255)
                                                                                   tGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATc                                                               >  1:1720497/1‑151 (MQ=255)
                                                                                   tGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATc                                                               >  4:1770179/1‑151 (MQ=255)
                                                                                     tATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGt                                                             >  4:61841/1‑151 (MQ=255)
                                                                                     tATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGt                                                             >  2:430581/1‑151 (MQ=255)
                                                                                      aTCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGt                                                             <  4:1437703/150‑1 (MQ=255)
                                                                                       tCGGTTACTAATTGACCTGATAGTTTTACTCTTGGGAGAATAATAGATATTTAACCCTTTGACGGAAACCAGCCAGTTTCTCTAGATGCTTGAATTTGATACCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGtt                                                            >  1:253058/1‑150 (MQ=255)
                                                                                        cGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTc                                                           >  1:582368/1‑150 (MQ=255)
                                                                                        cGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTc                                                           <  1:967895/150‑1 (MQ=255)
                                                                                          gTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATg                                                        >  4:1530388/1‑151 (MQ=255)
                                                                                            tACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATg                                                        >  4:1033288/1‑149 (MQ=255)
                                                                                              cTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATgaga                                                     <  4:1863615/150‑1 (MQ=255)
                                                                                               tCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATgaga                                                     >  1:983366/1‑149 (MQ=255)
                                                                                               tCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATgag                                                      <  4:502130/148‑1 (MQ=255)
                                                                                                    gatttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATaa                                              <  1:1019201/151‑1 (MQ=255)
                                                                                                     atttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAAc                                             <  4:704787/151‑1 (MQ=255)
                                                                                                      tttgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACg                                            >  4:1465739/1‑151 (MQ=255)
                                                                                                        tgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGt                                           >  3:311284/1‑150 (MQ=255)
                                                                                                        tgatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGt                                           <  1:765572/150‑1 (MQ=255)
                                                                                                         gatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAATCCAGACAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTAACCGGAATGGTTAATTCATGAACAAGGTGTTTTATCGTTCATGATAAGCATAACGTaa                                         <  1:743189/151‑1 (MQ=255)
                                                                                                         gatAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTa                                          <  4:1608555/150‑1 (MQ=255)
                                                                                                          atAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTaaa                                        <  1:1620088/151‑1 (MQ=255)
                                                                                                            aGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAgg                                      >  1:989539/1‑151 (MQ=255)
                                                                                                               tttACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGa                                    >  4:331412/1‑150 (MQ=255)
                                                                                                                 tACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGaaa                                  <  2:461006/150‑1 (MQ=255)
                                                                                                                 tACTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGaaa                                  <  3:438522/150‑1 (MQ=255)
                                                                                                                  aCTCTCGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAg                                <  1:1668948/151‑1 (MQ=255)
                                                                                                                      tcGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCg                            >  3:232066/1‑151 (MQ=255)
                                                                                                                       cGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCt                              <  1:92522/148‑1 (MQ=255)
                                                                                                                       cGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTc                             <  1:1476748/149‑1 (MQ=255)
                                                                                                                       cGGGAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTc                             <  4:1678380/149‑1 (MQ=255)
                                                                                                                        gggAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGAt                          <  1:1900906/151‑1 (MQ=255)
                                                                                                                         ggAGAATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGAtt                         <  4:1587268/151‑1 (MQ=255)
                                                                                                                            gaATAATAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGa                      <  1:985373/151‑1 (MQ=255)
                                                                                                                             aataatAGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGAc                     >  2:96565/1‑151 (MQ=255)
                                                                                                                                aataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGc                  >  1:1260727/1‑151 (MQ=255)
                                                                                                                                 ataGATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCa                 <  3:77334/151‑1 (MQ=255)
                                                                                                                                    gATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGaa              <  4:47243/151‑1 (MQ=255)
                                                                                                                                    gATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGa               >  1:565589/1‑150 (MQ=255)
                                                                                                                                    gATATTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGa               >  1:1237554/1‑150 (MQ=255)
                                                                                                                                     atatTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAt             >  4:1636588/1‑151 (MQ=255)
                                                                                                                                     atatTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAt             <  2:303751/151‑1 (MQ=255)
                                                                                                                                       atTTAATCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAAttt           >  3:143724/1‑151 (MQ=255)
                                                                                                                                           aaTCCATTGACGGAAACCCGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTc        >  2:22819/1‑150 (MQ=255)
                                                                                                                                           aaTCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTa                        <  3:259071/134‑1 (MQ=255)
                                                                                                                                           aaTCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTa                        >  2:259079/1‑134 (MQ=255)
                                                                                                                                           aaTCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCa       <  4:169196/151‑1 (MQ=255)
                                                                                                                                           aaTCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCa       <  4:1480077/151‑1 (MQ=255)
                                                                                                                                             tCCATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAgg     >  4:782543/1‑151 (MQ=255)
                                                                                                                                              ccATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAggg    <  4:1122446/151‑1 (MQ=255)
                                                                                                                                              ccATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAgg     <  4:153512/150‑1 (MQ=255)
                                                                                                                                              ccATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAgg     <  1:889316/150‑1 (MQ=255)
                                                                                                                                               cATTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAgggg   <  4:1113779/151‑1 (MQ=255)
                                                                                                                                                aTTGACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGt  >  3:313681/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ACATACTGGATTACACTGCAGGCACGCCTTATGAGAGAACGTGCCGCAGTGACGGGTTAATTATCTGAAAGAATTTGTGAGGCTGTATCGGTTACTCATTGATTTGATAGTTTTACTCTCGGGAGAATAATAGATATTTAATCCATTAACGGAAACCAGCCAGTTCCTTTCGATGCCTGAATTTGATCCCATAGTTTATCCGGAATGGTTAATTCATGAACAAGTTGTGTTATCGTTCATGAGAAGCATAACGTAAAGGGAAAAGCTCGATTAGACGGCAGAATTTGTCAGGGGT  >  NC_000913/313579‑313873

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: