Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,256,888 G→T 5.4% noncoding (433/451 nt) PP_mr14 ← cspA regulatory region

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,256,8880GT 5.4% 410.9 / 5.8 125noncoding (433/451 nt)PP_mr14cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base G (57/61);  new base T (2/5);  total (59/66)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 126 total aligned reads displayed.

CCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTT  >  NC_002947/1256743‑1257037
                                                                                                                                                 |                                                                                                                                                     
ccATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACCGTAGAAGCGCGCTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAgg                                                                                                                                                   <  2:1617585/150‑1 (MQ=38)
ccATCAGCGAAGACTGACATGCCATACCGGGTTTGCCCGCACACAGCGACTTCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAgg                                                                                                                                                   >  2:1722893/1‑150 (MQ=255)
 cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                                                                                                                  <  2:2333732/150‑1 (MQ=255)
 cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                                                                                                                  >  2:721597/1‑150 (MQ=255)
 cATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                                                                                                                  >  2:461401/1‑150 (MQ=255)
  aTCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCg                                                                                                                                                 <  2:2435178/150‑1 (MQ=255)
  aTCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCg                                                                                                                                                 >  1:1802596/1‑150 (MQ=255)
   tCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                <  1:2658700/150‑1 (MQ=255)
    cAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCggg                                                                                                                                               <  1:1424670/150‑1 (MQ=255)
      gCGAAGCCTGACCTGCCATACCCACTATGCGCGCACCCCCCGACTGCGCCGGACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtt                                                                                                                                             >  2:811838/1‑150 (MQ=37)
         aaGACTGACCTGCCATACCGAGTTTGCCCGCACACCCTGACGGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                                                                                                                          <  1:2336237/150‑1 (MQ=255)
           gACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttt                                                                                                                                        <  2:314052/150‑1 (MQ=255)
            aCTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTg                                                                                                                                       <  1:1550994/150‑1 (MQ=255)
            aCTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCATGTGGGCACTGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTg                                                                                                                                       <  1:372414/150‑1 (MQ=255)
             cTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGTCCTACTTGAAGCCTGTGCAGTCCCTGTCGGCAGGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                      <  1:687213/150‑3 (MQ=255)
               gACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttatgc                                                                                                                                    <  1:2953440/150‑5 (MQ=38)
               gACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                    <  1:461401/150‑1 (MQ=255)
                aCCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAg                                                                                                                                   >  1:1045315/1‑150 (MQ=255)
                   tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGGCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGGTTTTTTGAAAGCTa                                                                                                                                >  1:878116/1‑150 (MQ=255)
                   tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCg                                                                                                                                                 <  1:266046/133‑1 (MQ=255)
                   tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCg                                                                                                                                                 >  2:266046/1‑133 (MQ=255)
                   tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGtttt                                                                                                                                           >  2:411928/1‑139 (MQ=255)
                   tGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTa                                                                                                                                <  2:2971660/150‑1 (MQ=255)
                    gCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTAATTGAAGACTGAACAATCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAg                                                                                                                               >  1:1302951/1‑150 (MQ=255)
                     ccATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                >  2:1076703/1‑132 (MQ=255)
                       aTACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGaa                                                                                                                            >  1:1514371/1‑150 (MQ=255)
                        tACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAg                                                                                                                           <  1:1418651/150‑1 (MQ=255)
                          ccGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAATCCGCAAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCt                                                                                                                         >  1:1106802/1‑150 (MQ=255)
                          ccGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCt                                                                                                                         >  1:47194/1‑150 (MQ=255)
                           cGAGTTTGCCGGCACACAGCGACTGCGGCGTACTTGAAGCCTGGGCAGTCCCTGTGGGAGCTGGCGAGCCCTCGAAGAGGACCGCGCCTATTTCGAGGAGGCATAAAAAAACCGGCCAGAAGGCGGGTTTGTTTGAAAGCTAGGAAGCTa                                                                                                                        <  2:1056132/150‑1 (MQ=255)
                           cGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTCCGAGCAGGCATAAAAAAACCCGCCAGAAGGc                                                                                                                                                  <  1:661788/124‑1 (MQ=255)
                                    ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                >  2:2752090/1‑117 (MQ=255)
                                    ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCgg                                                                                                                                                <  1:2752090/117‑1 (MQ=255)
                                    ccGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACtt                                                                                                               <  2:2365076/150‑1 (MQ=255)
                                      gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAg                                                                                                             <  2:2705472/150‑1 (MQ=255)
                                      gCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAg                                                                                                             <  2:1302951/150‑1 (MQ=255)
                                        acacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttta                                                                                                                                       >  2:449425/1‑121 (MQ=255)
                                        acacaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttttta                                                                                                                                       <  1:449425/122‑2 (MQ=255)
                                          acaGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGa                                                                                                         >  1:1037825/1‑150 (MQ=255)
                                               gACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGa                                                                                                    >  1:2026808/1‑150 (MQ=255)
                                                 cTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                  <  1:176517/150‑1 (MQ=255)
                                                       cGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACttcttc                                                                                            >  1:655248/1‑150 (MQ=255)
                                                          aCTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGc                                                                                         <  2:2026808/150‑1 (MQ=255)
                                                              gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGttttt                                                                                                                                          >  1:1004821/1‑97 (MQ=255)
                                                              gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGa                                                                                                                             <  1:2408825/110‑1 (MQ=255)
                                                              gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGc                                                                                     <  1:2177532/150‑1 (MQ=255)
                                                              gAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGc                                                                                     >  1:1967784/1‑150 (MQ=255)
                                                                 gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATc                                                                                                                    >  2:2365681/1‑116 (MQ=255)
                                                                 gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATc                                                                                                                    <  1:2365681/116‑1 (MQ=255)
                                                                 gCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATg                                                                                  >  1:1641305/1‑150 (MQ=255)
                                                                    tgtgCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCt                                                                               <  2:2154476/150‑1 (MQ=255)
                                                                        cAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                    >  2:2136286/1‑93 (MQ=255)
                                                                        cAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                    <  1:2136286/93‑1 (MQ=255)
                                                                          gTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTgg                                                                         >  2:2579571/1‑150 (MQ=255)
                                                                             ccTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGATAGCAATCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcg                                                                      >  2:2255504/1‑150 (MQ=255)
                                                                              cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                      <  2:1669554/85‑3 (MQ=25)
                                                                              cTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgc                                                                     >  2:2584327/1‑150 (MQ=255)
                                                                               tgtgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATTGACGCAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCcgcg                                                                    >  1:1870161/1‑150 (MQ=255)
                                                                                 tgGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGTCCGCggg                                                                  <  1:721597/150‑1 (MQ=255)
                                                                                    gAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATc                                                                                                       <  2:1125331/110‑1 (MQ=255)
                                                                                    gAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATc                                                                                                       >  1:1125331/1‑110 (MQ=255)
                                                                                       cGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcga                                                            <  2:1037825/150‑1 (MQ=255)
                                                                                       cGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcga                                                            <  2:2120779/150‑1 (MQ=255)
                                                                                         ggCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaa                                                                                                                                     >  2:1492460/1‑75 (MQ=255)
                                                                                         ggCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGcgacg                                                          >  1:1117298/1‑150 (MQ=255)
                                                                                              gCCCTCAAAGAGGCCAGCGCCTATTTCGAGCAGGCATAACAAACCACGCCGTGGCTCGTGTTTTTTTGAAAGCTATGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGt                                                     <  1:689235/150‑1 (MQ=255)
                                                                                              gCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGt                                                     <  2:878116/150‑1 (MQ=255)
                                                                                                ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                      <  2:2596505/67‑3 (MQ=18)
                                                                                                ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCTTCTGGCGGGtttttttat                                                                                                                                      >  1:2596505/1‑65 (MQ=18)
                                                                                                ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGaa                                                                                                                            <  2:1707504/77‑1 (MQ=255)
                                                                                                ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGaa                                                                                                                            >  1:1707504/1‑77 (MQ=255)
                                                                                                ccGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTaa                                                   <  2:1080886/150‑1 (MQ=255)
                                                                                                  gcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACa                                                 <  2:1277113/150‑1 (MQ=255)
                                                                                                  gcgAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACa                                                 <  1:2803956/150‑1 (MQ=255)
                                                                                                         ggCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGc                                          <  1:37499/150‑1 (MQ=255)
                                                                                                         ggCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGc                                          <  1:2601265/150‑1 (MQ=255)
                                                                                                          gCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                    <  2:1220912/59‑1 (MQ=255)
                                                                                                          gCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGaaa                                                                                                                                    >  1:1220912/1‑59 (MQ=255)
                                                                                                            cAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCctt                                       <  2:1641305/150‑1 (MQ=255)
                                                                                                                gcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCttt                                   >  1:359408/1‑150 (MQ=255)
                                                                                                                gcCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTAAGCAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCttt                                   <  1:1722893/150‑1 (MQ=255)
                                                                                                                  cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGc                                               >  2:1988114/1‑136 (MQ=255)
                                                                                                                  cTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGc                                               <  1:1988114/136‑1 (MQ=255)
                                                                                                                   tATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAg                                <  1:2085617/150‑1 (MQ=255)
                                                                                                                         gAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGtttt                          >  1:102198/1‑150 (MQ=255)
                                                                                                                            cAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGaa                       <  2:1967784/150‑1 (MQ=255)
                                                                                                                               gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTTTGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGcc                    <  1:811838/150‑1 (MQ=255)
                                                                                                                               gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                  >  2:161289/1‑72 (MQ=255)
                                                                                                                               gCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAAc                                                                                                  <  1:161289/72‑1 (MQ=255)
                                                                                                                                    aaaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAg               >  2:1732784/1‑150 (MQ=255)
                                                                                                                                     aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTa                                                                <  2:2855434/100‑1 (MQ=255)
                                                                                                                                     aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAg                                                               >  1:2724531/1‑101 (MQ=255)
                                                                                                                                     aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAg                                                               <  2:2724531/101‑1 (MQ=255)
                                                                                                                                     aaaaaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGc              >  1:1388576/1‑150 (MQ=255)
                                                                                                                                      aaaaCCCGCCAGAAGGCGGGTTTTTTTAAAAACTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCt             >  1:80401/1‑150 (MQ=255)
                                                                                                                                        aaCCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTg           <  1:224834/150‑1 (MQ=255)
                                                                                                                                            cGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               >  1:2806704/1‑62 (MQ=255)
                                                                                                                                            cGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAActtc                                                                                               <  2:2806704/62‑1 (MQ=255)
                                                                                                                                             gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCtt            >  2:1844656/1‑144 (MQ=255)
                                                                                                                                             gCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATgg      <  2:522764/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
CCATCAGCGAAGACTGACCTGCCATACCGAGTTTGCCCGCACACAGCGACTGCGCCGTACTTGAAGCCTGTGCAGTCCCTGTGGGAGCGGGCGAGCCCGCGAAGAGGCCAGCGCCTATTTCGAGCAGGCATAAAAAAACCCGCCAGAAGGCGGGTTTTTTTGAAAGCTAGGAAGCTAGATCGACTTAGGCGATCTGAACTTCTTCAGCCTGCATGCCTTTCTGGCCGCGGGTAGCGACGAAAGTAACAGCCTGGCCTTCTTTCAGGGTTTTGAAGCCGTCAGCTTGGATGGCTTT  >  NC_002947/1256743‑1257037

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