Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,013,189 C→G C170S (TGC→TCC)  SBRG_Ecoli_galK ← SBRG_Ecoli_galK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,013,1890CG99.1% 1215.5 / ‑1.0 345C170S (TGC→TCC) SBRG_Ecoli_galKSBRG_Ecoli_galK
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base G (165/177);  total (166/179)
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.96e-01

Only 100 of 345 total aligned reads displayed.

ACCAGGGTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGC  >  NC_002947/1013042‑1013335
                                                                                                                                                   |                                                                                                                                                  
     ggTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACa                                                                                                                                             <  1:4923239/150‑1 (MQ=255)
        aCGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGcc                                                                                                                                          <  1:6838524/150‑1 (MQ=255)
           tttGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTAc                                                                                                                                       >  1:8299900/1‑150 (MQ=255)
              gAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACaaa                                                                                                                                    >  2:7402988/1‑150 (MQ=255)
               aaGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAAc                                                                                                                                   <  2:1356403/150‑1 (MQ=255)
                    aCTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGtt                                                                                                                              >  1:1623843/1‑150 (MQ=255)
                      tGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGtttt                                                                                                                            >  1:3636475/1‑150 (MQ=255)
                      tGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGtttt                                                                                                                            >  1:8675339/1‑150 (MQ=255)
                       gTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTc                                                                                                                           <  1:2452982/150‑1 (MQ=255)
                        ttGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCt                                                                                                                          >  2:6785837/1‑150 (MQ=255)
                         tgatgatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTg                                                                                                                         >  2:3138201/1‑150 (MQ=255)
                             gatgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTc                                                                                                                     <  1:6945404/150‑1 (MQ=255)
                               tgaCGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCt                                                                                                                   <  1:6476853/150‑1 (MQ=255)
                                   gacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGAcc                                                                                                               <  2:5420153/150‑1 (MQ=255)
                                   gacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGAcc                                                                                                               >  1:7351811/1‑150 (MQ=255)
                                   gacAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGAcc                                                                                                               <  2:2265891/150‑1 (MQ=255)
                                      aGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGGTCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGtt                                                                                                            <  1:8537301/150‑1 (MQ=255)
                                      aGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGtt                                                                                                            >  1:5461437/1‑150 (MQ=255)
                                        ccACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                                          >  2:398679/1‑150 (MQ=255)
                                        ccACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                                          <  2:5883783/150‑1 (MQ=255)
                                        ccACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                                          >  1:7641557/1‑150 (MQ=255)
                                        ccACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                                          <  1:7497048/150‑1 (MQ=255)
                                           cacCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAgcg                                                                                                       >  2:1545080/1‑150 (MQ=255)
                                               tttGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGAt                                                                                                   <  2:8516696/150‑1 (MQ=255)
                                                   ggCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTtgt                                                                                               <  2:5131055/150‑1 (MQ=255)
                                                    gCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTtgtg                                                                                              <  1:35563/150‑1 (MQ=255)
                                                    gCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTtgtg                                                                                              >  2:2549912/1‑150 (MQ=255)
                                                     cATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgc                                                                                             <  2:2732024/150‑1 (MQ=255)
                                                     cATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTgc                                                                                             <  1:1991342/150‑1 (MQ=255)
                                                        ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGCTTTGTGCGcc                                                                                          >  2:7831129/1‑150 (MQ=255)
                                                        ggAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGcc                                                                                          <  1:7298601/150‑1 (MQ=255)
                                                         gAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCg                                                                                         >  1:2585620/1‑150 (MQ=255)
                                                              tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                                          >  2:3200263/1‑128 (MQ=255)
                                                              tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTaa                                                                                                          <  1:3200263/128‑1 (MQ=255)
                                                              tGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAg                                                                                    >  1:2244885/1‑150 (MQ=255)
                                                               gCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGc                                                                                   >  2:2118369/1‑150 (MQ=255)
                                                               gCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGc                                                                                   <  1:2962431/150‑1 (MQ=255)
                                                                 tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAgc                                                                                                        >  1:5150431/1‑127 (MQ=255)
                                                                 tttGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCgg                                                                                 >  2:8093310/1‑150 (MQ=255)
                                                                     gTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAga                                                                             >  1:1813377/1‑150 (MQ=255)
                                                                         ccAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAgatgat                                                                         <  2:7361309/150‑1 (MQ=255)
                                                                               aGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgct                                                                   <  1:7849147/150‑1 (MQ=255)
                                                                                 cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctgc                                                                 <  1:3340764/150‑1 (MQ=255)
                                                                                 cGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctgc                                                                 <  2:9448131/150‑1 (MQ=255)
                                                                                  ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctgct                                                                <  2:8831998/150‑1 (MQ=255)
                                                                                  ggCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctgct                                                                <  2:69152/150‑1 (MQ=255)
                                                                                   gCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGctgctg                                                               >  2:9098068/1‑150 (MQ=255)
                                                                                      aTCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCaa                                                            >  1:9550276/1‑150 (MQ=255)
                                                                                      aTCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCaa                                                            <  2:1813377/150‑1 (MQ=255)
                                                                                        cGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATa                                                          >  1:9607130/1‑150 (MQ=255)
                                                                                         gATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCt                                                                                                                          >  1:4226923/1‑85 (MQ=255)
                                                                                         gATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGTGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATAc                                                         >  1:715931/1‑150 (MQ=255)
                                                                                               cAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTcc                                                   <  1:2665060/150‑1 (MQ=255)
                                                                                               cAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTcc                                                   >  2:7226337/1‑150 (MQ=255)
                                                                                                 aGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGa                                                 <  2:8382097/150‑1 (MQ=255)
                                                                                                  ggCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAgctgc                                                                 >  1:1917287/1‑133 (MQ=255)
                                                                                                   gCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGAcc                                               >  2:732539/1‑150 (MQ=255)
                                                                                                     aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgc                                             >  2:6054333/1‑150 (MQ=255)
                                                                                                     aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgc                                             <  2:566971/150‑1 (MQ=255)
                                                                                                     aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgc                                             >  2:187448/1‑150 (MQ=255)
                                                                                                     aTGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgc                                             <  2:7456840/150‑1 (MQ=255)
                                                                                                        aTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGAc                                          <  2:200223/150‑1 (MQ=255)
                                                                                                        aTCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGAc                                          >  1:5716104/1‑150 (MQ=255)
                                                                                                            ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTcc                                      <  2:3610284/150‑1 (MQ=255)
                                                                                                            ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTcc                                      <  2:3920016/150‑1 (MQ=255)
                                                                                                            ttcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTcc                                      >  1:6291684/1‑150 (MQ=255)
                                                                                                             tcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACcgc                                             <  2:813093/142‑1 (MQ=255)
                                                                                                             tcttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCa                                     <  1:8171186/150‑1 (MQ=255)
                                                                                                               ttGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGt                                   <  1:810987/150‑1 (MQ=255)
                                                                                                                  ccGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGaa                                <  2:3540138/150‑1 (MQ=255)
                                                                                                                  ccGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGaa                                <  2:7701104/150‑1 (MQ=255)
                                                                                                                  ccGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGaa                                <  2:904809/150‑1 (MQ=255)
                                                                                                                   cGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAg                               <  2:7008077/150‑1 (MQ=255)
                                                                                                                   cGAGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAg                               <  2:6191667/150‑1 (MQ=255)
                                                                                                                    gAGCACGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGc                              <  2:9571329/150‑1 (MQ=255)
                                                                                                                     aGCGCGGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCg                             >  1:8906843/1‑150 (MQ=255)
                                                                                                                       cgcgGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGga                           <  2:4344049/150‑1 (MQ=255)
                                                                                                                        gcgGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaa                          <  2:706139/150‑1 (MQ=255)
                                                                                                                         cgGAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaag                         >  1:5794344/1‑150 (MQ=255)
                                                                                                                          ggAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaaga                        >  1:3974225/1‑150 (MQ=255)
                                                                                                                          ggAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaaga                        >  2:8705087/1‑150 (MQ=255)
                                                                                                                          ggAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaaga                        >  2:8224812/1‑150 (MQ=255)
                                                                                                                          ggAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaaga                        >  2:7752996/1‑150 (MQ=255)
                                                                                                                           gAAATTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGTGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGgaagaa                       <  2:715931/150‑1 (MQ=255)
                                                                                                                            aaaTTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAAc                      <  1:4126860/150‑1 (MQ=255)
                                                                                                                             aaTTAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACt                     <  1:1036486/150‑1 (MQ=255)
                                                                                                                               ttAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTa                   <  2:5261930/150‑1 (MQ=255)
                                                                                                                                tAGCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTaa                  <  1:1184994/150‑1 (MQ=255)
                                                                                                                                tAGCTGACCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTaa                  <  1:9094564/150‑1 (MQ=255)
                                                                                                                                  gCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAAc                      >  1:4712033/1‑144 (MQ=255)
                                                                                                                                  gCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAAc                      <  2:4712033/144‑1 (MQ=255)
                                                                                                                                  gCTGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGGTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATACAGCGGCTGCACTACGGTTCCGACCGCGACTTCCAGTGAAGCGGACGAACTTAAcc                >  2:2960465/1‑150 (MQ=255)
                                                                                                                                  gCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAAcc                <  1:8848420/150‑1 (MQ=255)
                                                                                                                                  gCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAAcc                <  1:1421618/150‑1 (MQ=255)
                                                                                                                                    tGATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCg              <  2:3654335/150‑1 (MQ=255)
                                                                                                                                     gATCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCgg             >  2:1204303/1‑150 (MQ=255)
                                                                                                                                       tCCATGATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAg                                                                                                         <  1:226336/56‑1 (MQ=255)
                                                                                                                                            gATCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTg      >  1:9010094/1‑150 (MQ=255)
                                                                                                                                             aTCCCGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGc     <  1:6123721/150‑1 (MQ=255)
                                                                                                                                                ccGGAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCGATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGc  >  1:9212471/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ACCAGGGTACGTTTGAAGTTACTGTTGATGATGACGACAGCCACACCTTTGGGCATGGAAACTGCTTTGGTCCCCAGTGAGCGGCAATCGATCAGCAAGGCATGATCTTTCTTGCCGAGCGCGGAAATTAGCTGATCCATGATCCCGCAGTTACAGCCTACAAACTGGTTTTCTGCTTCCTGACCGTTAAGCGCGATTTGTGCGCCGTCCAGCGGCAGATGATAAAGCTGCTGCAATACGGTTCCGACCGCGACTTCCAGTGAAGCGGAAGAACTTAACCCGGCACCCTGCGGC  >  NC_002947/1013042‑1013335

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