Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 612,755 G→C 5.3% P477R (CCG→CGG)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947612,7550GC 5.3% 434.0 / 11.9 132P477R (CCG→CGG) dxs1‑deoxyxylulose‑5‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base G (56/69);  new base C (3/4);  total (59/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

Only 100 of 138 total aligned reads displayed.

GGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCG  >  NC_002947/612610‑612901
                                                                                                                                                 |                                                                                                                                                  
ggCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccgggcc                                                                                                                                                >  1:2057697/1‑150 (MQ=255)
 gCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGCTCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccgggccg                                                                                                                                               >  2:563616/1‑150 (MQ=255)
 gCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTATTGCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccgggccg                                                                                                                                               >  1:2290401/1‑150 (MQ=255)
   cTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                             <  1:1984066/150‑1 (MQ=255)
   cTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                             >  1:236998/1‑150 (MQ=255)
       gCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTAc                                                                                                                                         <  2:695828/150‑1 (MQ=255)
       gCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTAc                                                                                                                                         >  2:276912/1‑150 (MQ=255)
         cAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                             >  1:2142477/1‑144 (MQ=255)
         cAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                             <  2:2142477/144‑1 (MQ=255)
         cAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAg                                                                                                                                       <  1:557512/150‑1 (MQ=255)
           gTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCTCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAACGGCCCGGCCGCTGTAcggta                                                                                                                                     <  1:1684067/150‑5 (MQ=255)
                aCGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCg                                                                                                                                <  1:15300/150‑1 (MQ=255)
                  gCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGt                                                                                                                              >  2:1580737/1‑150 (MQ=255)
                       acacCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCCCATCGGCGCATTCGTTGGTGTTCCACTCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGa                                                                                                                         <  1:285749/150‑1 (MQ=255)
                           cagcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGa                                                                                                                         <  2:1436644/146‑1 (MQ=255)
                           cagcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGa                                                                                                                         >  1:1436644/1‑146 (MQ=255)
                            agcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                             <  1:1373914/125‑1 (MQ=255)
                            agcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                    >  1:1713468/1‑150 (MQ=255)
                            agcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTACAGTGATCACACGGGTGCACGGATCGGCGAATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                    <  1:756085/150‑1 (MQ=255)
                              cgatctcACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTCCAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCtt                                                                                                                  <  2:1693073/144‑1 (MQ=255)
                               agGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCCGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTg                                                                                                                 <  1:2172884/150‑1 (MQ=255)
                                gggCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTgc                                                                                                                <  1:1337762/150‑1 (MQ=255)
                                   ccACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTgcgcc                                                                                                             >  2:569514/1‑149 (MQ=255)
                                   ccACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCa                                                                                                             >  2:2192561/1‑150 (MQ=255)
                                    cACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                    >  2:1056808/1‑142 (MQ=255)
                                    cACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAACGGCCCGGCCGCTGTa                                                                                                                                          >  2:2325226/1‑120 (MQ=255)
                                    cACTTTTTCGCCCTGGCCGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCAGTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                    <  1:1056808/142‑1 (MQ=255)
                                     aCTTTTTCGCCCTGGCGGCGGACCACGCCCTTTCCGAGCTCCGATGGTTGCGAGTCTCCGCGGAGCGGCGCGTTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCAGTTGTACAGGTGGCCGGTGCTGCGCATCTTgagatcg                                                                                                           >  1:1693073/1‑144 (MQ=255)
                                      cTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgtt                                                                                                          <  2:632951/150‑1 (MQ=255)
                                       tttttCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttc                                                                                                         <  1:1722590/150‑1 (MQ=255)
                                        ttttCGCCCTGGCGGCGGACCACGCCCTTGCCGACCTCCAGTGGTTGCAGGTCCCCGCTGATCAGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcg                                                                                                        <  1:1052977/150‑1 (MQ=255)
                                          ttCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTTGACTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcgtt                                                                                                      <  1:1874254/150‑1 (MQ=255)
                                          ttCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcgtt                                                                                                      <  2:2037168/150‑1 (MQ=255)
                                          ttCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGACGTTGCCACGCGGGTATTGGATAGCTGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcgtt                                                                                                      <  1:2036865/150‑1 (MQ=255)
                                           tCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcgttc                                                                                                     <  2:1233286/150‑1 (MQ=255)
                                             gCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTgctc                                                                                                   <  1:563616/150‑1 (MQ=255)
                                              cccTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcg                                                                                                  >  2:2078650/1‑150 (MQ=255)
                                               ccTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcgt                                                                                                 >  1:1904277/1‑150 (MQ=255)
                                                 tGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCtcgtcg                                                                                               >  2:1424752/1‑150 (MQ=255)
                                                    cggcggACCACGCCCTTGCCGATCTCTAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTg                                                                                            <  2:2057697/150‑1 (MQ=255)
                                                     ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCggccgg                                                                                                                                                   >  2:767927/1‑94 (MQ=255)
                                                     ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgg                                                                                           >  2:2221303/1‑150 (MQ=255)
                                                     ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgg                                                                                           <  1:675633/150‑1 (MQ=255)
                                                       cggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                         <  1:195993/150‑1 (MQ=255)
                                                           ccACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtca                                                                                     >  2:1043552/1‑150 (MQ=255)
                                                             aCGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtcatc                                                                                   <  1:648178/150‑1 (MQ=255)
                                                                 ccTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCa                                                                               <  1:2078650/150‑1 (MQ=255)
                                                                  ctggCCGATTTCCAGTGGTTGCAGGCCTCCGCTCATCCGCGCATCCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTNCTGAGCNTCTTGCGCAGTTNGTTCTCGTCGCTGGGGGTCATCACCAg                                                                              <  1:2151883/147‑1 (MQ=255)
                                                                  cTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAg                                                                              <  1:78502/150‑1 (MQ=255)
                                                                  cTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAg                                                                              <  1:791519/150‑1 (MQ=255)
                                                                   ttGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCa                                                                                                                      <  2:1522367/109‑1 (MQ=255)
                                                                            tcCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCGGGGCCGTTGTATAGGTGGCCGGTGCTGGGCATCTTGCGCCGTTCGTTCTCGTCGCTGGGCGTCATCGCCAGCATGCCCggg                                                                    <  2:424350/150‑1 (MQ=255)
                                                                             ccAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                    <  1:201626/101‑1 (MQ=255)
                                                                             ccAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                    >  2:201626/1‑101 (MQ=255)
                                                                             ccAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGa                                                                   <  1:40165/150‑1 (MQ=255)
                                                                              cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTAGGGGTCATCACCAGCATGCCCGGGAt                                                                  >  1:623236/1‑150 (MQ=255)
                                                                              cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATACGGGCCGATACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTACGAAGTTCGTTCTCGTCGCTGGGGGTCATCAACAGCATGCCCGGGAt                                                                  >  2:444509/1‑150 (MQ=255)
                                                                                   gTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAgc                                                             >  2:642278/1‑150 (MQ=255)
                                                                                          gTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTgcg                                                                                                               <  1:446536/93‑1 (MQ=255)
                                                                                             tcCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTTGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGa                                                   >  2:1452406/1‑150 (MQ=255)
                                                                                               cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgtt                                                                                                          <  2:489764/93‑1 (MQ=255)
                                                                                               cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgtt                                                                                                          >  1:489764/1‑93 (MQ=255)
                                                                                               cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGc                                                                             >  1:655188/1‑122 (MQ=255)
                                                                                               cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACa                                                 <  2:1904277/150‑1 (MQ=255)
                                                                                               cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGAGGCCGGGTCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACa                                                 >  1:1118567/1‑150 (MQ=255)
                                                                                                gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAg                                                >  1:1636911/1‑150 (MQ=255)
                                                                                                gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAg                                                >  2:1793113/1‑150 (MQ=255)
                                                                                                gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAg                                                >  2:119425/1‑150 (MQ=255)
                                                                                                  tGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCAt                                                                           >  2:1753262/1‑121 (MQ=255)
                                                                                                    aTCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGc                                                                                              >  1:210444/1‑100 (MQ=255)
                                                                                                    aTCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGc                                                                                              <  2:210444/100‑1 (MQ=255)
                                                                                                      cGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTa                                          <  2:1930160/150‑1 (MQ=255)
                                                                                                       ggCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtcatca                                                                                  >  2:255614/1‑109 (MQ=255)
                                                                                                       ggCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtcatca                                                                                  <  1:255614/109‑1 (MQ=255)
                                                                                                        gcgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGc                                        >  1:975972/1‑150 (MQ=255)
                                                                                                        gcgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGc                                        >  1:1196731/1‑150 (MQ=255)
                                                                                                              tCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGCTACGACCTGTCGTACCTGcgct                                  <  1:37535/150‑4 (MQ=255)
                                                                                                                     ggTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                         >  1:1343349/1‑88 (MQ=255)
                                                                                                                     ggTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                         <  2:1343349/88‑1 (MQ=255)
                                                                                                                     ggTACCACGCGGGTACCGTACAGAGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGt                           >  1:468688/1‑150 (MQ=255)
                                                                                                                      gTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTc                          <  1:1793113/150‑1 (MQ=255)
                                                                                                                       tACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTAGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCg                         >  2:2303836/1‑150 (MQ=255)
                                                                                                                         ccACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCggg                       <  2:1754848/150‑1 (MQ=255)
                                                                                                                          cACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGc                      >  2:419284/1‑150 (MQ=255)
                                                                                                                           aCGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGcc                     <  1:2303836/150‑1 (MQ=255)
                                                                                                                            cgcgGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGAAGTTCGTTATCGTCTGTTGGGGTCATCACCAGTATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCg                    <  1:1450003/150‑1 (MQ=255)
                                                                                                                              cgGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTc                  <  1:721542/150‑1 (MQ=255)
                                                                                                                               gggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTcc                 >  1:1432597/1‑150 (MQ=255)
                                                                                                                                  tACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTGGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGt                           >  2:1384228/1‑137 (MQ=255)
                                                                                                                                  tACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTGGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGt                           <  1:1384228/137‑1 (MQ=255)
                                                                                                                                  tACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCg              >  2:664149/1‑150 (MQ=255)
                                                                                                                                  acctgTACAACGGCCCGGCCGCTGTACGGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCAc                                                                                 >  1:1628171/5‑83 (MQ=25)
                                                                                                                                  acctgTACAACGGCCCGGCCGCTGTACGGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCAc                                                                                 <  2:1628171/79‑1 (MQ=25)
                                                                                                                                   aCCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGc             <  1:1730103/150‑1 (MQ=255)
                                                                                                                                      gTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGcc            <  1:721600/148‑1 (MQ=255)
                                                                                                                                      gTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGcc            >  2:721600/1‑148 (MQ=255)
                                                                                                                                       tACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACCGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGAc         <  1:2221303/150‑1 (MQ=255)
                                                                                                                                       tACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGAc         >  2:214224/1‑150 (MQ=255)
                                                                                                                                       tACAACGGCCCGGCCGCTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATg                                                                 <  2:680450/94‑1 (MQ=255)
                                                                                                                                              gccgggccgTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCg  <  1:664149/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GGCCTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCG  >  NC_002947/612610‑612901

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: