Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,937,909 G→A A206T (GCG→ACG)  rbsK → ribokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,9090GA96.2% 1403.5 / 26.4 431A206T (GCG→ACG) rbsKribokinase
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base A (219/192);  total (231/200)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 431 total aligned reads displayed.

CCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAG  >  NC_000913/3937761‑3938056
                                                                                                                                                    |                                                                                                                                                   
ccATCAAAATAAGACTATTGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACg                                                                                                                                                   <  8:1463028/151‑1 (MQ=255)
ccATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACg                                                                                                                                                   <  6:851633/151‑1 (MQ=255)
ccATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACg                                                                                                                                                   <  2:101773/151‑1 (MQ=255)
 cATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACg                                                                                                                                                   >  5:884148/1‑150 (MQ=255)
 cATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGa                                                                                                                                                  <  3:1613622/151‑1 (MQ=255)
    cAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAgg                                                                                                                                               <  7:1297160/151‑1 (MQ=255)
      aaaTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAggcg                                                                                                                                             <  7:780546/151‑1 (MQ=255)
        aTAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACAGCTGGCGCTGGTGGACATTATTAAGACAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCATCGCAggcggg                                                                                                                                           >  5:1014967/1‑150 (MQ=255)
          aaGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGgcgc                                                                                                                                         <  1:1388310/151‑1 (MQ=255)
          aaGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGgcgc                                                                                                                                         >  5:1639117/1‑151 (MQ=255)
             aCTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAg                                                                                                                                       <  6:1684458/150‑1 (MQ=255)
              cTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGt                                                                                                                                     >  7:903947/1‑151 (MQ=255)
              cTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGt                                                                                                                                     >  3:1190955/1‑151 (MQ=255)
               tATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGt                                                                                                                                     >  6:713703/1‑150 (MQ=255)
                aTCGTTGCGGTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCACACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGTCCGCGCAGGTAc                                                                                                                                   >  4:185299/1‑151 (MQ=255)
                aTCGTTGCGCTTAACGCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATTAAGATGCAACGAAGGCGGCGCAGGt                                                                                                                                     <  3:1302338/149‑1 (MQ=255)
                   gTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGc                                                                                                                                <  5:988165/151‑1 (MQ=255)
                     tGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCAt                                                                                                                              >  3:1457473/1‑151 (MQ=255)
                      gcgcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCAt                                                                                                                              <  5:231085/150‑1 (MQ=255)
                      gcgcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATg                                                                                                                             >  6:1090923/1‑151 (MQ=255)
                        gcTTAACCGGGCTCGGGCTCGCGAAGTTCCTGCCGACCTGCTGGGGCTGGTGGTCATTATTTCGCCTATGGAAAGGGAAGCATAAAAGCTCACGGGTATTCGTGTTGAAATTGATGACGATGCATCGAAGGCGGCGCAGGTACTGCATGaa                                                                                                                           <  2:551217/151‑1 (MQ=255)
                        gcTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGa                                                                                                                            >  8:1043788/1‑150 (MQ=255)
                          ttAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGaaa                                                                                                                          >  3:681275/1‑150 (MQ=255)
                           tAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGaaaa                                                                                                                         <  1:1078080/150‑1 (MQ=255)
                            aaCCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGaaaa                                                                                                                         >  8:153038/1‑149 (MQ=255)
                             aCCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAgg                                                                                                                      >  6:495562/1‑151 (MQ=255)
                              cccGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGt                                                                                                                     <  5:1652826/151‑1 (MQ=255)
                               ccggctccggctCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTa                                                                                                                    <  2:1232782/151‑1 (MQ=255)
                                cggctccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTAt                                                                                                                   >  6:1294232/1‑151 (MQ=255)
                                 ggctccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGCTGCACCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATc                                                                                                                  >  5:503990/1‑151 (MQ=255)
                                  gctccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATcc                                                                                                                 <  4:645423/151‑1 (MQ=255)
                                    gccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGACACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCgt                                                                                                               <  4:1017247/150‑1 (MQ=255)
                                     ccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCgta                                                                                                              <  3:448250/151‑1 (MQ=255)
                                     ccggctcGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCg                                                                                                                <  3:178877/149‑1 (MQ=255)
                                           cgcgAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCgtactgt                                                                                                          <  1:1242285/149‑1 (MQ=255)
                                              gAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGtactctactg                                                                                                      >  5:579061/1‑150 (MQ=255)
                                                 attCCTGACGAACTGCAGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTa                                                                                                  <  8:1084329/150‑1 (MQ=255)
                                                  ttCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTa                                                                                                  >  4:367885/1‑150 (MQ=255)
                                                    ccTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACttt                                                                                              <  3:1140394/151‑1 (MQ=255)
                                                       gACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACttt                                                                                              >  5:1379130/1‑149 (MQ=255)
                                                       gACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTa                                                                                             <  4:1585556/150‑1 (MQ=255)
                                                       gACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAg                                                                                            <  2:399758/151‑1 (MQ=255)
                                                       gACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATCCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTa                                                                                             <  6:872174/150‑1 (MQ=255)
                                                        aCGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAgg                                                                                           >  5:69693/1‑151 (MQ=255)
                                                               gctgGCTCTGGTTGACATTATTACTCCAAACGAAACGTAAGCAGAAAAGCTCTCCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCg                                                                                     <  4:83479/150‑1 (MQ=255)
                                                               gctgGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTc                                                                                      <  2:340571/149‑1 (MQ=255)
                                                               gctctgggtggtggTCATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCTGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTc                                                                                      <  8:1166554/141‑1 (MQ=255)
                                                                   gcgcTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCgtggtg                                                                                <  7:1557720/151‑1 (MQ=255)
                                                                       tggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGTTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATg                                                                            >  2:1512279/1‑151 (MQ=255)
                                                                       tggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATg                                                                            >  6:967135/1‑151 (MQ=255)
                                                                       tggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATg                                                                            >  1:157910/1‑151 (MQ=255)
                                                                        ggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATgg                                                                           <  5:1130758/151‑1 (MQ=255)
                                                                        ggtggACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATgg                                                                           >  1:266450/1‑151 (MQ=255)
                                                                             aCATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCt                                                                        >  3:29773/1‑149 (MQ=255)
                                                                             aCATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCt                                                                        >  2:315499/1‑149 (MQ=255)
                                                                               attattACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCg                                                                    >  1:866581/1‑151 (MQ=255)
                                                                                 tattaCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGt                                                                   <  7:847256/150‑1 (MQ=255)
                                                                                     aCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAAt                                                               <  5:1469745/150‑1 (MQ=255)
                                                                                     aCGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATg                                                              <  4:27863/151‑1 (MQ=255)
                                                                                      cGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTAATCGTGTTCAAAATGATGAAGATGCAACGAAGGCGGCGCAGCTACTGCAGGAAAAAGCTCTCCGTACTGTACTGATTACTTTAGGAAGTCCAGGGTTAAGGGCTAGCGTGAATg                                                              >  7:627684/1‑150 (MQ=255)
                                                                                         cAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTTCTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAAtgatga                                                          <  2:941762/151‑1 (MQ=255)
                                                                                          aaacgaaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGaa                                                         <  2:1428979/151‑1 (MQ=255)
                                                                                          aaacgaaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGAGTGAAtgatga                                                          >  4:1332217/1‑150 (MQ=255)
                                                                                          aaacgaaacgCAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGaa                                                         <  2:594273/151‑1 (MQ=255)
                                                                                           aacgaaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAAtgatga                                                          >  7:1060270/1‑149 (MQ=255)
                                                                                           aacgaaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAg                                                        <  5:1342186/151‑1 (MQ=255)
                                                                                             cgaaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGt                                                      <  8:464029/151‑1 (MQ=255)
                                                                                              gaaacgGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCTTGGTGTATGGGCTAGCGTGAATGATGAAGGTc                                                     >  7:1017040/1‑151 (MQ=255)
                                                                                                     aaGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGt                                              >  6:616936/1‑151 (MQ=255)
                                                                                                      aGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGtt                                             >  8:520107/1‑151 (MQ=255)
                                                                                                         aGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCt                                          <  8:626935/151‑1 (MQ=255)
                                                                                                          gAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTg                                         >  6:343063/1‑151 (MQ=255)
                                                                                                          gAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTg                                         >  4:736825/1‑151 (MQ=255)
                                                                                                              aGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGAtt                                     <  4:1209671/151‑1 (MQ=255)
                                                                                                                gtCACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTcc                                   <  4:906524/150‑1 (MQ=255)
                                                                                                                  cACCGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTcc                                   <  3:355065/149‑1 (MQ=255)
                                                                                                                     cGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTg                              >  3:1422901/1‑151 (MQ=255)
                                                                                                                     cGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTg                              <  3:839300/151‑1 (MQ=255)
                                                                                                                     cGGTATTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTg                              >  1:444351/1‑151 (MQ=255)
                                                                                                                         aTTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAgg                          <  8:788097/151‑1 (MQ=255)
                                                                                                                         aTTCGTGTTGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTCGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAg                           <  8:1403026/150‑1 (MQ=255)
                                                                                                                                 tGAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGAt                  >  7:957140/1‑151 (MQ=255)
                                                                                                                                  gAAAATGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGAGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGAt                  >  8:317515/1‑150 (MQ=255)
                                                                                                                                   aaaaTGATGAGGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATAc                >  7:1452241/1‑151 (MQ=255)
                                                                                                                                   aaaaTGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATAc                <  7:444591/151‑1 (MQ=255)
                                                                                                                                   aaaaTGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATAc                <  2:183950/151‑1 (MQ=255)
                                                                                                                                     aaTGATGAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGATCCTGGATGCCGGGTGCGGGCTGTCGATccc               >  1:1134104/1‑150 (MQ=255)
                                                                                                                                      atgatgAAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCAt             <  5:313674/151‑1 (MQ=255)
                                                                                                                                       tgatgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCAtt            >  1:603011/1‑151 (MQ=255)
                                                                                                                                        gatgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGCTACCATtg           >  6:1460632/1‑151 (MQ=255)
                                                                                                                                         atgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATtgg          >  2:1110057/1‑150 (MQ=255)
                                                                                                                                         atgaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATtgc          <  2:442782/151‑1 (MQ=255)
                                                                                                                                           gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATtgctg        <  3:563722/151‑1 (MQ=255)
                                                                                                                                           gaAGATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATtgctg        <  1:1013665/151‑1 (MQ=255)
                                                                                                                                              gATGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGGGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCg     <  8:1689049/151‑1 (MQ=255)
                                                                                                                                               aTGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCgg    >  7:174664/1‑151 (MQ=255)
                                                                                                                                               aTGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCg     <  4:487401/150‑1 (MQ=255)
                                                                                                                                                tGCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCgg    >  2:1321823/1‑150 (MQ=255)
                                                                                                                                                 gCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGgag  <  2:85181/151‑1 (MQ=255)
                                                                                                                                                 gCAACGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGATGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGga   <  6:383288/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCATCAAAATAAGACTATCGTTGCGCTTAACCCGGCTCCGGCTCGCGAACTTCCTGACGAACTGCTGGCGCTGGTGGACATTATTACGCCAAACGAAACGGAAGCAGAAAAGCTCACCGGTATTCGTGTTGAAAATGATGAAGATGCAGCGAAGGCGGCGCAGGTACTGCATGAAAAAGGTATCCGTACTGTACTGATTACTTTAGGAAGTCGTGGTGTATGGGCTAGCGTGAATGGTGAAGGTCAGCGCGTTCCTGGATTCCGGGTGCAGGCTGTCGATACCATTGCTGCCGGAG  >  NC_000913/3937761‑3938056

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