Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,160,961 C→G 5.5% intergenic (+320/+548) PP_5739 → / ← PP_5402 hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,160,9610CG 5.5% 667.5 / 20.8 200intergenic (+320/+548)PP_5739/PP_5402hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit
Reads supporting (aligned to +/- strand):  ref base C (95/94);  new base G (4/7);  total (99/101)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 206 total aligned reads displayed.

ACCAAGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTCCG  >  NC_002947/6160812‑6161110
                                                                                                                                                     |                                                                                                                                                     
   aaGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCa                                                                                                                                                    <  1:4044015/150‑1 (MQ=255)
    aGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCaa                                                                                                                                                   <  1:5035672/150‑1 (MQ=255)
      ccGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAgg                                                                                                                                                 <  2:1436678/150‑1 (MQ=255)
       cGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAggc                                                                                                                                                <  1:36552/150‑1 (MQ=255)
        gCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAggcg                                                                                                                                               >  1:4724863/1‑150 (MQ=255)
            ttGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCtt                                                                                                                                           <  2:3274043/149‑1 (MQ=255)
               aTCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttt                                                                                                                                        <  2:1512486/150‑1 (MQ=255)
                tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttta                                                                                                                                       <  2:92377/150‑2 (MQ=255)
                tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttta                                                                                                                                       <  2:5186107/150‑2 (MQ=255)
                tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttta                                                                                                                                       <  2:1555729/150‑2 (MQ=255)
                tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttt                                                                                                                                       >  2:2600882/1‑150 (MQ=255)
                tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttt                                                                                                                                       >  2:2292275/1‑150 (MQ=255)
                tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttt                                                                                                                                       >  2:1694707/1‑150 (MQ=255)
                 ccATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCttttttt                                                                                                                                      >  2:3020882/1‑150 (MQ=255)
                 ccATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCttttttt                                                                                                                                      >  1:700455/1‑150 (MQ=255)
                 ccATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCttttttt                                                                                                                                      <  1:3496627/150‑1 (MQ=255)
                  cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttttt                                                                                                                                     <  1:1469455/150‑1 (MQ=255)
                  cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATGTCGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAAGTCGTTGGCTAAAAAGCCGCCCTCTCACGGCGGCtttttttt                                                                                                                                     <  2:1779856/150‑1 (MQ=255)
                    tttCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGc                                                                                                                                   >  2:5060965/1‑150 (MQ=255)
                        gCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCg                                                                                                                               >  1:485923/1‑150 (MQ=255)
                         cTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCgg                                                                                                                              >  2:1488317/1‑150 (MQ=255)
                          tttAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGa                                                                                                                             <  2:3222890/150‑1 (MQ=255)
                          tttAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGa                                                                                                                             >  1:4066807/1‑150 (MQ=255)
                             aGCGGCTAAATTTTTCGTTTGTTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTc                                                                                                                          <  1:651795/150‑1 (MQ=255)
                             aGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTc                                                                                                                          <  2:1490905/150‑1 (MQ=255)
                              gCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCg                                                                                                                         >  1:821040/1‑150 (MQ=255)
                               cGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGt                                                                                                                        <  2:4639743/150‑1 (MQ=255)
                                ggCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTa                                                                                                                       >  2:3314056/1‑150 (MQ=255)
                                  cTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTAGAGTCCTCACCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAAGTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATc                                                                                                                     <  2:1210054/150‑1 (MQ=255)
                                  cTAAATTTGTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATc                                                                                                                     <  2:4594967/150‑1 (MQ=255)
                                   tAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCg                                                                                                                    <  2:1872316/150‑1 (MQ=255)
                                   tAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCg                                                                                                                    >  1:3423360/1‑150 (MQ=255)
                                     aaTTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGt                                                                                                                  >  1:1250603/1‑150 (MQ=255)
                                             gTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACcgtcgt                                                                                                          <  1:2577338/150‑1 (MQ=255)
                                              tttGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGtg                                                                                                         >  2:1967688/1‑150 (MQ=255)
                                                  gTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGtgctgc                                                                                                     >  2:4649611/1‑150 (MQ=255)
                                                   ttGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGcc                                                                                                    <  1:3679094/150‑1 (MQ=255)
                                                    tGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCa                                                                                                   <  2:2108733/150‑1 (MQ=255)
                                                       gTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAgtg                                                                                                >  1:4329150/1‑150 (MQ=255)
                                                          attaattaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCg                                                                                             >  2:1872698/1‑150 (MQ=255)
                                                           ttaattaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGc                                                                                            >  2:1838739/1‑150 (MQ=255)
                                                           ttaattaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGc                                                                                            >  2:2741470/1‑150 (MQ=255)
                                                              attaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGa                                                                                         >  1:1166992/1‑150 (MQ=255)
                                                                    aTTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGa                                                                                   >  1:3127891/1‑150 (MQ=255)
                                                                    aTTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGa                                                                                   <  1:5443275/150‑1 (MQ=255)
                                                                      ttAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATc                                                                                 >  2:3342199/1‑150 (MQ=255)
                                                                      ttAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATc                                                                                 <  1:4971561/150‑1 (MQ=255)
                                                                         gCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACt                                                                              <  1:3342199/150‑1 (MQ=255)
                                                                           tataGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACtgt                                                                            <  2:4392223/150‑1 (MQ=255)
                                                                             taGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                          <  2:112877/150‑1 (MQ=255)
                                                                             taGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                          <  1:1071661/150‑1 (MQ=255)
                                                                             taGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                          <  1:1920218/150‑1 (MQ=255)
                                                                              atctcTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTATTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAg                                                                         <  1:1626295/146‑1 (MQ=255)
                                                                                tCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGa                                                                       >  2:2946477/1‑150 (MQ=255)
                                                                                  cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                     >  1:2943617/1‑150 (MQ=255)
                                                                                  cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                     >  1:3836279/1‑150 (MQ=255)
                                                                                  cTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAcc                                                                     >  2:5220052/1‑150 (MQ=255)
                                                                                   tCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCa                                                                    >  2:4508851/1‑150 (MQ=255)
                                                                                     aaCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                          <  1:2445847/142‑1 (MQ=255)
                                                                                         ttCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTAt                                                              <  2:4461370/150‑1 (MQ=255)
                                                                                            gACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGtgctgc                                                                                                     <  1:2220328/108‑1 (MQ=255)
                                                                                            gACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCAt                                                           <  2:3810603/150‑1 (MQ=255)
                                                                                              cTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGc                                                         <  1:5179089/150‑1 (MQ=255)
                                                                                                gTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGt                                                       <  2:3423360/150‑1 (MQ=255)
                                                                                                 tttCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTa                                                      >  1:3442938/1‑150 (MQ=255)
                                                                                                  ttCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTac                                                     <  2:5427938/150‑1 (MQ=255)
                                                                                                             gCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAaccacc                                          <  1:1329842/150‑1 (MQ=255)
                                                                                                              ccGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCt                                         <  1:5061107/150‑1 (MQ=255)
                                                                                                               cgacgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtg                                        >  1:2110991/1‑150 (MQ=255)
                                                                                                                gacgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtgt                                       <  2:4502902/150‑1 (MQ=255)
                                                                                                                gacgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtgt                                       <  1:4185158/150‑1 (MQ=255)
                                                                                                                gacgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTATTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtgt                                       <  1:3992409/150‑1 (MQ=255)
                                                                                                                  cgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGa                                     <  2:4061105/150‑1 (MQ=255)
                                                                                                                    aaCCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACt                                   <  1:5355395/150‑1 (MQ=255)
                                                                                                                       ctctCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGa                                                                                         >  1:1700962/1‑93 (MQ=255)
                                                                                                                        tctcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTCTCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTc                               >  2:3497894/1‑150 (MQ=255)
                                                                                                                        tctcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTc                               >  1:2615851/1‑150 (MQ=255)
                                                                                                                        tctcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGCACACTCAACCACCTGTGACTGTTc                               >  1:1306508/1‑150 (MQ=255)
                                                                                                                         ctcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCg                              <  2:5094879/150‑1 (MQ=255)
                                                                                                                         ctcAATTCGTTGGCTAAAAAGCCGCCCCCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCCGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCg                              <  2:2366053/150‑1 (MQ=255)
                                                                                                                          tcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGa                                                                                   >  2:807705/1‑96 (MQ=255)
                                                                                                                          tcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGa                                                                                   <  1:807705/96‑1 (MQ=255)
                                                                                                                           cAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGc                            <  1:1488799/150‑1 (MQ=255)
                                                                                                                            aaTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGcc                           >  1:593179/1‑150 (MQ=255)
                                                                                                                             aTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGACTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCt                          >  2:4933669/1‑150 (MQ=255)
                                                                                                                                cGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGtt                       <  1:5423586/150‑1 (MQ=255)
                                                                                                                                 gTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTg                      <  1:303176/150‑1 (MQ=255)
                                                                                                                                     gCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCaa                  <  2:4066807/150‑1 (MQ=255)
                                                                                                                                        aaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGa                                     <  1:3875642/128‑1 (MQ=255)
                                                                                                                                        aaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGa                                     >  2:3875642/1‑128 (MQ=255)
                                                                                                                                        aaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCa               >  1:2752108/1‑150 (MQ=255)
                                                                                                                                         aaaaGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGa                                                                          >  2:4903538/1‑90 (MQ=255)
                                                                                                                                         aaaaGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCAc              >  2:1593945/1‑150 (MQ=255)
                                                                                                                                              ccgccCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGc                    >  1:100534/1‑139 (MQ=255)
                                                                                                                                              ccgccCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGc                    <  2:100534/139‑1 (MQ=255)
                                                                                                                                                gccCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAgtc       <  1:5196419/150‑1 (MQ=255)
                                                                                                                                                    tctcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTa                                                      <  2:2123078/99‑1 (MQ=255)
                                                                                                                                                    tctcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTg                      >  1:4117661/1‑131 (MQ=255)
                                                                                                                                                    tctcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACGCAACCACCTGTGACTGTTCGGCCTGTTg                      <  2:4117661/131‑1 (MQ=255)
                                                                                                                                                     ctcAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTCCg  >  2:4210465/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCAAGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTCGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCAGTCGTCCG  >  NC_002947/6160812‑6161110

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