Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,059 Δ1 bp 7.6% intergenic (+35/+37) 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,0590T. 7.6% 210.8 / 11.1 66intergenic (+35/+37)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 T (30/31);  new base . (2/3);  total (32/34)
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.14e-01

CGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAT  >  NC_002947/4749918‑4750208
                                                                                                                                             |                                                                                                                                                       
cgcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGtt                                                                                                                                                        <  1:864731/143‑1 (MQ=255)
cgcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGtt                                                                                                                                                        >  2:864731/1‑143 (MQ=255)
    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                                                                                                                                             >  2:660772/1‑106 (MQ=255)
                                              ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑c                                                                                                                                             <  1:660772/106‑1 (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                                                                                                                                 <  2:707975/67‑1 (MQ=255)
                                                                                                 gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCttttt                                                                                                                                 >  1:707975/1‑67 (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                                                         <  1:1094897/135‑1 (MQ=255)
                                                                                                     aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                         >  2:1094897/1‑135 (MQ=255)
                                                                                                          aGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCgg                                     >  1:733227/1‑150 (MQ=255)
                                                                                                           ggCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGtt                                                        <  1:662629/130‑1 (MQ=255)
                                                                                                           ggCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGtt                                                        >  2:662629/1‑130 (MQ=255)
                                                                                                             ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  2:826836/93‑1 (MQ=255)
                                                                                                             ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           >  1:826836/1‑93 (MQ=255)
                                                                                                              cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                              <  2:1535425/57‑1 (MQ=255)
                                                                                                              cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                              >  1:1535425/1‑57 (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                                                                                                             <  2:940239/73‑1 (MQ=255)
                                                                                                               ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGc                                                                                                             >  1:940239/1‑73 (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                                                         >  1:711315/1‑114 (MQ=255)
                                                                                                                          ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGt                                                         <  2:711315/114‑1 (MQ=255)
                                                                                                                          ggCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGttt                     >  2:1203998/1‑150 (MQ=255)
                                                                                                                           gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                    <  1:234454/150‑1 (MQ=255)
                                                                                                                           gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                    >  2:287691/1‑150 (MQ=255)
                                                                                                                             aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                           >  1:495071/1‑93 (MQ=255)
                                                                                                                             aacaaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                           <  2:495071/93‑1 (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)
                                                                                                                                             |                                                                                                                                                       
CGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAT  >  NC_002947/4749918‑4750208

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