Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,174,764 G→T 7.4% intergenic (‑70/+49) 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,7640GT 7.4% 361.2 / 11.3 116intergenic (‑70/+49)glmU/atpCbifunctional N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase/ATP synthase subunit epsilon
Reads supporting (aligned to +/- strand):  ref base G (58/49);  new base T (3/6);  total (61/55)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

Only 100 of 118 total aligned reads displayed.

CCGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGA  >  NC_002947/6174616‑6174912
                                                                                                                                                    |                                                                                                                                                    
ccGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                     >  1:1755212/1‑150 (MQ=255)
ccGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                     <  2:503654/150‑1 (MQ=255)
ccGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGa                                                                                                                                                     <  2:560209/150‑1 (MQ=255)
     tGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACGGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTa                                                                                                                                                >  1:1219511/1‑150 (MQ=255)
     tGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTa                                                                                                                                                >  1:735896/1‑150 (MQ=255)
       cAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAAAAACGATATGAAGTGACATTACTGGGTACGAATCCTGGGCGGTTAGGTAAAGCTCGTGTGAGAGAACAGAAAAAGAAAAAGGGTAGCCGAGGCTagc                                                                                                                                              >  1:251850/1‑148 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                             >  2:1552247/1‑150 (MQ=255)
        aaCACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCAAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAAAAAAAGGGTAGCCGAGGCTAccc                                                                                                                                             >  1:640102/1‑150 (MQ=255)
         acacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCt                                                                                                                                            >  2:417698/1‑150 (MQ=255)
          cacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAAAAAAAGAAAAAGGGTAGCCTCGGCTACCCtt                                                                                                                                           <  2:866562/150‑1 (MQ=255)
          cacTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGCCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtt                                                                                                                                           <  1:505839/150‑1 (MQ=255)
           acTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCttt                                                                                                                                          >  2:2227168/1‑150 (MQ=255)
            cTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtttt                                                                                                                                         <  2:2366329/150‑1 (MQ=255)
            cTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCTCGGCTACCCtttt                                                                                                                                         <  2:2093274/150‑1 (MQ=255)
             tttGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTa                                                                                                                                        >  1:563603/1‑150 (MQ=255)
               tGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACt                                                                                                                                      <  1:1571229/150‑1 (MQ=255)
                                aTGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTAc                                                                                                                                       >  1:1486703/1‑132 (MQ=255)
                                aTGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTAc                                                                                                                                       <  2:1486703/132‑1 (MQ=255)
                                    gggTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                             <  2:1237771/138‑1 (MQ=255)
                                    gggTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCaa                                                                                                                             >  1:1237771/1‑138 (MQ=255)
                                     ggTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAg                                                                                                                <  2:1796178/150‑1 (MQ=255)
                                      gTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGc                                                                                                               >  2:1274541/1‑150 (MQ=255)
                                         ccTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                            <  1:209518/150‑1 (MQ=255)
                                          cTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCggg                                                                                                           <  2:1996712/150‑1 (MQ=255)
                                           ttGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAAcagcag                                                                                                                       <  1:2403406/137‑1 (MQ=255)
                                           ttGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAg                                                                                                                       >  2:2403406/1‑137 (MQ=255)
                                                     gcgAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACtt                                                                                                <  1:359427/150‑1 (MQ=255)
                                                         gaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGaa                                                                                                                        <  1:576574/122‑1 (MQ=255)
                                                         gaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGaa                                                                                                                        >  2:576574/1‑122 (MQ=255)
                                                         gaATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACtt                                                                                                >  1:1883028/1‑146 (MQ=255)
                                                           aTAACGATACCGAGTGACATTACTGGCGACCAATCCTGGGCGGTCAGTGAAAGACCGGGGCAGAAAACAGAAAAAGACAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCttt                                                                                          <  1:1372921/150‑1 (MQ=255)
                                                             aaCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGCAGTTATCAGCCGGGGCCGAGGACGTGCCTGTgc                                                                                        >  2:1720561/1‑150 (MQ=255)
                                                             aaCGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                        <  1:417698/150‑1 (MQ=255)
                                                               cGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTAt                                                                                                                   <  1:1300369/121‑1 (MQ=255)
                                                               cGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTAt                                                                                                                   >  2:1300369/1‑121 (MQ=255)
                                                                 atatCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCAt                                                                                    >  2:2406273/1‑150 (MQ=255)
                                                                      gAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCt                                                                               <  1:496293/150‑1 (MQ=255)
                                                                        gTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTgg                                                                             <  2:1219511/150‑1 (MQ=255)
                                                                          gACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGAc                                                                           >  1:972355/1‑150 (MQ=255)
                                                                           aCATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACg                                                                          <  2:11335/150‑1 (MQ=255)
                                                                           aCATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACg                                                                          >  1:1868311/1‑150 (MQ=255)
                                                                               tACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTAc                                                                      <  1:574327/150‑1 (MQ=255)
                                                                                  tGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCa                                                                   >  2:1493839/1‑150 (MQ=255)
                                                                                    gCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGc                                                                 >  1:814253/1‑150 (MQ=255)
                                                                                     cTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCt                                                                <  1:1274541/150‑1 (MQ=255)
                                                                                      tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTg                                                               >  2:104652/1‑150 (MQ=255)
                                                                                      tACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTg                                                               >  2:110344/1‑150 (MQ=255)
                                                                                       aCCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgcgc                                                                                      <  2:320508/126‑1 (MQ=255)
                                                                                       aCCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgcgc                                                                                      >  1:320508/1‑126 (MQ=255)
                                                                                        ccAATCCTGGGCGGTAAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTggg                                                             <  1:873075/150‑1 (MQ=255)
                                                                                         cAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATAAACCGGTGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGc                                                            >  2:677322/1‑150 (MQ=255)
                                                                                          aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                            <  2:314610/101‑1 (MQ=255)
                                                                                          aaTCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCgg                                                                                                            >  1:314610/1‑101 (MQ=255)
                                                                                              cTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGCCGGTACGCAGCTGGGCTGCgg                                                       >  2:849376/1‑150 (MQ=255)
                                                                                                    ggTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCgg                                                 >  1:141958/1‑150 (MQ=255)
                                                                                                     gTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGc                                                >  1:2278217/1‑150 (MQ=255)
                                                                                                          tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                        <  1:1163996/105‑1 (MQ=255)
                                                                                                          tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                        >  2:1163996/1‑105 (MQ=255)
                                                                                                          tGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTATGTAGCTGGGCTGCGGCCTCGGCCAGAc                                           <  1:677322/150‑1 (MQ=255)
                                                                                                            aaaGACCGGGGGTGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGCTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGt                                         <  1:1098115/150‑1 (MQ=255)
                                                                                                                 ccGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTTCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggc                                    <  2:2013600/150‑1 (MQ=255)
                                                                                                                   gggggAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTgcggcgg                                  <  1:2227168/150‑1 (MQ=255)
                                                                                                                            aaCAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGCGCTGCGGCCTCGGCCATGCGTGCGGCGGCGGCGCCGt                         >  2:21694/1‑150 (MQ=255)
                                                                                                                             aCAGAAAAATAAAAAGGGTAGTCTTTGCTACCGGTTTACTCAATCGCAACAGAAGTTATCAGCGGGGGCCTATTACTTTCCTTTACGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTTCGGCGGCGGCGCCGTa                        <  2:640102/150‑1 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCGGAGGCTACCCTTTTAATCAATTGCCACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTACGGCGGCGGCGCCGTAGt                      >  1:864068/1‑150 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGc                                                                                                               <  1:1825685/61‑1 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGc                                                                                                               >  2:1825685/1‑61 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTgc                                                                                        >  2:1487758/1‑84 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTgg                                                              >  1:1555797/1‑110 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTgg                                                              <  2:1555797/110‑1 (MQ=255)
                                                                                                                               agaaaaagaaaaAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGt                      >  2:1895836/1‑150 (MQ=255)
                                                                                                                                gaaaaagaaaaaGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAACTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTc                     >  2:2006018/1‑150 (MQ=255)
                                                                                                                                 aaaaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCg                    >  1:923481/1‑150 (MQ=255)
                                                                                                                                  aaaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGGCGa                   >  1:411153/1‑150 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAt                                                                                                     <  2:2359303/67‑1 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGAt                                                                                                     >  1:2359303/1‑67 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa                  <  2:1868311/150‑1 (MQ=255)
                                                                                                                                   aaagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGaa                  >  2:1373868/1‑150 (MQ=255)
                                                                                                                                    aagaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAg                                             <  1:161701/122‑1 (MQ=255)
                                                                                                                                     agaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGt                <  2:534696/150‑1 (MQ=255)
                                                                                                                                      gaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGc                                                         >  1:165872/1‑108 (MQ=255)
                                                                                                                                      gaaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGc                                                         <  2:165872/108‑1 (MQ=255)
                                                                                                                                       aaaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCACTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTcc              <  1:1373868/150‑1 (MQ=255)
                                                                                                                                        aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTa                                                                                                   <  2:1363219/64‑1 (MQ=255)
                                                                                                                                        aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTa                                                                                                   >  1:1363219/1‑64 (MQ=255)
                                                                                                                                        aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGt                <  2:2233315/147‑1 (MQ=255)
                                                                                                                                        aaaagGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGt                >  1:2233315/1‑147 (MQ=255)
                                                                                                                                           agGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACt                                                                                                 >  2:2368848/1‑63 (MQ=38)
                                                                                                                                           agGGTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACt                                                                                                 <  1:2368848/63‑1 (MQ=38)
                                                                                                                                             ggTAGCCTCGGCTACCCTTTTTCTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCt        >  1:1535395/1‑150 (MQ=255)
                                                                                                                                             ggTAGCCTCGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCt        >  2:119205/1‑150 (MQ=255)
                                                                                                                                               tAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCttt      <  2:563603/150‑1 (MQ=255)
                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCa    <  2:1755212/150‑1 (MQ=255)
                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCa    <  1:318861/150‑1 (MQ=255)
                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCa    <  1:1475840/150‑1 (MQ=255)
                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCa    >  2:318861/1‑150 (MQ=255)
                                                                                                                                                 gCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCa    <  2:141958/150‑1 (MQ=255)
                                                                                                                                                   cGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGa  >  1:1865594/1‑150 (MQ=255)
                                                                                                                                                    gAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc            <  2:33755/139‑1 (MQ=255)
                                                                                                                                                    gAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCcgc            >  1:33755/1‑139 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCGGGTGCAACACTTTGGGCAGCGCCGAGCGCATGCGGGTGCCTTGGCCGGCGGCGAGAATAACGATATCGAGTGACATTACTGGCTACCAATCCTGGGCGGTCAGTGAAAGACCGGGGGAGAAAACAGAAAAAGAAAAAGGGTAGCCGAGGCTACCCTTTTACTCAATCGCAACAGAAGTTATCAGCCGGGGCCGATTACTTGCCTTTGCGCATTTGCTGGACGGTACGCAGCTGGGCTGCGGCCTCGGCCAGACGTGCGGCGGCGGCGCCGTAGTCGAAGTCCGCGCCTTTCAGA  >  NC_002947/6174616‑6174912

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: