Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 612,761 T→C 5.4% N475S (AAC→AGC)  dxs ← 1‑deoxyxylulose‑5‑phosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947612,7610TC 5.4% 393.1 / 8.1 131N475S (AAC→AGC) dxs1‑deoxyxylulose‑5‑phosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (55/68);  new base C (3/4);  total (59/72)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.46e-01

Only 100 of 139 total aligned reads displayed.

CTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCGGCGCGAT  >  NC_002947/612613‑612908
                                                                                                                                                    |                                                                                                                                                    
cTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                                     >  1:236998/1‑150 (MQ=255)
    gCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTAc                                                                                                                                                 >  2:276912/1‑150 (MQ=255)
    gCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTAc                                                                                                                                                 <  2:695828/150‑1 (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)
               gCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTCCAGGTCGCCCGt                                                                                                                                      <  2:1560995/150‑1 (MQ=255)
                    acacCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCCCATCGGCGCATTCGTTGGTGTTCCACTCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGa                                                                                                                                 <  1:285749/150‑1 (MQ=255)
                        cagcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGa                                                                                                                                 >  1:1436644/1‑146 (MQ=255)
                        cagcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGa                                                                                                                                 <  2:1436644/146‑1 (MQ=255)
                         agcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                                     <  1:1373914/125‑1 (MQ=255)
                         agcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGtt                                                                                                                                                     >  2:1373914/1‑125 (MQ=255)
                         agcagGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                            >  1:1713468/1‑150 (MQ=255)
                           cgatctcACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTCCAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCtt                                                                                                                          <  2:1693073/144‑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)
                                 cACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAACGGCCCGGCCGCTGTa                                                                                                                                                  <  1:2325226/120‑1 (MQ=255)
                                 cACTTTTTCGCCCTGGCCGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCAGTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                            <  1:1056808/142‑1 (MQ=255)
                                  aCTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgt                                                                                                                   >  1:582073/1‑150 (MQ=255)
                                  aCTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGACCGGTGCTGAGCAGCTTGCGCAgt                                                                                                                   >  1:1935471/1‑150 (MQ=255)
                                    tttttCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttc                                                                                                                 >  2:37535/1‑150 (MQ=255)
                                     ttttCGCCCTGGCGGCGGACCACGCCCTTGCCGACCTCCAGTGGTTGCAGGTCCCCGCTGATCAGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcg                                                                                                                <  1:1052977/150‑1 (MQ=255)
                                       ttCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcgtt                                                                                                              <  2:2037168/150‑1 (MQ=255)
                                        tCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgttcgttc                                                                                                             <  2:1233286/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)
                                                  ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgg                                                                                                   <  1:1666315/150‑1 (MQ=255)
                                                  ggcggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgg                                                                                                   <  1:514718/150‑1 (MQ=255)
                                                    cggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                                 <  1:195993/150‑1 (MQ=255)
                                                    cggACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTgggg                                                                                                 <  1:2192561/150‑1 (MQ=255)
                                                        ccACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtca                                                                                             >  2:2130859/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)
                                                                  gCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCAt                                                                                   <  1:261704/150‑1 (MQ=255)
                                                                         tcCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCGGGGCCGTTGTATAGGTGGCCGGTGCTGGGCATCTTGCGCCGTTCGTTCTCGTCGCTGGGCGTCATCGCCAGCATGCCCggg                                                                            <  2:424350/150‑1 (MQ=255)
                                                                          ccAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATc                                                                                                                            <  1:201626/101‑1 (MQ=255)
                                                                           cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTg                                                                                                    <  1:2063478/124‑1 (MQ=255)
                                                                           cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTAGGGGTCATCACCAGCATGCCCGGGAt                                                                          >  1:623236/1‑150 (MQ=255)
                                                                           cAGTGGTTGCAGGTCTCCGCTGATCGGCGCATACGGGCCGATACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTACGAAGTTCGTTCTCGTCGCTGGGGGTCATCAACAGCATGCCCGGGAt                                                                          >  2:444509/1‑150 (MQ=255)
                                                                               ggTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAg                                                                      <  2:2015922/150‑1 (MQ=255)
                                                                                gTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAgc                                                                     >  2:642278/1‑150 (MQ=255)
                                                                                     agttctgCGCTGATCGCCGCATTCGGGCCGGTACCACGCGGGTAGCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTCCGTAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAgg                                                               <  2:1006324/147‑1 (MQ=255)
                                                                                       gTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTgcg                                                                                                                       <  1:446536/93‑1 (MQ=255)
                                                                                            cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAgtt                                                                                                                  <  2:489764/93‑1 (MQ=255)
                                                                                            cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGc                                                                                     <  2:655188/122‑1 (MQ=255)
                                                                                            cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACa                                                        <  2:1904277/150‑1 (MQ=255)
                                                                                            cGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGAGGCCGGGTCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACa                                                        >  1:1118567/1‑150 (MQ=255)
                                                                                             gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAg                                                       >  2:1793113/1‑150 (MQ=255)
                                                                                             gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAg                                                       <  2:1713468/150‑1 (MQ=255)
                                                                                             gCTGATCGGCGCATTCGGGCCGGTACCACGCGGGGACCGTACAGCGGGCGGGCCGGTGAACAGGTGGGCGGTGATGCGCATGTTGCGCAGTACGTTGTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGCGGCAGGCATGGTACGACAg                                                       >  2:721542/1‑150 (MQ=255)
                                                                                               tGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCAt                                                                                   <  1:1753262/121‑1 (MQ=255)
                                                                                                 aTCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGc                                                                                                      <  2:210444/100‑1 (MQ=255)
                                                                                                 aTCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGc                                                                                                      >  1:210444/1‑100 (MQ=255)
                                                                                                   cGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTa                                                 <  2:1930160/150‑1 (MQ=255)
                                                                                                    ggCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGtcatca                                                                                          >  2:255614/1‑109 (MQ=255)
                                                                                                     gcgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGc                                               >  1:1196731/1‑150 (MQ=255)
                                                                                                     gcgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGc                                               >  2:221227/1‑150 (MQ=255)
                                                                                                     gcgcATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGc                                               >  1:975972/1‑150 (MQ=255)
                                                                                                        cATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACGGGTCGTAGCTGc                                            <  2:1149333/150‑1 (MQ=255)
                                                                                                        cATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGc                                            <  2:1447900/150‑1 (MQ=255)
                                                                                                        cATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGc                                            <  1:657432/150‑1 (MQ=255)
                                                                                                           tCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGCTACGACCTGTCGTACCTGcgct                                         <  1:37535/150‑4 (MQ=255)
                                                                                                             gggCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCg                                       >  2:874110/1‑150 (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)
                                                                                                                   gTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCAAGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTc                                 >  1:2007822/1‑150 (MQ=255)
                                                                                                                    tACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTAGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCg                                >  2:2303836/1‑150 (MQ=255)
                                                                                                                      ccACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCggg                              <  2:582073/150‑1 (MQ=255)
                                                                                                                      ccACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCggg                              <  2:1754848/150‑1 (MQ=255)
                                                                                                                         gcgtggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCg                           >  1:1705922/5‑150 (MQ=255)
                                                                                                                         cgcgGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGAAGTTCGTTATCGTCTGTTGGGGTCATCACCAGTATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCg                           <  1:1450003/150‑1 (MQ=255)
                                                                                                                           cgGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTc                         <  1:721542/150‑1 (MQ=255)
                                                                                                                            gggTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTcc                        >  1:1432597/1‑150 (MQ=255)
                                                                                                                               tACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTGGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGt                                  <  1:1384228/137‑1 (MQ=255)
                                                                                                                               acctgTACAACGGCCCGGCCGCTGTACGGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCAc                                                                                         <  2:1628171/79‑1 (MQ=25)
                                                                                                                               acctgTACAACGGCCCGGCCGCTGTACGGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCAc                                                                                         >  1:1628171/5‑83 (MQ=25)
                                                                                                                                aCCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGc                    <  1:1730103/150‑1 (MQ=255)
                                                                                                                                   gTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGcc                   <  1:721600/148‑1 (MQ=255)
                                                                                                                                    tACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACCGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGAc                <  1:2221303/150‑1 (MQ=255)
                                                                                                                                    tACAACGGCCCGGCCGCTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATg                                                                         <  2:680450/94‑1 (MQ=255)
                                                                                                                                    tACAACGGCCCGGCCGCTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATg                                                                         >  1:680450/1‑94 (MQ=255)
                                                                                                                                       aGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAg             <  1:419284/150‑1 (MQ=255)
                                                                                                                                           gccgggccgTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCg         <  1:664149/150‑1 (MQ=255)
                                                                                                                                                gccgTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCGGcgcg    >  1:217355/1‑150 (MQ=255)
                                                                                                                                                  cgTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCTGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCGGCGcgat  >  1:1661636/1‑150 (MQ=255)
                                                                                                                                                   gTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTc                         >  1:890884/1‑126 (MQ=255)
                                                                                                                                                    ttGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCGGCGcga   >  2:802216/1‑147 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTCGGCCAGTTGTACGCCGAACACCAGCAGGGCCACTTTTTCGCCCTGGCGGCGGACCACGCCCTTGCCGATCTCCAGTGGTTGCAGGTCTCCGCTGATCGGCGCATTCGGGCCGGTACCACGCGGGTACCGTACAGCGGCCGGGCCGTTGTACAGGTGGCCGGTGCTGAGCATCTTGCGCAGTTCGTTCTCGTCGCTGGGGGTCATCACCAGCATGCCCGGGATGCAGCGCAGGTACGACAGGTCGTAGCTGCCGGCGTGGGTCGGGCCGTCCTCGCCGACCAGGCCGGCGCGAT  >  NC_002947/612613‑612908

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