Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,700,770 G→A 6.6% intergenic (+82/‑85) hslU → / → PP_5002 protease HslVU ATPase subunit/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,700,7700GA 6.6% 276.4 / 6.9 91intergenic (+82/‑85)hslU/PP_5002protease HslVU ATPase subunit/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (41/44);  new base A (2/4);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGCCTCA  >  NC_002947/5700621‑5700912
                                                                                                                                                      |                                                                                                                                              
cggcgTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCgggg                                                                                                                                                <  1:1460785/150‑1 (MQ=255)
cggcgTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGGTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCgggg                                                                                                                                                <  1:2032371/150‑1 (MQ=255)
  gCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgc                                                                                                                                              >  1:1198048‑M2/1‑148 (MQ=255)
  gCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgc                                                                                                                                              >  1:940966‑M2/1‑148 (MQ=255)
  gCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgc                                                                                                                                              >  1:3766111‑M2/1‑148 (MQ=255)
     tATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcggg                                                                                                                                           <  1:1311098‑M2/150‑6 (MQ=255)
     tATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcggg                                                                                                                                           >  1:4938381‑M2/1‑145 (MQ=255)
     tATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcggg                                                                                                                                           >  1:2203543‑M2/1‑145 (MQ=255)
      aTGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggc                                                                                                                                          >  2:3373072‑M2/1‑144 (MQ=255)
      aTGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggc                                                                                                                                          >  2:3373113‑M2/1‑144 (MQ=255)
         tgAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttg                                                                                                                                       >  2:2246875‑M2/1‑141 (MQ=255)
          gAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgc                                                                                                                                      >  1:461153‑M2/1‑140 (MQ=255)
          gAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgc                                                                                                                                      >  2:418626‑M2/1‑140 (MQ=255)
           aaCAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcc                                                                                                                                     >  1:734502‑M2/1‑139 (MQ=255)
            aCAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTGGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgccc                                                                                                                                    >  1:3354560‑M2/1‑138 (MQ=255)
                ccACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcaagcccgccc                                                                                                                                >  1:1548838‑M2/1‑134 (MQ=255)
                  aCCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTGGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAag                                                                                                                              <  1:3113938/150‑1 (MQ=255)
                   ccTGGGTGAGCTGGCACAGAACGAAGACCTGTCGGGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTTTGGTGATCCGTGGGGGGGCGGGCTTGCCCGCGAaga                                                                                                                             <  2:3741497/150‑1 (MQ=255)
                                gaacaTAACGAACACCTGTCGCTTTACTTTCTTTAAGACCGGTTTCGACCGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGGGGGCTTGCCCGCGAAGAGGCCCTTGAAtcc                                                                                                                <  1:815300/148‑1 (MQ=255)
                                  acaGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGACAAGCCCCCTCCAGCACGCGGCACCACAACCCTTGGGCCTGTGGTGTTACGTGCGGGGGCGGGCGTGCCCGCGAAGAGGCCCTTGAAtacct                                                                                                              >  1:1292165/1‑150 (MQ=255)
                                     gAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAtccctccc                                                                                                           <  2:398350/150‑1 (MQ=255)
                                              gcTGTCGGTTTACATTGTTTATGACGGTTTTCGAATGTTTGGCGGGAATGCCGGCTCCTGCAGGGGCCACCACAACCGGTGGGCCTGTGGTGTTGCGTGGGGGGGGGGGGTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCt                                                                                                  <  1:4846080/149‑1 (MQ=255)
                                                      gTTACTTTCTTTAAGACGGGTTTGGACTGTTTCGCGGGCATGCCCGCTCCTGCACGCGCCACCACAACCATTGGGCGTGTGGTTTTCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCtg                                                                                          <  2:231172/150‑1 (MQ=255)
                                                          cATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCt                                                                                      <  2:3335446/150‑1 (MQ=255)
                                                              cTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTccat                                                                                  <  2:1721480/150‑1 (MQ=255)
                                                               ctgtactccGGTTTCCGACAGTTGGGGGGGCACCCCGGCCCGGGCCCGCCCCACCAAACCCCGGGGGCCGGGGGTGACCGGGGGGGGGGCGGGGTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcat                                                                                  <  1:2887004/143‑1 (MQ=255)
                                                                gTATGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCGTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                                                                                <  2:3434903/150‑1 (MQ=255)
                                                                gTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCCCACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                                                                                <  2:4710398/150‑1 (MQ=255)
                                                                gTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTtcatca                                                                                <  1:4871708/150‑1 (MQ=255)
                                                                    gACCGGTTTCGACTGTTGCGCGGGCAAGCCCGCTCCTGCACGCGCCACCCCAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTc                                                                            <  1:4060187/150‑1 (MQ=255)
                                                                    gACCGGTTTCGACAGTTTCGCGGGCAAGCCCCCGCCAGCACCCGCCACCACAAACCAAGGGGCGGTGGTGTTCCCGGCGGGGGCGGGCGTGCCCCCGCAAAGGCCCATTTATCCCACCCCCAAAAACAGGAAACAGTGCGtcaaaaccac                                                                            >  2:1302565/1‑143 (MQ=255)
                                                                      ccGGTTTCGACAGGTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACCACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTccc                                                                          <  2:5063572/150‑1 (MQ=255)
                                                                        ggTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCga                                                                        <  2:433589/150‑1 (MQ=255)
                                                                              cGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcg                                                                                                                                             <  1:4477155‑M2/76‑4 (MQ=255)
                                                                              cGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcg                                                                                                                                             >  2:4477155‑M2/1‑73 (MQ=255)
                                                                               gccAGTTTCGCGGGCACGCCCGCACCGCCACGCGCCACCACACCCCTTGGGCCTGTGGTTAACCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGAtt                                                                  <  1:3603861/148‑1 (MQ=255)
                                                                                  aGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCCCCACCACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATCAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTg                                                               <  2:4404405/150‑1 (MQ=255)
                                                                                    tttCGCGGGCAAGCCCGCACCTTCACGCGCCACCAAAAACCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCCAAGAGGCCCTTGAATCCCTCCCCCATAAGCAGGAAGCTGTGCTTCAACAGCTCCCGAGAGAATCTggc                                                             >  2:2924422/1‑148 (MQ=255)
                                                                                      tCGCGGGCAAGCCCGCTCCTGCCCGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCa                                                           <  2:1428562/150‑1 (MQ=255)
                                                                                              acgCCCGCTCCTCCCCGCGCCACCACACCCCGTGGGCGTGTGGTTAGCGGTGGGGGGGCGGGGTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgc                                                   <  1:88075/148‑1 (MQ=255)
                                                                                              aaGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcac                                                   <  2:2011563‑M2/150‑94 (MQ=255)
                                                                                              aaGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctg                                                        <  2:851049‑M2/145‑89 (MQ=255)
                                                                                              aaGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctg                                                        >  1:851049‑M2/1‑57 (MQ=255)
                                                                                                gCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcct                                                 <  2:4046068/150‑1 (MQ=255)
                                                                                                gCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcct                                                 >  1:4256572/1‑150 (MQ=255)
                                                                                                gCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCAAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcct                                                 >  1:3502301/1‑150 (MQ=255)
                                                                                                gCCCGCTCCTGCACGCGCCACCACAACCCTTGCGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcct                                                 <  2:4046019/150‑1 (MQ=255)
                                                                                                 cccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctg                                                <  2:3940708/150‑1 (MQ=255)
                                                                                                 cccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctg                                                >  1:1000150/1‑150 (MQ=255)
                                                                                                  ccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatac                                               <  2:3617077‑M2/150‑98 (MQ=255)
                                                                                                  ccGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctgc                                               >  2:3309343/1‑150 (MQ=255)
                                                                                                   cGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacg                                              <  2:1607712‑M2/150‑99 (MQ=255)
                                                                                                    gCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgc                                             <  1:3487822‑M2/150‑100 (MQ=255)
                                                                                                    gCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCgcctgcct                                             >  2:4962075/1‑150 (MQ=255)
                                                                                                      tCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccg                                           <  2:1095031‑M2/150‑102 (MQ=255)
                                                                                                        cTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCg                                         >  1:2263390/1‑150 (MQ=255)
                                                                                                         tGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGc                                        >  2:1349650/1‑150 (MQ=255)
                                                                                                         tGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGc                                        >  2:3674108/1‑150 (MQ=255)
                                                                                                          gCACGCGCCACCACAACCCTTGGGCCTGGGGTGATCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGcc                                       <  1:1626255/150‑1 (MQ=255)
                                                                                                           cACGCGCCACCACACCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAAGCCCTCCCCCAGAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCa                                      <  2:777466/150‑1 (MQ=255)
                                                                                                              gcgcCACCACAACCCTTGGGCCTGGGGTGGCGCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgaccggtctgcgttctgtgccagctcacccaggtggcggttcacatatgccgcgtcgatc                                   <  2:2644180‑M2/150‑110 (MQ=255)
                                                                                                                gcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAAc                                 >  1:931075/1‑150 (MQ=255)
                                                                                                                gcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctg                                 <  1:2022197‑M2/150‑112 (MQ=255)
                                                                                                                gcCACCACAACCCTTGGGCCTGTGGTGATCCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctg                                 <  2:1047240‑M2/150‑112 (MQ=255)
                                                                                                                   accacAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTg                              >  1:910036/1‑150 (MQ=255)
                                                                                                                   accacAACCCTTGGGCCGGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTg                              <  2:1487900/150‑1 (MQ=255)
                                                                                                                     cacaCCCCTGGGGCCTGGGGTGTGCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGca                            <  1:125030/150‑1 (MQ=255)
                                                                                                                     cacaACCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctggatcg                            <  2:3664785‑M2/150‑117 (MQ=255)
                                                                                                                      ccacccgttGGGCGGGGGGTGATCCGGGGGGGGGGGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGcac                           <  1:3276093/143‑1 (MQ=255)
                                                                                                                      acaACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGcac                           >  1:428951/1‑150 (MQ=255)
                                                                                                                        aaCCCTTGGGCCTGTGGTGATCCGTGCGGGGgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctggatcggcg                         <  1:2261477‑M2/150‑120 (MQ=255)
                                                                                                                        aaCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaa                         >  1:3668122/1‑150 (MQ=255)
                                                                                                                        aaCCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaa                         >  2:1843422/1‑150 (MQ=255)
                                                                                                                        aaCCCATGGGCCTGTGGTGATCCGTGCGGGGgcgggcttcccagccaaactgtcgaaaccggtcttacagaaagtaacgcgacaggtcttcgttctttgccagctcacccaggtggctgttcacatacgccgcgacgatctggatcggcg                         <  1:4394824‑M2/150‑120 (MQ=255)
                                                                                                                         aCCCTGGGGCCTGTGGTGATCCGTGGGGGGGGGGGCTGGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACaaa                        <  1:2572226/150‑1 (MQ=255)
                                                                                                                            cTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGcc                     >  2:4456850/1‑150 (MQ=255)
                                                                                                                            cTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGcc                     >  1:2322458/1‑150 (MQ=255)
                                                                                                                                gggttgtgGTGGCGCGTGCAGGAGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCa                 >  1:383158/5‑150 (MQ=38)
                                                                                                                                 gCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCaa                <  2:3082967/150‑1 (MQ=255)
                                                                                                                                 gCCTGTGGTGATCCGTGCAGGAgcgggcttgcccgcgaaactgtcgaaaccggtcttacagaatgtaacgcgacaggtcttcgttctgtgccagctcacccaggtggctgttcacatacgccgcgtcgatctggatcggcgcttcgtcat                >  2:2022197‑M2/1‑22 (MQ=255)
                                                                                                                                  gttgtgGTGGCGCGTGCAGGAGCGGGCTTGCCCGCGAAAAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAg               <  2:1538978/148‑1 (MQ=38)
                                                                                                                                   cTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGa              >  1:2861042/1‑150 (MQ=255)
                                                                                                                                    tgtgGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCga                                                                        >  2:62483/1‑91 (MQ=255)
                                                                                                                                    tgtgGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCga                                                                        <  1:62483/91‑1 (MQ=255)
                                                                                                                                    tgtgGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGAc             >  1:3413308/1‑150 (MQ=255)
                                                                                                                                      tggtgATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGAccc           >  1:219806/1‑150 (MQ=255)
                                                                                                                                      tggtgATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGAccc           <  2:4070991/150‑1 (MQ=255)
                                                                                                                                        gtgATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTc         >  2:1228310/1‑150 (MQ=255)
                                                                                                                                         tgATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCa        >  1:2043220/1‑150 (MQ=255)
                                                                                                                                            tCCGTGGGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGcc     <  2:504372/150‑1 (MQ=255)
                                                                                                                                            tCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGcc     >  2:4576405/1‑150 (MQ=255)
                                                                                                                                             ccGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCACATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGACTCCAAGACCCTCAGCCt    >  2:759415/1‑150 (MQ=255)
                                                                                                                                               gTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGCCTCa  <  2:1096165/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
CGGCGTATGTGAACAGCCACCTGGGTGAGCTGGCACAGAACGAAGACCTGTCGCGTTACATTCTGTAAGACCGGTTTCGACAGTTTCGCGGGCAAGCCCGCTCCTGCACGCGCCACCACAACCCTTGGGCCTGTGGTGATCCGTGCGGGGGCGGGCTTGCCCGCGAAGAGGCCCTTGAATCCCTCCCCCATAAGCTGGAAGCTGTGCTTCATCAGCTCCCGAGAGATTCTGCCCATGGCCCGCCTGCCTACCGCCATCAACCTGCACAAAGCCTCCAAGACCCTCAGCCTCA  >  NC_002947/5700621‑5700912

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