Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,257,515 T→A 5.8% T161S (ACC→TCC) ‡ PP_2854 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,257,5150TA 5.8% 287.0 / 6.7 100T161S (ACC→TCC) ‡PP_2854hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (51/43);  new base A (2/4);  total (53/47)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.70e-01

AAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGC  >  NC_002947/3257369‑3257661
                                                                                                                                                  |                                                                                                                                                  
aaGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTg                                                                                                                                                 <  1:2431687/150‑1 (MQ=255)
       gTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                          <  1:2416200/150‑1 (MQ=255)
        tCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcgg                                                                                                                                         >  1:959295/1‑150 (MQ=255)
        tCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcgg                                                                                                                                         >  2:363667/1‑150 (MQ=255)
         cGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCggcggc                                                                                                                                        <  1:1174451/150‑1 (MQ=255)
         cGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCggcggc                                                                                                                                        <  1:1683471/150‑1 (MQ=255)
         cGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCggcggc                                                                                                                                        <  2:2176413/150‑1 (MQ=255)
          gAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCACGGTCTGGCCGGCGGcc                                                                                                                                       >  2:1224840/1‑150 (MQ=255)
            gTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCac                                                                                                                                     >  1:527291/1‑149 (MQ=255)
                tGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAgaccc                                                                                                                                 >  1:70888/1‑146 (MQ=255)
                 gcagcaGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCt                                                                                                                                <  2:1750190/150‑1 (MQ=255)
                  cagcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAgaccctg                                                                                                                               <  2:70888/150‑7 (MQ=255)
                   agcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTgg                                                                                                                              >  1:733112/1‑150 (MQ=255)
                       ggCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCtt                                                                                                                          >  1:678315/1‑150 (MQ=255)
                       ggCAACCTTGGGCCAGGCGTTGAACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGATGGAAGACACACTGTGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTggtcgg                                                                                                                          <  2:58448/150‑5 (MQ=255)
                        gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAgaccctg                                                                                                                               >  2:1404147/1‑138 (MQ=255)
                        gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAgaccctg                                                                                                                               <  1:1404147/144‑7 (MQ=255)
                        gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTggtcggg                                                                                                                         <  2:1174348/150‑6 (MQ=255)
                        gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTg                                                                                                                         >  2:2122444/1‑150 (MQ=255)
                         ccACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGc                                                                                                                        <  2:1040475/150‑1 (MQ=255)
                          cACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCACCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCt                                                                                                                       >  1:1980449/1‑150 (MQ=255)
                             cTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                    >  2:2389177/1‑150 (MQ=255)
                             cTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                    >  2:167544/1‑150 (MQ=255)
                             cTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                    >  2:2408358/1‑150 (MQ=255)
                             cTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                    >  2:1646815/1‑150 (MQ=255)
                             cTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                                    >  1:1396986/1‑150 (MQ=255)
                              ttGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgttg                                                                                                                   >  1:1056489/1‑150 (MQ=255)
                              ttGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgttg                                                                                                                   <  2:959295/150‑1 (MQ=255)
                                     aGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTAcggcg                                                                                                            >  1:1147744/1‑150 (MQ=255)
                                       gCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCggg                                                                                                          >  2:266990/1‑150 (MQ=255)
                                         gTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTg                                                                                                        >  2:770828/1‑150 (MQ=255)
                                                gctgctGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGtt                                                                                                 <  1:266990/150‑1 (MQ=255)
                                                 ctgctgCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTg                                                                                                >  1:1718336/1‑150 (MQ=255)
                                                  tgctgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                               >  1:333471/1‑150 (MQ=255)
                                                   gctgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGc                                                                                              <  2:2449354/150‑1 (MQ=255)
                                                    ctgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCt                                                                                             <  1:1282415/150‑1 (MQ=255)
                                                     tgcCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCtt                                                                                            >  1:2221821/1‑150 (MQ=255)
                                                      gccatGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCAACGGCACGGCGGGTGCCTGGTTGGCttt                                                                                           <  1:1520107/146‑1 (MQ=255)
                                                       cccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACCCCGCTGGAAGACACCCTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                          <  2:1890651/150‑1 (MQ=255)
                                                       cccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                          >  2:419249/1‑150 (MQ=255)
                                                         cTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGt                                                                                        >  1:675077/1‑150 (MQ=255)
                                                          tGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTc                                                                                       >  1:1096779/1‑150 (MQ=255)
                                                          tGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTc                                                                                       >  1:226381/1‑150 (MQ=255)
                                                             ttgttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACg                                                                                    <  2:333471/150‑1 (MQ=255)
                                                               gttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGtt                                                                                  <  2:226381/150‑1 (MQ=255)
                                                                ttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTg                                                                                 <  1:419249/150‑1 (MQ=255)
                                                                 tgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGa                                                                                <  2:1396986/150‑1 (MQ=255)
                                                                    tGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGatat                                                                             >  2:131489/1‑150 (MQ=255)
                                                                     ggCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCGTTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATc                                                                            <  2:1091498/150‑1 (MQ=255)
                                                                      gCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                               <  2:2068588/130‑1 (MQ=255)
                                                                      gCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTgg                                                                                               >  1:2068588/1‑130 (MQ=255)
                                                                       cATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGa                                                                          <  2:1286591/150‑1 (MQ=255)
                                                                        aTTGATGTTGAGCTGCACACCCAGACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGATCTgg                                                                                                                                                >  2:1190364/1‑79 (MQ=255)
                                                                        aTTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTgg                                                                                                                                                <  1:1190364/79‑1 (MQ=255)
                                                                            aTGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttg                                                                     <  2:1980449/150‑1 (MQ=255)
                                                                             tGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgt                                                                    <  2:675077/150‑1 (MQ=255)
                                                                                  aGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTaa                                                               <  1:1611511/150‑1 (MQ=255)
                                                                                  aGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTaa                                                               <  1:167544/150‑1 (MQ=255)
                                                                                     tGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcg                                                            >  1:486484/1‑150 (MQ=255)
                                                                                     tGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCACCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcg                                                            >  1:486725/1‑150 (MQ=255)
                                                                                      gCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGttgttg                                                                  >  2:612176/1‑143 (MQ=255)
                                                                                      gCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGttgttg                                                                  <  1:612176/143‑1 (MQ=255)
                                                                                          acccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCa                                                                                      <  2:461355/119‑1 (MQ=255)
                                                                                          acccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCa                                                                                      >  1:461355/1‑119 (MQ=255)
                                                                                          acccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgtt                                                       >  1:1618000/1‑150 (MQ=255)
                                                                                            ccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgttgt                                                     >  2:1313171/1‑150 (MQ=255)
                                                                                            ccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgttgt                                                     <  2:678315/150‑1 (MQ=255)
                                                                                             cacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgttgtt                                                    <  1:2131949/150‑1 (MQ=255)
                                                                                                 ccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCa                                                >  1:1216143/1‑150 (MQ=255)
                                                                                                 ccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCa                                                >  1:34516/1‑150 (MQ=255)
                                                                                                  cGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCACACATGGATTACCTGCGCAGACATTGTTGCGGGGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCAc                                               >  1:1485966/1‑150 (MQ=255)
                                                                                                      ggAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcggcg                                                         >  1:1509324/1‑136 (MQ=255)
                                                                                                      ggAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcggcg                                                         <  2:1509324/136‑1 (MQ=255)
                                                                                                      ggAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGc                                           >  1:1513353/1‑150 (MQ=255)
                                                                                                       gAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCt                                          >  2:820464/1‑150 (MQ=255)
                                                                                                           acacacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCAc                                      >  2:1679951/1‑150 (MQ=255)
                                                                                                           acacacTCAGCGAGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCAc                                      >  2:1685584/1‑150 (MQ=255)
                                                                                                             acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCCGCTGCTGCACgg                                    >  1:958083/1‑150 (MQ=255)
                                                                                                             acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACgg                                    <  2:733112/150‑1 (MQ=255)
                                                                                                              cacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                   >  1:1245420/1‑118 (MQ=255)
                                                                                                              cacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                                   <  2:1245420/118‑1 (MQ=255)
                                                                                                               acTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                          >  1:339012/1‑94 (MQ=255)
                                                                                                               acTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                          <  2:339012/94‑1 (MQ=255)
                                                                                                               acTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGac                                  >  1:187312/1‑150 (MQ=255)
                                                                                                                 tGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGTCAATGTTCTGGTTGGCTTTGGTCACGTTGATATGGACTTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGAACGGagca                                >  1:1118064/1‑150 (MQ=255)
                                                                                                                      gcgcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTgg                           <  2:2068550/150‑1 (MQ=255)
                                                                                                                             cgcCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                    <  2:2221821/150‑1 (MQ=255)
                                                                                                                              gcCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACAcc                   >  1:1403936/1‑150 (MQ=255)
                                                                                                                                cATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                    <  1:1310071/147‑1 (MQ=255)
                                                                                                                                cATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTacac                    >  2:1310071/1‑147 (MQ=255)
                                                                                                                                   tttCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCgg              <  2:2323863/150‑1 (MQ=255)
                                                                                                                                       cccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTTATCGTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggt          <  1:2122444/150‑1 (MQ=255)
                                                                                                                                        ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGTCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtg         <  2:527291/150‑1 (MQ=255)
                                                                                                                                        ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtg         <  1:131489/150‑1 (MQ=255)
                                                                                                                                        ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGAAGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtg         >  2:2160660/1‑150 (MQ=255)
                                                                                                                                          gACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtggt       >  1:706884/1‑150 (MQ=255)
                                                                                                                                            ccAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCa     <  1:1679951/150‑1 (MQ=255)
                                                                                                                                            ccAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCa     <  1:1685584/150‑1 (MQ=255)
                                                                                                                                             cAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAg    >  2:1651052/1‑150 (MQ=255)
                                                                                                                                               gggTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGGGGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGc  <  1:874159/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTGGTCAGGC  >  NC_002947/3257369‑3257661

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