Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,275,093 G→C 8.4% intergenic (‑284/+28) 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,0930GC 8.4% 258.2 / 12.0 83intergenic (‑284/+28)PP_5462/PP_5463hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (39/37);  new base C (4/3);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGT  >  NC_002947/1274944‑1275235
                                                                                                                                                     |                                                                                                                                               
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                                                                                                                                                <  2:645023/146‑1 (MQ=255)
     gAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATgg                                                                                                                                                >  1:645023/1‑146 (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                                                                                                                                  <  1:1516985/100‑1 (MQ=255)
                                                                gTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAggg                                                                                                                                  >  2:1516985/1‑100 (MQ=255)
                                                                 tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGc                                                                                                   >  1:144821/1‑130 (MQ=255)
                                                                 tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGc                                                                                                   <  2:144821/130‑1 (MQ=255)
                                                                 tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTa                                                                               <  2:1568403/150‑1 (MQ=255)
                                                                 tCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTa                                                                               <  2:904904/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                                                                                                                                    <  2:1741728/83‑1 (MQ=255)
                                                                               gTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAg                                                                                                                                    >  1:1741728/1‑83 (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                                                                                                                      >  2:1596475/1‑80 (MQ=255)
                                                                                                tgatgGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAgg                                                                                                                      <  1:1596475/80‑1 (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                                              >  2:1126216/1‑110 (MQ=255)
                                                                                                                                          aGCCCGCCCATCGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATc                                              <  1:1126216/110‑1 (MQ=255)
                                                                                                                                              cgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGt  >  2:918681/1‑150 (MQ=255)
                                                                                                                                              cgccgATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGt  >  1:571709/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
AGTGCGAATGGGCATTTTGGTCTATCCGGCGGAGAGTGAATGTCGAAGCCTGCTTCCCGGTTACGTCGGCTCCAATAAGGTGCGGGAGAGCGGTAATGATGGGGGTGATGGTGGTAAAGGTGAAAGGAGAATGAGCAAAGCCCGCCGATGGGCGGGCTTCAGGGGGGGGCAGAAGGGCTATCGTACTACGGTCGCTTCTCGTCCACCGAAAGGTAGGAAGCACCCGCCCGCACTCCTCATCAGGAATCGCACTTTCAGGTGACTGCCTTTGGATAAGTTGTATCCAGCCTGT  >  NC_002947/1274944‑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: