Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,343,668 T→G 5.6% L287R (CTG→CGG) ‡ pgi‑II → glucose 6‑phosphate isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,343,6680TG 5.6% 351.9 / 11.2 124L287R (CTG→CGG) ‡pgi‑IIglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base T (53/64);  new base G (2/5);  total (55/69)
Fisher's exact test for biased strand distribution p-value = 4.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 125 total aligned reads displayed.

AGCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGC  >  NC_002947/5343519‑5343816
                                                                                                                                                     |                                                                                                                                                    
   cGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGc                                                                                                                                                   >  1:2565287/1‑150 (MQ=35)
        tGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCg                                                                                                                                              >  2:911560/1‑150 (MQ=35)
        ggTACCGCCACTTCATTGCGTTATCCAGCAACAAGGCCGTGGCTGTAGGGTGGGGTATCCGGGAAGAAAACATTTCCCCAATTGGGGACTGGGTGGGGGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCg                                                                                                                                              <  2:599706/149‑1 (MQ=11)
           cccGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGCa                                                                                                                                           <  2:1833125/149‑1 (MQ=11)
           aCCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGCa                                                                                                                                           <  1:1093762/150‑1 (MQ=11)
           aCCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGCa                                                                                                                                           <  1:308348/150‑1 (MQ=11)
           aCCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCa                                                                                                                                           >  1:1996028/1‑150 (MQ=35)
             cGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGCag                                                                                                                                          <  2:1158368/149‑2 (MQ=11)
             cGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGc                                                                                                                                                   >  2:1327553/1‑140 (MQ=37)
             cGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGc                                                                                                                                                   <  1:1327553/140‑1 (MQ=37)
              gCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATc                                                                                                                                               <  1:2769586/143‑1 (MQ=34)
               ccACTTCATTCCCGTATCAAGCAAAAGGGCCGTTGCGGTGGCGTTCGGTATCCGCGAAGAAACCATCTCCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGc                                                                                                                                       <  1:2646388/150‑1 (MQ=11)
                  cTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCaa                                                                                                                                    >  1:1407984/1‑150 (MQ=35)
                   ttCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAAc                                                                                                                                   >  2:2674905/1‑150 (MQ=35)
                      tttGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCGGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTc                                                                                                                                <  2:2617827/149‑1 (MQ=21)
                      aTTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGACTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTc                                                                                                                                >  1:796934/1‑150 (MQ=25)
                        ggCCGTATCCAGCAACAAGCCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCaa                                                                                                                              <  1:1804538/149‑1 (MQ=25)
                         gCCGTATCCAGCAACAAGGCCGTTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAg                                                                                                                             <  2:833842/150‑1 (MQ=21)
                                 cAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGCAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCAAACTTCAAGGAActgct                                                                                                                     >  1:881517/1‑150 (MQ=25)
                                 cAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgct                                                                                                                     >  2:1389320/1‑150 (MQ=35)
                                 cAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgct                                                                                                                     <  1:1761576/150‑1 (MQ=34)
                                  aGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgctg                                                                                                                    >  1:1038773/1‑150 (MQ=35)
                                   gCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGt                                                                                                                   >  2:1606933/1‑150 (MQ=34)
                                     aacaaGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTcc                                                                                                                 >  1:2116670/1‑150 (MQ=35)
                                     aaaaaGGCCGTTGCTTTAGCGTTCGTTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTcc                                                                                                                 <  2:796934/147‑1 (MQ=255)
                                      acaaGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCg                                                                                                                >  1:601675/1‑150 (MQ=35)
                                        aaGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgc                                                                                                              >  2:1353777/1‑150 (MQ=35)
                                        aaGCCCGCTGCTTTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGCACTGGGTGGGCGGGCACTACTCGGTGTTGTCGGCCATCGGCCTGGCGATCCCCCTTGCTATCTGCACTGCCAACTTCAAGGAACTGCTGTCCGgc                                                                                                              <  1:1241156/150‑1 (MQ=11)
                                         aGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgcg                                                                                                             >  2:343099/1‑150 (MQ=37)
                                            ccGCTGCTGTAGCGTTCGGTATCCGCAAAGAAAACATCTTCCCAATGTGGGACTGGCTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCt                                                                                                          <  2:1159744/150‑1 (MQ=255)
                                                  ctgTGGCGTTCGGTATCCGCAAAGAAAACACCTTCCCATGTTTGGACTGGTTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAGGAACTGCTGTCCGGCGCCTACACCAt                                                                                                   <  2:1567402/150‑1 (MQ=11)
                                                   tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAgg                                                                                                                            >  1:1563173/1‑125 (MQ=37)
                                                   tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAgg                                                                                                                            <  2:1563173/125‑1 (MQ=35)
                                                   tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCAt                                                                                                   <  2:688951/150‑1 (MQ=38)
                                                    gTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTAGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATg                                                                                                  >  1:928181/1‑150 (MQ=38)
                                                          ttCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAg                                                                                            <  1:2737588/150‑1 (MQ=38)
                                                           tCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGc                                                                                           >  1:2763898/1‑150 (MQ=39)
                                                               tATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACtt                                                                                       >  2:2291328/1‑150 (MQ=39)
                                                                  ccGCGAAGAAAACATCTCCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCa                                                                                    <  2:2452760/150‑1 (MQ=38)
                                                                    gcgAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTa                                                                                                         <  1:778883/127‑1 (MQ=39)
                                                                    gcgAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGa                                                                                  >  1:3008969/1‑150 (MQ=39)
                                                                    gcgAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGa                                                                                  >  1:668643/1‑150 (MQ=39)
                                                                    gcgAAAAAAACTTCTTCCAAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCAGGCGCCTACACCATGGACCAGCACTTCCAGa                                                                                  <  1:2098979/150‑1 (MQ=255)
                                                                      gaagaaAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGAcc                                                                                >  1:2814381/1‑150 (MQ=39)
                                                                            aaCATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgc                                                                                                              >  2:848057/1‑114 (MQ=39)
                                                                            aaCATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgc                                                                                                              <  1:848057/114‑1 (MQ=38)
                                                                              cATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                        <  2:3008969/150‑1 (MQ=38)
                                                                              cATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                        <  2:204777/150‑1 (MQ=38)
                                                                               aTCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAg                                                                                            <  2:2798174/129‑1 (MQ=39)
                                                                               aTCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAg                                                                                            >  1:2798174/1‑129 (MQ=40)
                                                                               aTCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTc                                                                       <  1:1311671/150‑1 (MQ=37)
                                                                               aTCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTc                                                                       >  2:1161457/1‑150 (MQ=39)
                                                                                tCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCg                                                                      >  2:3070316/1‑150 (MQ=39)
                                                                                 tttCCCAAGGTGGGACTGGGGGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGa                                                                     <  2:2763898/149‑1 (MQ=255)
                                                                                 cTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGa                                                                     <  2:1733249/150‑1 (MQ=38)
                                                                                   tCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTa                                                                                                         >  2:2825500/1‑112 (MQ=255)
                                                                                   tCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTa                                                                                                         <  1:2825500/112‑1 (MQ=40)
                                                                                   tCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACa                                                                   >  1:2217562/1‑150 (MQ=39)
                                                                                     ccAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAg                                                                 <  2:183145/150‑1 (MQ=38)
                                                                                      cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCt                                                                                                          <  1:1138785/108‑1 (MQ=40)
                                                                                      cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGAc                                                                                               >  2:2322213/1‑119 (MQ=40)
                                                                                      cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                >  1:682065/1‑150 (MQ=38)
                                                                                      cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                <  2:1967722/150‑1 (MQ=38)
                                                                                      cAAGGTGGGACGGGGTGGGCGGGCTCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                <  2:711444/150‑1 (MQ=255)
                                                                                        aTGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAAc                                                              <  2:668643/150‑1 (MQ=35)
                                                                                             ggACTGGGTGGGAGGGCGCTACTCGATGTGGCTCGCCATCGGCCTGCAGATCGCCCTGGCCATCGGCACTGCCAACTACAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGcc                                                         <  2:2173626/150‑1 (MQ=11)
                                                                                                   ggtgggCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATg                                                                                                  <  2:1972186/103‑1 (MQ=38)
                                                                                                         cgggcgCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCTATCTGCTCTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGc                                             >  2:431555/1‑150 (MQ=16)
                                                                                                         cgggcgCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGc                                             <  1:2883773/150‑1 (MQ=33)
                                                                                                             cgcTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCctg                                         >  2:566564/1‑150 (MQ=32)
                                                                                                                 aCTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTGTGGCTCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                                     <  1:2674905/150‑1 (MQ=12)
                                                                                                                   tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGt                                                      <  2:2219230/131‑1 (MQ=33)
                                                                                                                   tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGt                                                      >  1:2219230/1‑131 (MQ=32)
                                                                                                                   tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGc                                   <  1:2291328/150‑1 (MQ=32)
                                                                                                                      cTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtg                                <  1:343099/150‑1 (MQ=32)
                                                                                                                       tgtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                               <  1:1775410/150‑1 (MQ=32)
                                                                                                                        gtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                              <  1:2233133/150‑1 (MQ=32)
                                                                                                                        gtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                              <  2:2231157/150‑1 (MQ=32)
                                                                                                                          ggTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGt                            <  1:566564/150‑1 (MQ=32)
                                                                                                                           gtcgtcCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGTAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCAACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTa                           <  1:1540378/150‑1 (MQ=12)
                                                                                                                            tCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGATCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTAt                          <  2:601675/150‑1 (MQ=17)
                                                                                                                            tCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTAt                          >  1:1066733/1‑150 (MQ=32)
                                                                                                                            tCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTAt                          >  1:953088/1‑150 (MQ=32)
                                                                                                                              ggCGATCGGCAGGCCGATGGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATgg                        >  2:2988877/1‑150 (MQ=255)
                                                                                                                              ggCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCa                                                                                             >  1:2435735/1‑81 (MQ=35)
                                                                                                                              ggCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATgg                        <  1:1353777/150‑1 (MQ=32)
                                                                                                                              ggCCAGCGGCCTGCCGATCGCCCTGGCCATCGGAACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATgg                        <  2:2804772/150‑1 (MQ=17)
                                                                                                                                  ttggcagtGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAAc                    >  1:2056704/8‑150 (MQ=11)
                                                                                                                                     gcagtGCCGATGGCCAGGGCGATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTc                 <  2:2056704/146‑1 (MQ=11)
                                                                                                                                      gCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGAc                                                                                               <  1:2773164/71‑1 (MQ=35)
                                                                                                                                      gCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGAc                                                                                               >  2:2773164/1‑71 (MQ=35)
                                                                                                                                      gCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCt                >  2:1067321/1‑150 (MQ=32)
                                                                                                                                        cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCTGCAAGACCATGCCGCTGCTACCACCACTGATCGTCGTGTGGTATGGCAACTTCTgg              <  1:1606933/150‑1 (MQ=11)
                                                                                                                                             gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                                     >  2:778831/1‑122 (MQ=33)
                                                                                                                                             gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTAATGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGc         >  1:194994/1‑150 (MQ=17)
                                                                                                                                             gATCGCCCTGGCCATCGGCACTGCCAACTTCAACGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGc         <  1:980049/150‑1 (MQ=17)
                                                                                                                                                  cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTGCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGACTATGGCAAGGCAAGTTCTGTGGCCCAAACa    >  2:1183834/1‑150 (MQ=11)
                                                                                                                                                  cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    <  2:1407984/150‑1 (MQ=32)
                                                                                                                                                  cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    <  2:2116670/150‑1 (MQ=32)
                                                                                                                                                  cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    <  1:381888/150‑1 (MQ=32)
                                                                                                                                                     |                                                                                                                                                    
AGCCGAGCTGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGC  >  NC_002947/5343519‑5343816

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