Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,030,478 C→T 8.5% intergenic (+435/‑259) pcaY → / → mgtA aromatic acid chemoreceptor/ATP‑dependent magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,030,4780CT 8.5% 261.1 / 12.0 94intergenic (+435/‑259)pcaY/mgtAaromatic acid chemoreceptor/ATP‑dependent magnesium transporter
Reads supporting (aligned to +/- strand):  ref base C (40/46);  new base T (4/4);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.53e-01

TTGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAACCATGGATTGCGC  >  NC_002947/3030331‑3030625
                                                                                                                                                   |                                                                                                                                                   
ttGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGgcgc                                                                                                                                                   <  2:815537/150‑1 (MQ=255)
ttGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGgcgc                                                                                                                                                   <  2:1786929/150‑1 (MQ=255)
 tGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGcgcg                                                                                                                                                  <  2:1542050/150‑1 (MQ=255)
     ggCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCtt                                                                                                                                              >  1:1813865/1‑150 (MQ=255)
        aaGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                                           <  1:339080/150‑1 (MQ=255)
        aaGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                                           >  2:1426750/1‑150 (MQ=255)
        aaGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                                           >  1:1530541/1‑150 (MQ=255)
           cTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCa                                                                                                                                        >  1:277563/1‑150 (MQ=255)
            tCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCAc                                                                                                                                       >  1:893867/1‑150 (MQ=255)
            tCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCAc                                                                                                                                       <  2:1530541/150‑1 (MQ=255)
             cATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACg                                                                                                                                      >  1:675149/1‑150 (MQ=255)
                  aCCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGa                                                                                                                                 <  2:799100/150‑1 (MQ=255)
                    cTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTACACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGc                                                                                                                               >  1:308008/1‑150 (MQ=255)
                     tCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGcc                                                                                                                              >  2:821591/1‑150 (MQ=255)
                     tCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTCCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGcc                                                                                                                              >  2:494334/1‑150 (MQ=255)
                         ttCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCg                                                                                                                          >  1:1201887/1‑150 (MQ=255)
                            aTTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCAc                                                                                                                       >  1:1361303/1‑150 (MQ=255)
                             ttCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACt                                                                                                                      >  2:1521512/1‑150 (MQ=255)
                              tCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACtt                                                                                                                     <  2:675149/150‑1 (MQ=255)
                                gTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGACAGCCCTCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTtct                                                                                                                   <  2:51477/150‑1 (MQ=255)
                                 tttCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                                           >  2:121141/1‑125 (MQ=255)
                                 tttCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCAt                                                                                                                                           <  1:121141/125‑1 (MQ=255)
                                   tCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCAc                                                                                                                >  1:1325209/1‑150 (MQ=255)
                                     acaTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTg                                                                                                              <  1:1253017/150‑1 (MQ=255)
                                       aTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGTACGCTTGAGCCCCCGCACTTCTCACTGaa                                                                                                            >  1:372882/1‑150 (MQ=255)
                                       aTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTtc                                                                                                                    <  2:343645/142‑1 (MQ=255)
                                       aTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTtc                                                                                                                    >  1:343645/1‑142 (MQ=255)
                                       aTCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTgga                                                                                                            >  1:966928/1‑148 (MQ=255)
                                        tCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAAt                                                                                                           <  1:1426750/150‑1 (MQ=255)
                                           ccAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaa                                                                                                        >  2:1291788/1‑150 (MQ=255)
                                            cAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                                       <  1:402296/150‑1 (MQ=255)
                                            cAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                                       <  1:1291788/150‑1 (MQ=255)
                                                  gCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCATAAACCAAGGATTTGCGTGGGGCTTGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATAtt                                                                                                 >  1:1669406/1‑150 (MQ=255)
                                                  gCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATAtt                                                                                                 <  1:1521512/150‑1 (MQ=255)
                                                     tGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGaa                                                                                              <  1:277777/150‑1 (MQ=255)
                                                     tGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGaa                                                                                              <  1:821591/150‑1 (MQ=255)
                                                       ccTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACg                                                                                            >  2:1583858/1‑150 (MQ=255)
                                                        cTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGTGGTGCTGGCGCGTGTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACgg                                                                                           >  2:1038192/1‑150 (MQ=255)
                                                            cgccgcGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                                       <  1:755416/134‑1 (MQ=255)
                                                            cgccgcGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                                       >  2:755416/1‑134 (MQ=255)
                                                               gcggccAGGCGCGGGCCTGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTg                                                                                   >  2:372777/6‑150 (MQ=255)
                                                                  ggCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTa                                                                                 <  1:436561/150‑1 (MQ=255)
                                                                   gCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTAc                                                                                <  2:34233/150‑1 (MQ=255)
                                                                     aGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTc                                                                              <  1:1410931/150‑1 (MQ=255)
                                                                        cccGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATAtt                                                                                                 <  1:1233862/128‑1 (MQ=255)
                                                                        cccGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATAtt                                                                                                 >  2:1233862/1‑128 (MQ=255)
                                                                        cccGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGaa                                                                           <  2:1480725/150‑1 (MQ=255)
                                                                            cgcCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACttt                                                                       <  2:1813865/150‑1 (MQ=255)
                                                                               cTGGCCGCGGCGTCTGCACGTCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCa                                                                    <  1:1918523/150‑1 (MQ=255)
                                                                                tGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGATACTGTCTCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCAt                                                                   <  2:58343/150‑1 (MQ=255)
                                                                                   ccGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGtt                                                                >  2:1530236/1‑150 (MQ=255)
                                                                                      cggcggCTGCACGGCTACCTTGACCACCACCAGCACCCCCCCAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGTATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTaaa                                                             <  2:852642/150‑1 (MQ=255)
                                                                                      cggcgTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTaaa                                                             <  1:985923/150‑1 (MQ=255)
                                                                                           tCTGCACGTCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTc                                                        <  1:407921/150‑1 (MQ=255)
                                                                                           tCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTTATTTGACAAGTTc                                                        >  2:1544740/1‑150 (MQ=255)
                                                                                             tGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTcgc                                                      >  2:1839941/1‑150 (MQ=255)
                                                                                             tGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTcgc                                                      <  1:494334/150‑1 (MQ=255)
                                                                                               cACGGCTACCTTGACCACCACCAGCCCCCCCCCAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGACCGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCgg                                                    <  2:592139/150‑1 (MQ=255)
                                                                                                 cGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGTTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGa                                                  <  2:617855/150‑1 (MQ=255)
                                                                                                  ggCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGt                                                                                  >  1:843785/1‑117 (MQ=255)
                                                                                                  ggCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGNGAt                                                 <  2:1201887/150‑1 (MQ=255)
                                                                                                  ggCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGAt                                                 >  1:392978/1‑150 (MQ=255)
                                                                                                  ggCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTAGCTTTGt                                                                                  <  2:843785/117‑1 (MQ=255)
                                                                                                     tACCTTGCCCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAAc                                              <  1:1475910/150‑1 (MQ=255)
                                                                                                       ccTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTa                                            >  2:34551/1‑150 (MQ=255)
                                                                                                            accaccaccAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTGCTGTAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTcc                                       >  1:404738/1‑150 (MQ=255)
                                                                                                             ccaccaccaGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCgg                                                    >  2:59445/1‑136 (MQ=255)
                                                                                                             ccaccaccaGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCgg                                                    <  1:59445/136‑1 (MQ=255)
                                                                                                              caccaccaGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCt                                     <  1:1839941/150‑1 (MQ=255)
                                                                                                                  gccaGCACCCCGCAACTACTTGGGGGGTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACt                                                                                                               >  1:1329349/2‑72 (MQ=18)
                                                                                                                  gccaGCACCCCGCAACTACTTGGGGGGTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACt                                                                                                               <  2:1329349/71‑1 (MQ=18)
                                                                                                                     aGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTa                              <  2:277563/150‑1 (MQ=255)
                                                                                                                      gCACCCCGCAACTACTTGGGGGGTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTAt                             >  1:1289473/1‑150 (MQ=255)
                                                                                                                      gCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTAt                             <  2:392978/150‑1 (MQ=255)
                                                                                                                      gCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTAt                             <  2:893867/150‑1 (MQ=255)
                                                                                                                        aCCCCGCAACTACTTGGGGGGTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGa                           >  2:822556/1‑150 (MQ=255)
                                                                                                                        aCCCCGCAACTACTTGGGGGGTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGa                           >  2:822776/1‑150 (MQ=255)
                                                                                                                        aCCCCGCAACTACTTGGGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGa                           >  1:812678/1‑150 (MQ=255)
                                                                                                                        aCCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGa                           >  1:691371/1‑150 (MQ=255)
                                                                                                                           ccccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAAGTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGaaaa                        >  2:1002751/1‑150 (MQ=255)
                                                                                                                           ccccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGTTCTTTTTCTATTTCCCTGTGTTTATGaaaa                        <  1:990215/150‑1 (MQ=255)
                                                                                                                           ccccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGaaaa                        >  1:834873/1‑150 (MQ=255)
                                                                                                                            cccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCttt                                                                                    >  1:1640766/1‑89 (MQ=255)
                                                                                                                            cccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCttt                                                                                    <  2:1640766/89‑1 (MQ=255)
                                                                                                                            cccAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAg                       >  2:1640848/1‑150 (MQ=255)
                                                                                                                                    gtTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCaaa               <  2:1121215/150‑1 (MQ=255)
                                                                                                                                    gtTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTAGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCaaa               <  1:481789/150‑1 (MQ=255)
                                                                                                                                    gtTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTAACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCaaa               <  1:1544740/150‑1 (MQ=255)
                                                                                                                                    cttggggggTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCaaa               <  1:822556/145‑1 (MQ=255)
                                                                                                                                    cttggggggTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCaaa               <  1:822776/145‑1 (MQ=255)
                                                                                                                                      tGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAAcc             <  2:1029080/150‑1 (MQ=255)
                                                                                                                                       ggggggTGCTGGTGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAACCa            <  2:1289473/148‑1 (MQ=255)
                                                                                                                                         ggggTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAACCATg          >  1:1244550/1‑150 (MQ=255)
                                                                                                                                              gCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAACCATGGATTg     <  1:1196628/150‑1 (MQ=255)
                                                                                                                                                tGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                                       >  2:1860749/1‑50 (MQ=255)
                                                                                                                                                tGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATaaaa                                                                                                       <  1:1860749/50‑1 (MQ=255)
                                                                                                                                                 ggCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAACCATGGATTgcgc  <  1:372777/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTGGAGGCAAGCTCATGAACCTCACTTCATTCGTTTCACATCGCCAGCCTGCGTGCCTTGCGCCGCGGCAGGCCCGCGCCTGGCCGCGGCGTCTGCACGGCTACCTTGACCACCACCAGCACCCCCCAAGTAGTTGCGGGGTGCTGGCGCGTCTTCATGCACGCTTGAGCCCCCGCACTTCTCACTGAATAAAAGATATTGAACGGTATCTTTGTACTCGAACTTTTCATGTTAAAAGTTCGCGGGATAACTATTTCCCTGTGTTTATGAAAAGGAATCAAACCATGGATTGCGC  >  NC_002947/3030331‑3030625

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: