Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 542,084 G→C 5.3% L107L (CTG→CTC rpoC → DNA‑directed RNA polymerase subunit beta'

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947542,0840GC 5.3% 694.8 / 20.7 207L107L (CTG→CTCrpoCDNA‑directed RNA polymerase subunit beta'
Reads supporting (aligned to +/- strand):  ref base G (101/95);  new base C (6/5);  total (107/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 209 total aligned reads displayed.

TGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGG  >  NC_002947/541935‑542232
                                                                                                                                                     |                                                                                                                                                    
tGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTg                                                                                                                                                      <  1:4322093/150‑1 (MQ=255)
   gcCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGcc                                                                                                                                                   >  1:1907620/1‑150 (MQ=255)
     cAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTc                                                                                                                                                 >  1:3289865/1‑150 (MQ=255)
       aGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGc                                                                                                                                               <  2:3280279/150‑1 (MQ=255)
        gATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGcc                                                                                                                                              <  1:2380258/150‑1 (MQ=255)
        gATCTTGGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGACGTGAGCGCATGGCGCACATCGAGCTGGCCTCGcc                                                                                                                                              <  2:3039400/150‑1 (MQ=255)
         aTCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCg                                                                                                                                             >  2:3030511/1‑150 (MQ=255)
          tCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCgg                                                                                                                                            >  2:2493568/1‑150 (MQ=255)
          tCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCgg                                                                                                                                            <  1:2043379/150‑1 (MQ=255)
            tttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGtt                                                                                                                                          <  2:547713/150‑1 (MQ=255)
             ttGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTg                                                                                                                                         <  2:2566159/150‑1 (MQ=255)
                gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGccc                                                                                                                                      >  2:4203598/1‑150 (MQ=255)
                gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGccc                                                                                                                                      >  1:3027314/1‑150 (MQ=255)
                gCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGc                                                                                                                                        <  2:3426912/148‑1 (MQ=255)
                  ccAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCcac                                                                                                                                    >  2:1489356/1‑150 (MQ=255)
                    aGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACAt                                                                                                                                  <  2:1406721/150‑1 (MQ=255)
                      tCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCt                                                                                                                                >  1:1819052/1‑150 (MQ=255)
                       cAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTg                                                                                                                               <  2:1839402/150‑1 (MQ=255)
                        aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                              >  1:4271731/1‑150 (MQ=255)
                        aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                              >  1:1409322/1‑150 (MQ=255)
                        aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                              >  1:1409342/1‑150 (MQ=255)
                        aaGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTgg                                                                                                                              >  2:2326480/1‑150 (MQ=255)
                           gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                           >  1:992125/1‑150 (MQ=255)
                           gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                           <  2:4785147/150‑1 (MQ=255)
                           gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                           <  2:5258731/150‑1 (MQ=255)
                           gACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                           <  2:1907620/150‑1 (MQ=255)
                            aCTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTcc                                                                                                                          <  2:1322332/150‑1 (MQ=255)
                              tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                                                                                                        >  1:262886/1‑150 (MQ=255)
                              tACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTg                                                                                                                        >  1:953812/1‑150 (MQ=255)
                                cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                                                                                      >  2:1267421/1‑150 (MQ=255)
                                cGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                                                                                      <  1:4829244/150‑1 (MQ=255)
                                 gAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAg                                                                                                                     >  2:3646930/1‑150 (MQ=255)
                                     gCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGc                                                                                                                 <  1:3247122/150‑1 (MQ=255)
                                       cTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTg                                                                                                               <  1:4406910/150‑1 (MQ=255)
                                       cTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTg                                                                                                               <  1:1246145/150‑1 (MQ=255)
                                        tgtgCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGc                                                                                                              <  2:3772530/150‑1 (MQ=255)
                                            cGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTc                                                                                                          <  2:2068474/150‑1 (MQ=255)
                                             ggCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTcc                                                                                                         <  1:4921872/150‑1 (MQ=255)
                                             ggCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCAGCCGTcc                                                                                                         >  2:3165712/1‑150 (MQ=255)
                                              gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCATATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                                                                                        >  1:1268839/1‑150 (MQ=255)
                                              gCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTccc                                                                                                        >  1:4915232/1‑150 (MQ=255)
                                                aaGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGt                                                                                                      >  2:4714640/1‑150 (MQ=255)
                                                   tACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATc                                                                                                   >  2:836671/1‑150 (MQ=255)
                                                     cAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGaa                                                                                                                      <  2:424465/129‑1 (MQ=255)
                                                       aGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCt                                                                                               >  1:1686735/1‑150 (MQ=255)
                                                           ccTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCt                                                                                           >  1:4189230/1‑150 (MQ=255)
                                                             tCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGa                                                                                         >  2:3923379/1‑150 (MQ=255)
                                                               aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                                                                                       <  2:5360591/150‑1 (MQ=255)
                                                               aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                                                                                       >  2:975218/1‑150 (MQ=255)
                                                               aaGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATg                                                                                       >  2:3390688/1‑150 (MQ=255)
                                                                  cACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGAc                                                                                    <  2:4045201/150‑1 (MQ=255)
                                                                   aCCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGAct                                                                                   <  2:1819052/150‑2 (MQ=255)
                                                                          cgTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCt                                                                            >  1:4694145/1‑150 (MQ=255)
                                                                              aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATcgc                                                                                                        >  2:2128623/1‑116 (MQ=255)
                                                                              aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATcgc                                                                                                        <  1:2128623/118‑3 (MQ=255)
                                                                              aTCTGCGAAAAGTGCGGCGTTGAAGTCACGCTGTCGCGCGTTCGTCGTGAGCGCATGGGCCATATTGAGCTCGCAGCACCTGTTGCCCATATCTGGTTCTTGAAGTCGCTGCCATCGCGCATCGGCACGCTGCTGGACATGACGCtcaag                                                                        <  1:3035398/150‑5 (MQ=255)
                                                                               tCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATAGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTg                                                                       >  2:59152/1‑150 (MQ=255)
                                                                                  gCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGata                                                                    <  1:583762/150‑1 (MQ=255)
                                                                                        aGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAgc                                                              <  2:1594088/150‑1 (MQ=255)
                                                                                          tGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAgcgc                                                            <  2:3342372/150‑1 (MQ=255)
                                                                                          tGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAgcgc                                                            >  1:681082/1‑150 (MQ=255)
                                                                                                gTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGctc                                                      >  1:1491510/1‑150 (MQ=255)
                                                                                                gTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGctc                                                      >  1:1491468/1‑150 (MQ=255)
                                                                                                 ttGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                           >  1:971628/1‑80 (MQ=255)
                                                                                                 ttGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTc                                                                                                                           <  2:971628/80‑1 (MQ=255)
                                                                                                        tGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGctctcgaagta                                              <  2:453381/150‑8 (MQ=255)
                                                                                                          cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                            <  1:2272304/150‑1 (MQ=255)
                                                                                                          cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                            <  2:440131/150‑1 (MQ=255)
                                                                                                          cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                            <  2:439960/150‑1 (MQ=255)
                                                                                                          cccTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCgaga                                            <  2:2467225/150‑1 (MQ=255)
                                                                                                             tGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCt                                         <  1:1929259/150‑1 (MQ=255)
                                                                                                             tGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCATGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCt                                         <  1:1154073/150‑1 (MQ=255)
                                                                                                                 gAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGt                                     >  1:3818962/1‑150 (MQ=255)
                                                                                                                     gTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGtt                                 >  2:3189262/1‑150 (MQ=255)
                                                                                                                     gTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGtt                                 <  1:1320561/150‑1 (MQ=255)
                                                                                                                      ttCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTa                                <  1:4510155/150‑1 (MQ=255)
                                                                                                                       tcgtcgTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTAt                               <  1:2053098/150‑1 (MQ=255)
                                                                                                                       tcgtcgTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTAt                               <  1:2316846/150‑1 (MQ=255)
                                                                                                                        cgtcgtGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATc                              <  2:4950282/150‑1 (MQ=255)
                                                                                                                          tcgtGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGa                            >  1:2326517/1‑150 (MQ=255)
                                                                                                                               aGCGCATGGCGCACATAGAGCAGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCgg                       >  2:3652446/1‑150 (MQ=255)
                                                                                                                                gcgcATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATc                              <  1:1229363/142‑1 (MQ=255)
                                                                                                                                gcgcATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATc                              >  2:1229363/1‑142 (MQ=255)
                                                                                                                                 cgcATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                                                                                                <  1:2397890/75‑1 (MQ=255)
                                                                                                                                 cgcATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGc                                                                                                >  2:2397890/1‑75 (MQ=255)
                                                                                                                                    aTGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATg                  >  2:5035546/1‑150 (MQ=255)
                                                                                                                                     tGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGa                 <  1:5276080/150‑1 (MQ=255)
                                                                                                                                         gcACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTAc             >  1:838572/1‑150 (MQ=255)
                                                                                                                                          cacaTCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTAcc            <  2:2651492/150‑1 (MQ=255)
                                                                                                                                          cacaTCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTAcc            >  2:4140705/1‑150 (MQ=255)
                                                                                                                                              tCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTgg        <  1:1370938/150‑1 (MQ=255)
                                                                                                                                              tCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTgg        <  1:857443/150‑1 (MQ=255)
                                                                                                                                               cGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGa       <  1:3646930/150‑1 (MQ=255)
                                                                                                                                                 aGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGaaa     >  2:1024780/1‑150 (MQ=255)
                                                                                                                                                    tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg  <  2:982902/150‑1 (MQ=255)
                                                                                                                                                    tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg  <  1:5046678/150‑1 (MQ=255)
                                                                                                                                                    tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg  <  1:3956552/150‑1 (MQ=255)
                                                                                                                                                    tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg  <  1:3469049/150‑1 (MQ=255)
                                                                                                                                                    tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg  >  1:2999007/1‑150 (MQ=255)
                                                                                                                                                    tGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAgg  <  1:1198517/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGCGCCAAGATCTTTGGCCCAGTCAAGGACTACGAGTGCCTGTGCGGCAAGTACAAGCGCCTCAAGCACCGCGGCGTAATCTGCGAGAAGTGCGGCGTTGAAGTTGCCCTGGCGAAGGTTCGTCGTGAGCGCATGGCGCACATCGAGCTGGCCTCGCCGGTTGCCCACATCTGGTTCCTGAAGTCGCTGCCGTCCCGTATCGGCTTGCTGATGGACATGACCCTGCGTGATATCGAGCGCGTGCTCTACTTCGAGAGCTATGTCGTTATCGACCCGGGCATGACTACCCTGGAAAAGG  >  NC_002947/541935‑542232

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: