Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,068:2 +A 5.7% intergenic (+44/+28) PP_4203 → / ← mqsA ubiquinone oxidoreductase electron transfer flavoprotein/type II toxin‑antitoxin MqsRA system antitoxin/DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,750,0682.A 5.7% 292.9 / 10.8 70intergenic (+44/+28)PP_4203/mqsAubiquinone oxidoreductase electron transfer flavoprotein/type II toxin‑antitoxin MqsRA system antitoxin/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (32/34);  new base A (2/2);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

GAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGACC  >  NC_002947/4749922‑4750213
                                                                                                                                                    |                                                                                                                                                 
gAACTGCGTGCACTGCAAGACTTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGCCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccg                                                                                                                                                <  1:1013444/150‑6 (MQ=255)
         gCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTCAAGGCGCTGCCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCtttt                                                                                                                                       <  2:1205327/150‑1 (MQ=38)
            cTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgt                                                                                                                                    <  2:1283254/150‑1 (MQ=255)
             tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                                   <  1:1466214/150‑1 (MQ=255)
               cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTgg                                                                                                                                 <  2:58664/150‑1 (MQ=255)
               cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGAGCAGGTCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTgg                                                                                                                                 <  1:1335455/150‑1 (MQ=255)
                           cATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCa                                                                                                                     <  2:960863/150‑1 (MQ=255)
                            aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                                    >  2:413069/1‑150 (MQ=255)
                               aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCAGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                                 <  1:355604/150‑1 (MQ=255)
                                     ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGc                                                                                                           >  2:397863/1‑150 (MQ=255)
                                       ggCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCAc                                                                                                         >  2:234454/1‑150 (MQ=255)
                                          ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑c                                                                                                                                                  <  1:660772/106‑1 (MQ=255)
                                          ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑c                                                                                                                                                  >  2:660772/1‑106 (MQ=255)
                                            aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCt                                                                                                    <  1:413069/150‑1 (MQ=255)
                                                     ccTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                           >  2:1408772/1‑150 (MQ=255)
                                                       tGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                         <  1:1408772/150‑1 (MQ=255)
                                                          gTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGa                                                                                      <  2:1354262/150‑1 (MQ=255)
                                                             aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   <  1:330315/118‑1 (MQ=255)
                                                             aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   >  2:330315/1‑118 (MQ=255)
                                                               cccTGAAGGCGCTGGCGGGCCGAACTATGCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CGGGGGGCTTTTTGTTGGCTGCCATTAGCTGGCAGTACACACTTCCTCAAAGAGTACTgtagtggg                                                                                 >  2:671549/1‑143 (MQ=255)
                                                               cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                            >  2:1266949/1‑91 (MQ=255)
                                                               cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                            <  1:1266949/91‑1 (MQ=255)
                                                                 cTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCACCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGAt                                                                               >  1:476446/1‑150 (MQ=255)
                                                                      ggCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCAACTATCAGCCAACATGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagc                                                                          >  2:1503222/1‑150 (MQ=255)
                                                                           tGGCGGGCCGAACTATCTAAACGTGTAAGGCCCCTGCCCTTGCGGCAACAAACAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCATGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                     <  1:1318213/150‑1 (MQ=255)
                                                                           tGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGAGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                     <  2:1134728/150‑1 (MQ=255)
                                                                               gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAAGCACTACGCGTCTGCTCGATGAATTTGTTATGGCGATCCAGCAGCTTGAGc                                                                 >  1:1598889/1‑150 (MQ=255)
                                                                                  ccGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                              >  2:653683/1‑150 (MQ=255)
                                                                                            cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                   <  2:1452996/135‑1 (MQ=255)
                                                                                            cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                   >  1:1452996/1‑135 (MQ=255)
                                                                                            cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCCCCAACg                                                    >  2:613170/1‑150 (MQ=255)
                                                                                             gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCttttt                                                                                                                                      >  1:707975/1‑67 (MQ=255)
                                                                                             gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCttttt                                                                                                                                      <  2:707975/67‑1 (MQ=255)
                                                                                             gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACt                                                                                                        <  2:1445967/97‑1 (MQ=255)
                                                                                             gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACt                                                                                                        >  1:1445967/1‑97 (MQ=255)
                                                                                                cATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCAc                                                <  1:653683/150‑1 (MQ=255)
                                                                                                 aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                              >  2:1094897/1‑135 (MQ=255)
                                                                                                 aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                              <  1:1094897/135‑1 (MQ=255)
                                                                                                      aGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCgg                                          >  1:733227/1‑150 (MQ=255)
                                                                                                       ggCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGtt                                                             >  2:662629/1‑130 (MQ=255)
                                                                                                       ggCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGtt                                                             <  1:662629/130‑1 (MQ=255)
                                                                                                         ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                >  1:826836/1‑93 (MQ=255)
                                                                                                         ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                <  2:826836/93‑1 (MQ=255)
                                                                                                          cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                                   >  1:1535425/1‑57 (MQ=255)
                                                                                                          cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                                   <  2:1535425/57‑1 (MQ=255)
                                                                                                          cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGc                                                                                                                  >  1:278571/1‑74 (MQ=255)
                                                                                                          cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGc                                                                                                                  <  2:278571/74‑1 (MQ=255)
                                                                                                          cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                      <  1:1084414/150‑1 (MQ=255)
                                                                                                           ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGc                                                                                                                  >  1:940239/1‑73 (MQ=255)
                                                                                                           ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGc                                                                                                                  <  2:940239/73‑1 (MQ=255)
                                                                                                             tGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTc                                   >  1:1413579/1‑150 (MQ=255)
                                                                                                              gCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGAAGTTTCACCAACGCCACTGGCGGGGTGGTCt                                  >  1:83033/1‑150 (MQ=255)
                                                                                                                ccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTg                                <  1:397863/150‑1 (MQ=255)
                                                                                                                      ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                              <  2:711315/114‑1 (MQ=255)
                                                                                                                      ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                              >  1:711315/1‑114 (MQ=255)
                                                                                                                      ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGttt                          >  2:1203998/1‑150 (MQ=255)
                                                                                                                       gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                         >  2:287691/1‑150 (MQ=255)
                                                                                                                       gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                         <  1:234454/150‑1 (MQ=255)
                                                                                                                         aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                                <  2:495071/93‑1 (MQ=255)
                                                                                                                         aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                                >  1:495071/1‑93 (MQ=255)
                                                                                                                            aaaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATaa                    <  2:733227/150‑1 (MQ=255)
                                                                                                                              aaaGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                         >  1:1348678/1‑79 (MQ=18)
                                                                                                                              aaaGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                         <  2:1348678/79‑1 (MQ=18)
                                                                                                                               aaGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGCCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg                 >  1:1205868/1‑150 (MQ=38)
                                                                                                                                    cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa            >  2:1347002/1‑150 (MQ=255)
                                                                                                                                        gTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCa        <  2:476446/150‑1 (MQ=255)
                                                                                                                                         ttCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGTCATGGTCATCCAGGAGGCACTGCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAt       <  2:1413579/150‑1 (MQ=255)
                                                                                                                                            tcATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTGCTGGATCGCCATc                                                                              <  1:952278/76‑1 (MQ=255)
                                                                                                                                             cATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCACGCAGTACGCACTTCCTCAAAGAGTTCTGGATCGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAc   <  1:671549/150‑1 (MQ=255)
                                                                                                                                              aTGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                         >  1:74186/1‑63 (MQ=255)
                                                                                                                                              aTGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                         <  2:74186/63‑1 (MQ=255)
                                                                                                                                              aTGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAcc  >  1:164318/1‑150 (MQ=255)
                                                                                                                                              aTGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACAAACGCCACTGGCGGGGTGTCCTGGCCGGTTTCATAACGGGAGAACGCATTGAcc  >  2:129621/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGACC  >  NC_002947/4749922‑4750213

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: