Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,337,829 G→C 6.9% S330S (TCG→TCC) ‡ PP_2934 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,337,8290GC 6.9% 398.9 / 15.0 126S330S (TCG→TCC) ‡PP_2934hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (59/58);  new base C (5/4);  total (64/62)
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.42e-01

Only 100 of 127 total aligned reads displayed.

GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACC  >  NC_002947/3337680‑3337968
                                                                                                                                                     |                                                                                                                                            
gggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCg                                                                                                                                              >  2:160559/1‑150 (MQ=255)
 ggCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGc                                                                                                                                             <  2:3414/150‑1 (MQ=255)
  gCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGAAGCCAGACACGACAGCGTTCTTCAGCCAGACGCTGGGGTCGcc                                                                                                                                            >  2:1123934/1‑150 (MQ=255)
   cATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCa                                                                                                                                           <  1:359834/150‑1 (MQ=255)
     tttCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgc                                                                                                                                         >  1:72298/1‑150 (MQ=255)
       tCGAATGCATTGCCCATTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAgcgc                                                                                                                                       <  2:262484/150‑1 (MQ=255)
        cGTATTCATGGCCCAAGGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGNAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGCGACCCCAGCGcc                                                                                                                                      <  1:1442420/150‑1 (MQ=255)
         gTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCC‑GGCGCTGGCGACCCCAGCGcctg                                                                                                                                    <  1:534673/150‑3 (MQ=255)
          tATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑gg                                                                                                                                   >  1:1521198/1‑150 (MQ=255)
            ttCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCt                                                                                                                                 >  1:2029577/1‑150 (MQ=255)
              cATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGa                                                                                                                               <  2:3089/150‑1 (MQ=255)
                  gCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGt                                                                                                                           >  2:1318422/1‑150 (MQ=255)
                    ccACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAATGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTc                                                                                                                         >  2:1069523/1‑150 (MQ=255)
                     cACTGTCATGCCGAGCAGTACGCGCGAATGCCGGCGCTGGGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGTGCGTTCa                                                                                                                        <  2:1344389/150‑1 (MQ=255)
                     cACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                        <  1:743448/150‑1 (MQ=255)
                      aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑g                                                                                                                                    <  2:1599488/137‑1 (MQ=255)
                      aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑g                                                                                                                                    >  1:1599488/1‑137 (MQ=255)
                      aCTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAg                                                                                                                       <  1:1123934/150‑1 (MQ=255)
                       cTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGc                                                                                                                      <  2:2036110/150‑1 (MQ=255)
                         gCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTg                                                                                                                    <  1:228855/150‑1 (MQ=255)
                             tGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGAc                                                                                                                >  2:1817228/1‑150 (MQ=255)
                                 gAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGgc                                                                                                            >  1:1089931/1‑150 (MQ=255)
                                       tatgcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGGC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccg                                                                                                      >  1:953352/4‑150 (MQ=255)
                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                        >  2:43016/1‑129 (MQ=255)
                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCa                                                                                                                        <  1:43016/129‑1 (MQ=255)
                                          gcgcgAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   >  2:899390/1‑150 (MQ=255)
                                              gAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccagc                                                                                               >  1:649909/1‑150 (MQ=255)
                                                 tGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                            >  2:222648/1‑150 (MQ=255)
                                                   ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCTATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                          <  1:2068151/150‑1 (MQ=255)
                                                   ccGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACaa                                                                                          >  2:1503000/1‑150 (MQ=255)
                                                    cGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                         >  2:1630088/1‑150 (MQ=255)
                                                     ggCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACg                                                                                        >  1:915163/1‑150 (MQ=255)
                                                          tgtgCAAGGACGCCGGCGGAGACCAGCCGCTGTATATCCATCATTAGCTACAGCGCGAACCCACAGCTGGCTGCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                   <  1:44472/150‑1 (MQ=255)
                                                                tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGAcc                                                                             >  2:2325104/3‑150 (MQ=255)
                                                                tcgACGCCGTCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGACGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACGACGCTAGGCCGAcc                                                                             >  2:1126018/3‑150 (MQ=255)
                                                                       cGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAAAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACGACGTTAGGCCGACCCCGGCCa                                                                      >  2:1239473/1‑150 (MQ=255)
                                                                       cGACGGCGTCGACCGCCGGTACATACACGAATCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCa                                                                      >  1:65048/1‑150 (MQ=255)
                                                                        gACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCaa                                                                     >  1:2032291/1‑150 (MQ=255)
                                                                         aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCggactac                                                                    <  1:1069523/150‑6 (MQ=255)
                                                                         aCGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                    <  2:654210/150‑1 (MQ=255)
                                                                              gTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAg                                                               >  2:1867665/1‑150 (MQ=255)
                                                                               tCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGACAAGGCGAGc                                                              >  1:1844670/1‑150 (MQ=255)
                                                                                  aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCaaa                                                           <  2:1844670/150‑3 (MQ=255)
                                                                                  aCCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAgg                                                           <  2:2032291/150‑1 (MQ=255)
                                                                                   ccgccgGTATATCCATACTTCGCTGCAGCGAGAACCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTGCCCAGCGGG‑GGCTGGGCGTACACCCGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAggg                                                          <  2:2407781/150‑1 (MQ=255)
                                                                                    cgccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGc                                                         >  2:989937/1‑150 (MQ=255)
                                                                                      ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAg                                                                                    <  1:1324228/121‑1 (MQ=255)
                                                                                      ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAg                                                                                    >  2:1324228/1‑121 (MQ=255)
                                                                                      ccgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAg                                                       >  2:1698564/1‑150 (MQ=255)
                                                                                       cgGTATCTCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGc                                                      <  1:1565677/150‑1 (MQ=255)
                                                                                       cgGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGc                                                      <  2:915163/150‑1 (MQ=255)
                                                                                        ggTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGcc                                                     <  2:836347/150‑1 (MQ=255)
                                                                                          tataTCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCAc                                                   <  2:410280/150‑1 (MQ=255)
                                                                                             atCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                <  2:1521198/150‑1 (MQ=255)
                                                                                             atCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                >  1:404464/1‑150 (MQ=255)
                                                                                              tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCt                                                                                                                     >  2:1382383/1‑80 (MQ=255)
                                                                                              tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCt                                                                                                                     <  1:1382383/80‑1 (MQ=255)
                                                                                              tCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCag                                               >  2:409727/1‑150 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAAc                                                                                         >  2:76482/1‑106 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                             >  2:908780/1‑150 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                             <  1:1838017/150‑1 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCagaa                                             >  1:2286183/1‑150 (MQ=255)
                                                                                                catcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCa                                                <  1:1631759/147‑1 (MQ=255)
                                                                                                 atcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAgg                                                                                   <  2:908418/111‑1 (MQ=255)
                                                                                                  tcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                >  1:588132/1‑97 (MQ=255)
                                                                                                  tcatTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCgccag                                                                                                <  2:588132/97‑1 (MQ=255)
                                                                                                    atTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATc                                         <  1:1698564/150‑1 (MQ=255)
                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                             >  1:830791/1‑81 (MQ=255)
                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACggg                                                                                                             <  2:830791/81‑1 (MQ=255)
                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                    <  2:1836175/122‑1 (MQ=255)
                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAg                                                                    >  1:1836175/1‑122 (MQ=255)
                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAg                                                            <  2:1864116/130‑1 (MQ=255)
                                                                                                     ttCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAg                                                            >  1:1864116/1‑130 (MQ=255)
                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  1:899390/150‑1 (MQ=255)
                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  1:156242/150‑1 (MQ=255)
                                                                                                       cGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCACTTGAGACGTGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAg                                      <  2:83296/150‑1 (MQ=255)
                                                                                                             aGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAc                                <  1:660959/150‑1 (MQ=255)
                                                                                                              tgtagaAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGTCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGCCAGCCACGCAGAAGATCCAGGCATAcc                               >  2:1150834/4‑150 (MQ=255)
                                                                                                              gCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcc                               <  2:1511829/150‑1 (MQ=255)
                                                                                                                  gaAACGGCAGCGTTCTTCAGCCATGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAg                           <  1:2089365/150‑1 (MQ=255)
                                                                                                                   aaaCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCT‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAg                           <  2:2188945/149‑1 (MQ=255)
                                                                                                                     aCGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCCGCAACAACGTTCACCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGgtg                        >  1:286449/1‑150 (MQ=255)
                                                                                                                      cGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGcc                                                                                                                                      >  2:655008/1‑40 (MQ=255)
                                                                                                                       ggCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCcac                                                                                            <  1:1218018/80‑1 (MQ=255)
                                                                                                                         cAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                    <  2:1905953/150‑1 (MQ=255)
                                                                                                                          aGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCg                   >  1:142830/1‑150 (MQ=255)
                                                                                                                           tgaacgctcAGCC‑GGCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGc                  >  2:1976105/8‑150 (MQ=37)
                                                                                                                           gCGTTCTTCAGCCAGGCGATGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGc                  <  1:1817228/150‑1 (MQ=255)
                                                                                                                             gTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCa                >  1:1783284/1‑150 (MQ=255)
                                                                                                                                gctcAGCC‑GGCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACgg             >  1:296251/3‑150 (MQ=38)
                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   >  2:2186076/1‑59 (MQ=255)
                                                                                                                                     gCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGcgcgccgcgc                                                                                                   <  1:2186076/59‑1 (MQ=255)
                                                                                                                                      ccAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa       >  1:809130/1‑150 (MQ=255)
                                                                                                                                       ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAgg                                     >  1:1875595/3‑120 (MQ=38)
                                                                                                                                       ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa       >  2:1703486/3‑150 (MQ=255)
                                                                                                                                       ccggCGCTGGCGACCCCAGCGCCTGGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCa       >  2:1703363/3‑150 (MQ=255)
                                                                                                                                       cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                             >  1:94999/1‑127 (MQ=255)
                                                                                                                                       cAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATAcccc                             <  2:94999/127‑1 (MQ=255)
                                                                                                                                                   tCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                    >  1:886259/1‑124 (MQ=255)
                                                                                                                                                   tCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTAcc                    <  2:886259/124‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GGGCATTTCGTATTCATGGCCCACTGCGATGCCGAGCAGTACGCGCGAATGCCGGCGCTGTGCAAGGACGCCGACGGCGTCGACCGCCGGTATATCCATCATTCGCTGCAGCGAGAAACGGCAGCGTTCTTCAGCCAGGCGCTGGGGTCGCCAGCGCC‑GGCTGAGCGTTCAGCTGCGACGGGCGCGCCGCGCCAGCCACAACGTTAGGCCGACCCCGGCCAAGGCGAGCAGGGCAGCCACGCAGAAGATCGAGGCATACCCCAGGTGTACCGCCACGGCGCCCATCACC  >  NC_002947/3337680‑3337968

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