Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,564,769 A→G 5.5% intergenic (‑865/‑33) vgrG‑III ← / → glgA type VI secretion system protein/glycogen synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,564,7690AG 5.5% 229.3 / 2.1 73intergenic (‑865/‑33)vgrG‑III/glgAtype VI secretion system protein/glycogen synthase
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base G (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.53e-01

CCAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCAC  >  NC_002947/4564620‑4564915
                                                                                                                                                     |                                                                                                                                                  
ccAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                    <  1:131596/150‑1 (MQ=255)
ccAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGACGGAGCGGCGATCGCGCCGTGCGTGCGCGCACCGTACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                    <  1:925967/150‑1 (MQ=255)
    ccACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCAc                                                                                                                                                >  2:753715/1‑150 (MQ=255)
      aCGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCAccc                                                                                                                                              >  1:576396/1‑150 (MQ=255)
        ggCCGCCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGa                                                                                                                                            >  2:508971/1‑150 (MQ=255)
         gCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCgtg                                                                                                                                           <  1:753715/150‑3 (MQ=255)
                  ccaacAGCCCAATCGCTTTCGGCCGGAACGCCGACCGCGCTCTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGCAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                  <  1:1030833/150‑1 (MQ=255)
                    aacAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAggg                                                                                                                                >  2:240444/1‑150 (MQ=255)
                         cccAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGga                                                                                                                           >  1:701537/1‑150 (MQ=255)
                         cccAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGga                                                                                                                           >  2:689017/1‑150 (MQ=255)
                         cccAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGga                                                                                                                           <  2:52240/150‑1 (MQ=255)
                            aaTCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtgg                                                                                                                                          >  1:1282759/1‑129 (MQ=255)
                            aaTCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtgg                                                                                                                                          <  2:1282759/132‑4 (MQ=255)
                            aaTCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                        <  1:1076771/134‑1 (MQ=255)
                            aaTCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCaa                                                                                                                                        >  2:1076771/1‑134 (MQ=255)
                                  tttCGACCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTAt                                                                                                                  >  1:112857/1‑150 (MQ=255)
                                       gCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                    >  2:535211/1‑111 (MQ=255)
                                       gCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCa                                                                                                                                                    <  1:535211/111‑1 (MQ=255)
                                       gCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCa                                                                                                             >  2:1072485/1‑150 (MQ=255)
                                        ccGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAg                                                                                                            >  2:519853/1‑150 (MQ=255)
                                         cGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGt                                                                                                           <  2:112857/150‑1 (MQ=255)
                                          ggAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTg                                                                                                          <  1:689017/150‑1 (MQ=255)
                                                  gATCGCGCCGTGCGTGCGCGCACGAGGGAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTc                                                                                                  >  1:814406/1‑150 (MQ=255)
                                                  gATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTc                                                                                                  >  1:392143/1‑150 (MQ=255)
                                                       cgcCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCAc                                                                                                                                                <  2:324469/99‑1 (MQ=255)
                                                       cgcCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCAc                                                                                                                                                >  1:324469/1‑99 (MQ=255)
                                                         ccGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGTCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTc                                                                                           <  2:169477/150‑1 (MQ=255)
                                                                     gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                  <  1:1119777/99‑1 (MQ=255)
                                                                     gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAg                                                                                                                                  >  2:1119777/1‑99 (MQ=255)
                                                                     gcACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCtt                                                                               <  1:985979/150‑1 (MQ=255)
                                                                      cACGACACAGGGGCATGGAGCGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                              <  2:842709/150‑1 (MQ=255)
                                                                      cACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                              <  1:519853/150‑1 (MQ=255)
                                                                      cACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCAGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                              <  1:577439/150‑1 (MQ=255)
                                                                           cacaGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCACGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGCCGAGCCTCACGTAGACGCTTTCAAAc                                                                         <  1:814388/150‑1 (MQ=255)
                                                                             caGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCg                                                                       <  1:1072485/150‑1 (MQ=255)
                                                                             caGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCg                                                                       <  1:249690/150‑1 (MQ=255)
                                                                              aGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                              <  2:536258/126‑1 (MQ=255)
                                                                              aGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGc                                                                                              >  1:536258/1‑126 (MQ=255)
                                                                              aGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCAc                                                                                         >  2:431203/1‑131 (MQ=255)
                                                                              aGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCAc                                                                                         <  1:431203/131‑1 (MQ=255)
                                                                                  gCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtg                                                                                                                                           >  2:411013/1‑75 (MQ=255)
                                                                                  gCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCTGGCACCCgtg                                                                                                                                           <  1:411013/77‑3 (MQ=255)
                                                                                  gCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATaa                                                                  >  2:410219/1‑150 (MQ=255)
                                                                                  gCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATaa                                                                  <  2:576396/150‑1 (MQ=255)
                                                                                   cATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAAc                                                                 <  1:508971/150‑1 (MQ=255)
                                                                                    aTGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAAcc                                                                >  1:357364/1‑150 (MQ=255)
                                                                                     tGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcg                                                               >  2:1163486/1‑150 (MQ=255)
                                                                                      ggAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcgc                                                              <  2:1069527/150‑1 (MQ=255)
                                                                                      ggAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACcgc                                                              >  2:320807/1‑150 (MQ=255)
                                                                                             aCTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                              >  1:682757/1‑127 (MQ=255)
                                                                                             aCTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCttt                                                                              <  2:682757/127‑1 (MQ=255)
                                                                                                  gcagcaTTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCAcc                                                  >  1:1170570/1‑150 (MQ=255)
                                                                                                      caTTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCgg                                              >  2:1128145/1‑150 (MQ=255)
                                                                                                         tCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGAtt                                           >  2:1269631/1‑150 (MQ=255)
                                                                                                           agcaagcaGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTc                                         <  2:392143/150‑1 (MQ=255)
                                                                                                              aagcaGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGc                                                          <  2:1267687/130‑1 (MQ=255)
                                                                                                              aagcaGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGc                                                          >  1:1267687/1‑130 (MQ=255)
                                                                                                               gcactggagaACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAg                                                           >  2:688800/7‑128 (MQ=255)
                                                                                                               gcactggagaACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAg                                                           <  1:688800/122‑1 (MQ=255)
                                                                                                                    agagaACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACa                                <  2:357364/150‑1 (MQ=255)
                                                                                                                     gagaACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACag                               <  2:1170570/150‑1 (MQ=255)
                                                                                                                        aaCCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCACTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACaggca                            >  2:1068593/1‑150 (MQ=255)
                                                                                                                         aCCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACaggcaa                           <  1:240444/150‑1 (MQ=255)
                                                                                                                          ccAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCAcccc                                                >  2:384128/1‑128 (MQ=255)
                                                                                                                          ccAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCAcccc                                                <  1:384128/128‑1 (MQ=255)
                                                                                                                           cAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACaggcaagg                         >  1:83682/1‑150 (MQ=255)
                                                                                                                                 gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTc                                                                                                  <  1:1074683/71‑1 (MQ=255)
                                                                                                                                 gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTc                                                                                                  >  2:1074683/1‑71 (MQ=255)
                                                                                                                                 gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCAcc                                                  >  1:325539/1‑119 (MQ=255)
                                                                                                                                 gAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCAcc                                                  <  2:325539/119‑1 (MQ=255)
                                                                                                                                    aCATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAAc                                                                         >  2:416420/1‑93 (MQ=255)
                                                                                                                                    aCATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAAc                                                                         <  1:416420/93‑1 (MQ=255)
                                                                                                                                          ccACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAgc          >  1:768027/1‑150 (MQ=255)
                                                                                                                                          ccACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGACACGCGACGCCCAgc          >  2:357735/1‑150 (MQ=255)
                                                                                                                                            tggggTGCCAGCAACCGACGAAGCTGAGGGCCGTAGAGAAGTATGAACAGTGCAGCTGTCGAGCGGCACGGAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACACCGGATGTAGCCACGACAGGCACGGCACCAGTCGCACAgcgc        >  2:212165/3‑150 (MQ=255)
                                                                                                                                                gTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGc    <  2:701537/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGcac  >  2:385902/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CCAGCCACGGCCACCAACCCAACAGCCCAATCGCTTTCGGCCGGAACGCCGATCGCGCCGTGCGTGCGCGCACGACACAGGGGCATGGAACGCACTACGCAGCATTCAGCAAGCAGAGAGAACCAGCACGAAACATCGCCACGGGTGCCAGCACCCGACCAAGCATAGGGACGGAGAGAAGTATGATCAGTGCCGCTGTCGAGCCTCACGTAGACGCTTTCAAACCGGATAACCGCGAGCCGCTCACCCCGGATTTCGCCACGACAGGCAAGGCACCCGGCGCCCAGCGCCAGCAC  >  NC_002947/4564620‑4564915

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