Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,343,672 C→G 5.9% A288A (GCC→GCG 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,6720CG 5.9% 350.6 / 12.1 121A288A (GCC→GCGpgi‑IIglucose 6‑phosphate isomerase
Reads supporting (aligned to +/- strand):  ref base C (50/62);  new base G (2/5);  total (53/68)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.90e-01

Only 100 of 123 total aligned reads displayed.

TGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGC  >  NC_002947/5343527‑5343816
                                                                                                                                                 |                                                                                                                                                
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)
     cGCCACTTCATTGCCGTATCAAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATGGCCAGGGCGATCGGCag                                                                                                                                          >  1:1158368/1‑148 (MQ=11)
      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)
                tGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCaa                                                                                                                              <  1:2778715/150‑1 (MQ=34)
                 gCCGTATCCAGCAACAAGGCCGTTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAg                                                                                                                             <  2:833842/150‑1 (MQ=21)
                       tCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActg                                                                                                                       >  1:1733249/1‑150 (MQ=35)
                         cAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGCAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCAAACTTCAAGGAActgct                                                                                                                     >  1:881517/1‑150 (MQ=25)
                         cAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgct                                                                                                                     <  1:1761576/150‑1 (MQ=34)
                          aGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAActgctg                                                                                                                    >  1:1038773/1‑150 (MQ=35)
                             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)
                                 aGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgcg                                                                                                             >  2:343099/1‑150 (MQ=37)
                                          ctgTGGCGTTCGGTATCCGCAAAGAAAACACCTTCCCATGTTTGGACTGGTTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAGGAACTGCTGTCCGGCGCCTACACCAt                                                                                                   <  2:1567402/150‑1 (MQ=11)
                                           tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAgg                                                                                                                            <  2:1563173/125‑1 (MQ=35)
                                           tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAgg                                                                                                                            >  1:1563173/1‑125 (MQ=37)
                                           tgTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCAt                                                                                                   <  2:688951/150‑1 (MQ=38)
                                                  ttCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAg                                                                                            <  1:911560/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)
                                                            gcgAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGa                                                                                  >  1:668643/1‑150 (MQ=39)
                                                              gaagaaAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGAcc                                                                                >  1:2814381/1‑150 (MQ=39)
                                                                 gaaAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGcc                                                                             >  2:1775410/1‑149 (MQ=39)
                                                                 gaaAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCa                                                                             >  2:1540378/1‑150 (MQ=37)
                                                                    aaCATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgc                                                                                                              >  2:848057/1‑114 (MQ=39)
                                                                    aaCATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGgc                                                                                                              <  1:848057/114‑1 (MQ=38)
                                                                      cATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                        <  2:204777/150‑1 (MQ=38)
                                                                      cATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGtt                                                                        <  2:3008969/150‑1 (MQ=38)
                                                                       aTCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAg                                                                                            <  2:2798174/129‑1 (MQ=39)
                                                                       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                                                                                                         <  1:2825500/112‑1 (MQ=40)
                                                                           tCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTa                                                                                                         >  2:2825500/1‑112 (MQ=255)
                                                                           tCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACa                                                                   >  1:2217562/1‑150 (MQ=39)
                                                                             ccAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAg                                                                 <  2:183145/150‑1 (MQ=38)
                                                                              cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCt                                                                                                          >  2:1138785/1‑108 (MQ=40)
                                                                              cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCt                                                                                                          <  1:1138785/108‑1 (MQ=40)
                                                                              cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGAc                                                                                               <  1:2322213/119‑1 (MQ=39)
                                                                              cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                >  1:682065/1‑150 (MQ=38)
                                                                              cAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGa                                                                >  1:2978387/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)
                                                                                      gACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCg                                                        >  1:555575/1‑150 (MQ=35)
                                                                                           ggtgggCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATg                                                                                                  >  1:1972186/1‑103 (MQ=38)
                                                                                           ggtgggCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCAGCACTGCCAACTTCATTGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCt                                                   >  1:2675513/1‑150 (MQ=255)
                                                                                                 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                                                      >  1:2219230/1‑131 (MQ=32)
                                                                                                           tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGt                                                      <  2:2219230/131‑1 (MQ=33)
                                                                                                           tCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGc                                   <  2:2217562/150‑1 (MQ=32)
                                                                                                              cTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtg                                <  1:343099/150‑1 (MQ=32)
                                                                                                               tgtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                               <  2:1242011/150‑1 (MQ=32)
                                                                                                               tgtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCgtgt                               <  1:1775410/150‑1 (MQ=32)
                                                                                                                gtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                              <  2:2231157/150‑1 (MQ=32)
                                                                                                                gtgGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGtgtg                              <  1:2233133/150‑1 (MQ=32)
                                                                                                                    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                                                                                             <  2:2435735/81‑1 (MQ=35)
                                                                                                                      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                                                                                               >  2:2773164/1‑71 (MQ=35)
                                                                                                                              gCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGAc                                                                                               <  1:2773164/71‑1 (MQ=35)
                                                                                                                              gCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCt                >  2:1067321/1‑150 (MQ=32)
                                                                                                                                cTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCTGCAAGACCATGCCGCTGCTACCACCACTGATCGTCGTGTGGTATGGCAACTTCTgg              <  1:1606933/150‑1 (MQ=11)
                                                                                                                                     gATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTgg                                     <  1:778831/122‑1 (MQ=33)
                                                                                                                                     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:2116670/150‑1 (MQ=32)
                                                                                                                                          cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    <  1:381888/150‑1 (MQ=32)
                                                                                                                                          cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    <  2:1407984/150‑1 (MQ=32)
                                                                                                                                          cccTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACa    <  2:1996028/150‑1 (MQ=32)
                                                                                                                                            cTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGc  <  1:1389320/150‑1 (MQ=32)
                                                                                                                                                 cATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTgg                                              <  1:538772/101‑1 (MQ=33)
                                                                                                                                                 cATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTgg                                              >  2:538772/1‑101 (MQ=33)
                                                                                                                                                 |                                                                                                                                                
TGTACCGCCACTTCATTGCCGTATCCAGCAACAAGGCCGCTGCTGTAGCGTTCGGTATCCGCGAAGAAAACATCTTCCCAATGTGGGACTGGGTGGGCGGGCGCTACTCGCTGTGGTCGGCCATCGGCCTGCCGATCGCCCTGGCCATCGGCACTGCCAACTTCAAGGAACTGCTGTCCGGCGCCTACACCATGGACCAGCACTTCCAGACCGCACCGTTCGACAAGAACATGCCGGTGCTACTGGCCCTGCTGGGCGTGTGGTATGGCAACTTCTGGGACGCAAACAGC  >  NC_002947/5343527‑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: