Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,174,765 A→C 5.1% intergenic (‑71/+48) glmU ← / ← atpC bifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,174,7650AC 5.1% 281.2 / 6.0 99intergenic (‑71/+48)glmU/atpCbifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base A (56/38);  new base C (2/3);  total (58/41)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

CCGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCA  >  NC_002947/6174616‑6174910
                                                                                                                                                     |                                                                                                                                                 
ccGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                   >  1:2398546/1‑150 (MQ=255)
 cGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAg                                                                                                                                                  <  2:1623743/150‑1 (MQ=255)
  gggTGCAACACGTTGGGCAGAGCTGAGGACAGGCGGGACCCATGGCCGGTAGTGAGAATAACGATTTAGAGGGATATTATTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAgg                                                                                                                                                 <  2:343814/150‑1 (MQ=255)
     tGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGCCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCtt                                                                                                                                              <  2:2102004/150‑2 (MQ=255)
       cAACACTTTGTGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAcc                                                                                                                                            <  1:2715408/150‑1 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAg                                                                                                                                                  >  1:2979749/1‑143 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAg                                                                                                                                                  <  2:2979749/143‑1 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                           >  2:1040369/1‑150 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                           >  2:2122483/1‑150 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                           >  1:2101943/1‑150 (MQ=255)
         acacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCt                                                                                                                                          >  2:1855410/1‑150 (MQ=255)
            cTTGGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtttt                                                                                                                                       <  1:1214513/150‑1 (MQ=255)
             tttGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCttttt                                                                                                                                      >  1:464851/1‑149 (MQ=255)
             tttGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCGGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      <  1:3045844/150‑1 (MQ=255)
                 ggCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtttttcttt                                                                                                                                  <  2:464851/150‑6 (MQ=255)
                   cAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCAGAGGCTACCCTTTTACTAAAt                                                                                                                                >  2:2324549/1‑150 (MQ=255)
                    aGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATc                                                                                                                               <  1:762777/150‑1 (MQ=255)
                         cGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAAc                                                                                                                          <  1:2122483/150‑1 (MQ=255)
                         cGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAAc                                                                                                                          >  1:1120386/1‑150 (MQ=255)
                             cgcATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGaa                                                                                                                      <  1:2864218/150‑1 (MQ=255)
                              gcATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAg                                                                                                                     <  1:2747805/150‑1 (MQ=255)
                               cATGCGGGTGCCTTGGCCGGCGGCGGGAATAACGATATCGAGTGACATTACTGGCTACCAGTCCTGGGCGGTCAGTGAAACAGCGGGGGAGAAAACAGAAAAAGAAAAAGGATAGCAGAGACTAACGATTTACTCAATCGCAGCAGAAgg                                                                                                                    >  1:1727357/1‑149 (MQ=255)
                                aTGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGtt                                                                                                                   >  1:1294513/1‑150 (MQ=255)
                                    gggTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCa                                                                                                               >  1:3028282/1‑150 (MQ=255)
                                          cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      >  1:1337311/1‑121 (MQ=255)
                                          cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      <  2:1337311/121‑1 (MQ=255)
                                          cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCggg                                                                                                         >  1:2634114/1‑150 (MQ=255)
                                               ccGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGa                                                                                                    <  1:2910081/150‑1 (MQ=255)
                                                cggcggCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGACCACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAt                                                                                                   <  2:2122513/150‑1 (MQ=255)
                                                cggcggCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAt                                                                                                   >  2:2822339/1‑150 (MQ=255)
                                                  gcggcgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTa                                                                                                 <  2:3028282/150‑1 (MQ=255)
                                                      cgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTg                                                                                             <  2:2398546/150‑1 (MQ=255)
                                                       gagaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCaa                                                                                                                                 <  2:2949610/113‑1 (MQ=255)
                                                       gagaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCaa                                                                                                                                 >  1:2949610/1‑113 (MQ=255)
                                                              aCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGc                                                                                                       >  1:988508/1‑132 (MQ=255)
                                                              aCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGc                                                                                                       <  2:988508/132‑1 (MQ=255)
                                                                     cGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGc                                                                              <  1:1273910/150‑1 (MQ=255)
                                                                      gactcaCATTACTGGCTATCAATCCTTGAGGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                             <  2:1259897/145‑1 (MQ=255)
                                                                      gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                             >  2:1374723/1‑150 (MQ=255)
                                                                           aCATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACg                                                                        >  1:641772/1‑150 (MQ=255)
                                                                                aCTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACg                                                                   >  2:2428380/1‑150 (MQ=255)
                                                                                    gCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGc                                                               >  2:486398/1‑150 (MQ=255)
                                                                                      tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTg                                                             >  2:2073048/1‑150 (MQ=255)
                                                                                        ccAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      <  2:1996315/75‑1 (MQ=255)
                                                                                        ccAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      >  1:1996315/1‑75 (MQ=255)
                                                                                        ccAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTggg                                                           >  1:705119/1‑150 (MQ=255)
                                                                                         cAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGc                                                          <  2:3025133/150‑1 (MQ=255)
                                                                                         cAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGc                                                          >  2:2597666/1‑150 (MQ=255)
                                                                                          aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCt                                                         >  2:1944611/1‑150 (MQ=255)
                                                                                          aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCt                                                         >  2:1269611/1‑150 (MQ=255)
                                                                                              cTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                         <  1:494707/130‑1 (MQ=255)
                                                                                              cTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCgg                                                     >  1:2948880/1‑150 (MQ=255)
                                                                                              cTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAAGCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                         >  2:494707/1‑130 (MQ=255)
                                                                                                      tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCa                                                                                                               <  2:462772/84‑1 (MQ=255)
                                                                                                      tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCa                                                                                                               >  1:462772/1‑84 (MQ=255)
                                                                                                      tCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGcc                                             >  1:642261/1‑150 (MQ=255)
                                                                                                       cAGTGAAAGGCCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTTCGGCCTCGGCCa                                            <  1:2730391/150‑1 (MQ=255)
                                                                                                       cAGTGAAAGACCGGGGGAGAAAACAGCAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCa                                            <  1:1074426/150‑1 (MQ=255)
                                                                                                            aaaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      >  2:97213/1‑55 (MQ=255)
                                                                                                            aaaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                      <  1:97213/55‑1 (MQ=255)
                                                                                                            aaaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGt                                       >  1:3057570/1‑150 (MQ=255)
                                                                                                            aaaGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGt                                       >  2:2345930/1‑150 (MQ=255)
                                                                                                               gACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcg                                    >  2:241137/1‑150 (MQ=255)
                                                                                                                aCCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcgg                                   >  1:2269866/1‑150 (MQ=255)
                                                                                                                    ggggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGACCAGACGTgcggcggc                               >  1:2169980/1‑150 (MQ=255)
                                                                                                                     gggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGc                                                       <  1:861182/125‑1 (MQ=255)
                                                                                                                     gggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcggcg                              >  2:1773526/1‑150 (MQ=255)
                                                                                                                     gggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGAATTTGCTGGACGGTACGCAGCTGGGCTGc                                                       >  2:861182/1‑125 (MQ=255)
                                                                                                                      ggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                           <  1:87034/56‑1 (MQ=255)
                                                                                                                      ggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                           >  2:87034/1‑56 (MQ=255)
                                                                                                                          aaaaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGcc                         >  1:605419/1‑150 (MQ=255)
                                                                                                                           aaaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCg                        >  2:61175/1‑150 (MQ=255)
                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                          >  2:2934093/1‑67 (MQ=255)
                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                          <  1:2934093/67‑1 (MQ=255)
                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGt                       >  1:1546513/1‑150 (MQ=255)
                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGt                       >  1:2787601/1‑150 (MQ=255)
                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGt                       >  1:2617973/1‑150 (MQ=255)
                                                                                                                              cAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGACCAGACGTGCGGCGGCGGCGCCGTAg                     <  2:2169980/150‑1 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                      <  1:3038154/116‑1 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCg                                                      >  2:3038154/1‑116 (MQ=255)
                                                                                                                                gaaaaagaaaaaGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTc                   <  1:2428380/150‑1 (MQ=255)
                                                                                                                                 aaaaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCg                  >  1:1565795/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGAc                                         <  2:651217/125‑1 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGAc                                         >  1:651217/1‑125 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa                <  2:2617973/150‑1 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa                <  2:1294513/150‑1 (MQ=255)
                                                                                                                                    aagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCAGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAg                     <  2:775859/144‑1 (MQ=255)
                                                                                                                                    aagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCAGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAg                     >  1:775859/1‑144 (MQ=255)
                                                                                                                                      gtaaaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGAGGCGGCGGCGCCGTAGTCGAAGTc             >  1:1033191/3‑150 (MQ=255)
                                                                                                                                      gaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTc             >  1:525725/1‑150 (MQ=255)
                                                                                                                                      gaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGACCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTc             >  2:3070568/1‑150 (MQ=255)
                                                                                                                                       aaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAGGTTATTAGCCGGGGCCGATTACTTTGCTTTGAGAATTTACTGGACGCTGCGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTATTCGAAGTcc            >  1:2993431/1‑150 (MQ=255)
                                                                                                                                        aaaagGGTAGCCGAGGCTACCCTTTTTCTCAATCGCAACAGACGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcg           >  2:1123087/1‑150 (MQ=255)
                                                                                                                                        aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCATCAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcg           >  1:2478915/1‑150 (MQ=255)
                                                                                                                                         aaagGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc          <  2:1062365/150‑1 (MQ=255)
                                                                                                                                         aaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc          <  1:1374723/150‑1 (MQ=255)
                                                                                                                                          aagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgcg         >  1:2101965/1‑150 (MQ=255)
                                                                                                                                          aagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgcg         >  1:2115052/1‑150 (MQ=255)
                                                                                                                                           agGGTAGCCGAGTCTACCCTTTTACTCAATCGAAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCgcgc        <  1:2324549/150‑1 (MQ=255)
                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCa  >  1:878585/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCA  >  NC_002947/6174616‑6174910

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: