Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,275,090 G→C 8.8% intergenic (‑281/+31) PP_5462 ← / ← PP_5463 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,275,0900GC 8.8% 267.8 / 16.3 91intergenic (‑281/+31)PP_5462/PP_5463hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (43/40);  new base C (4/4);  total (47/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

TCCAGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGT  >  NC_002947/1274941‑1275235
                                                                                                                                                     |                                                                                                                                                  
tCCAGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCccgccg                                                                                                                                                    <  2:1386698/150‑1 (MQ=255)
tCCAGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCccgccg                                                                                                                                                    <  2:654788/150‑1 (MQ=255)
  cAGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGTTTGCTTCGCGGATACGTCGGCTCCAATAAGTTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGAAAAGGTGAAAGGCGGAGGAGCATAGCCCGCCGAt                                                                                                                                                  >  2:408998/1‑150 (MQ=255)
   aGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTTCAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATg                                                                                                                                                 <  1:408998/150‑1 (MQ=255)
   aGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATg                                                                                                                                                 <  1:210163/150‑1 (MQ=255)
        gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGGTGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgg                                                                                                                                            >  2:1829097/1‑150 (MQ=255)
        gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgg                                                                                                                                            >  2:1286086/1‑150 (MQ=255)
        gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgg                                                                                                                                            >  1:1862531/1‑150 (MQ=255)
        gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgg                                                                                                                                            >  2:1426713/1‑150 (MQ=255)
        gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgg                                                                                                                                                >  1:645023/1‑146 (MQ=255)
        gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgg                                                                                                                                                <  2:645023/146‑1 (MQ=255)
          aTGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgggcgggc                                                                                                                                          <  1:481490/150‑1 (MQ=255)
           tGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATCGGCGGGCt                                                                                                                                         <  2:1603019/150‑1 (MQ=255)
               cATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                    <  1:1629000/150‑1 (MQ=255)
                aTTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGAAGGAGGAGGAGTTAAAAGGAGAAGGAGCAAAGCCCGCCGATGGGCGGGCTTCAgg                                                                                                                                   >  1:1868593/1‑150 (MQ=255)
                   ttGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATCGGCGGGCtttgctcat                                                                                                                                 <  1:1743488/150‑8 (MQ=255)
                    tGGTGTATCCGGCGGAGAGTGAATGTCGAAGGCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAgggggg                                                                                                                               <  2:952496/150‑1 (MQ=255)
                    tGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAgggggg                                                                                                                               <  2:540477/150‑1 (MQ=255)
                      gTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTACGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGGGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAgggggggg                                                                                                                             >  1:346699/1‑150 (MQ=255)
                         tATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAg                                                                                                                          >  1:699007/1‑150 (MQ=255)
                           tCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAGTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGaa                                                                                                                        >  1:1736822/1‑150 (MQ=255)
                           tCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGACGGAAGAGCAGGTGTTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGaa                                                                                                                        >  1:329145/1‑150 (MQ=255)
                            ccGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCAATAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAg                                                                                                                       <  1:1060594/150‑1 (MQ=255)
                             cggcggAGAGTGAATGTCGAAGTCTGCTTGCTGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAgg                                                                                                                      <  1:479078/150‑1 (MQ=255)
                                   agagTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCACTAACCTGCCGGGGAGTGCAAGCTGTCCTGTTGGTGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATc                                                                                                                <  1:225080/150‑1 (MQ=255)
                                     agTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGt                                                                                                              >  1:1568403/1‑150 (MQ=255)
                                       tGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGAGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGtac                                                                                                            >  1:4635/1‑150 (MQ=255)
                                        gAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGtact                                                                                                           >  1:830965/1‑150 (MQ=255)
                                        gAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCGTCAGGGGGGGGCAGAAGGGCTAGCGtact                                                                                                           >  2:1906263/1‑150 (MQ=255)
                                         aaTGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGtacta                                                                                                          <  2:1862531/150‑1 (MQ=255)
                                           tGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACg                                                                                                        <  1:1699621/150‑1 (MQ=255)
                                           tGTCGAAGCCTGCTTCCCGGTTACATCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACg                                                                                                        >  2:1286398/1‑150 (MQ=255)
                                              cGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTc                                                                                                     <  2:1694917/150‑1 (MQ=255)
                                              cGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTc                                                                                                     <  2:937748/150‑1 (MQ=255)
                                              cGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCCATCAtactacagac                                                                                                     <  2:472031/150‑5 (MQ=255)
                                               gAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCg                                                                                                    >  1:932022/1‑150 (MQ=255)
                                                aaGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGc                                                                                                   <  2:95861/150‑1 (MQ=255)
                                                 aGCCTGCTTCCCGGTTACGTCTGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCt                                                                                                  >  1:1477132/1‑150 (MQ=255)
                                                   ccTGCTTCCCGGTTACGTCGGCGCAAAGAAGGTGCGGGAGAGAGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtc                                                                                                <  2:346699/150‑1 (MQ=255)
                                                     tGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTtctc                                                                                              <  2:1297702/150‑1 (MQ=255)
                                                             ggTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCg                                                                                      >  1:904904/1‑150 (MQ=255)
                                                                   gTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAggg                                                                                                                                  >  2:1516985/1‑100 (MQ=255)
                                                                   gTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAggg                                                                                                                                  <  1:1516985/100‑1 (MQ=255)
                                                                    tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGc                                                                                                   >  1:144821/1‑130 (MQ=255)
                                                                    tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGc                                                                                                   <  2:144821/130‑1 (MQ=255)
                                                                    tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTa                                                                               <  2:904904/150‑1 (MQ=255)
                                                                    tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTa                                                                               <  2:1568403/150‑1 (MQ=255)
                                                                     cGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAg                                                                              <  2:329145/150‑1 (MQ=255)
                                                                     cGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTGCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAg                                                                              >  1:962349/1‑150 (MQ=255)
                                                                         tCCAATAAGGTGCGGGAGAGCGGTAGTGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCCATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTGCTCGTCCACCGAAAGGTAGGAAg                                                                          <  2:948095/150‑1 (MQ=255)
                                                                         tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGtactac                                                                                                         >  1:171381/1‑119 (MQ=255)
                                                                         tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGtactac                                                                                                         <  2:171381/119‑1 (MQ=255)
                                                                         tCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAg                                                                          <  2:1617064/150‑1 (MQ=255)
                                                                           cAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCa                                                                        <  2:1868593/150‑1 (MQ=255)
                                                                            aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAgg                                                                                                                      >  2:1212437/1‑103 (MQ=255)
                                                                            aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAgg                                                                                                                      <  1:1212437/103‑1 (MQ=255)
                                                                            aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAc                                                                       <  2:699007/150‑1 (MQ=255)
                                                                            aaTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGCAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAc                                                                       >  1:1045401/1‑150 (MQ=255)
                                                                             aTAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAcc                                                                      >  1:1395213/1‑150 (MQ=255)
                                                                              tAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAccc                                                                     >  1:1408290/1‑150 (MQ=255)
                                                                                  gTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                    >  1:1741728/1‑83 (MQ=255)
                                                                                  gTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                    <  2:1741728/83‑1 (MQ=255)
                                                                                  gTGCGGGAGAGCGGTAAAGATAGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGCTGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTAGTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCAcccgcac                                                                 >  2:617413/1‑148 (MQ=255)
                                                                                         agagCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCActcc                                                          >  2:183317/1‑150 (MQ=255)
                                                                                              ggTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGGAGGAGCAAAGCCCGCCGATGGGCGGGCTCAATGTGGGGGCAGAAGCGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCAAGCACGCCtgatc                                                     >  1:1572006/1‑150 (MQ=255)
                                                                                               gTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCtcatca                                                    <  1:1352145/150‑1 (MQ=255)
                                                                                                tAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGaaa                                                                                   >  2:376141/1‑119 (MQ=255)
                                                                                                tAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAg                                                   >  2:929094/1‑150 (MQ=255)
                                                                                                gaaTGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGaaa                                                                                   <  1:376141/118‑1 (MQ=255)
                                                                                                   tgatgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAgg                                                                                                                      <  1:1596475/80‑1 (MQ=255)
                                                                                                   tgatgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAgg                                                                                                                      >  2:1596475/1‑80 (MQ=255)
                                                                                                    gatgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAAt                                               >  1:926066/1‑150 (MQ=255)
                                                                                                     gtgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                              <  2:1477132/149‑1 (MQ=255)
                                                                                                     atgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                              <  2:830965/150‑1 (MQ=255)
                                                                                                               tggtggTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAg                                    >  2:282928/1‑150 (MQ=255)
                                                                                                                  tggtAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTg                                 >  2:395091/1‑150 (MQ=255)
                                                                                                                   ggtAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAAt                                               <  2:747170/135‑1 (MQ=255)
                                                                                                                   ggtAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAAt                                               >  1:747170/1‑135 (MQ=255)
                                                                                                                     tAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGaa                                                                           >  2:653544/1‑105 (MQ=255)
                                                                                                                     tAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGaa                                                                           <  1:653544/105‑1 (MQ=255)
                                                                                                                            gAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTg                       >  1:458097/1‑150 (MQ=255)
                                                                                                                            gAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTg                       >  1:458145/1‑150 (MQ=255)
                                                                                                                                ggAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATa                   >  1:642361/1‑150 (MQ=255)
                                                                                                                                    aaTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAgg                                                                             >  1:103874/1‑88 (MQ=255)
                                                                                                                                    aaTGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAgg                                                                             <  2:103874/88‑1 (MQ=255)
                                                                                                                                       gAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTa            <  2:1395213/150‑1 (MQ=255)
                                                                                                                                         ctgaaGCCCGCCCATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATc          >  1:420138/4‑150 (MQ=255)
                                                                                                                                          tgaaGCCCGCCCATCGGCGGGCTTTAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATcc         >  1:877958/3‑150 (MQ=255)
                                                                                                                                           gaaGCCCGCCCATCGGCGGGCTTTAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCa        >  1:821377/2‑150 (MQ=255)
                                                                                                                                           aaaGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCa        <  2:1736822/150‑1 (MQ=255)
                                                                                                                                             aGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGc      <  1:183317/150‑1 (MQ=255)
                                                                                                                                             aGCCCGCCCATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                              <  1:1126216/110‑1 (MQ=255)
                                                                                                                                             aGCCCGCCCATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                              >  2:1126216/1‑110 (MQ=255)
                                                                                                                                                 cgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGt  >  1:571709/1‑150 (MQ=255)
                                                                                                                                                 cgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGt  >  2:918681/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TCCAGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGT  >  NC_002947/1274941‑1275235

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