Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,060 Δ1 bp 5.2% intergenic (+36/+36) 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,0600T. 5.2% 310.0 / 9.7 97intergenic (+36/+36)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 (42/50);  new base . (2/3);  total (44/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.39e-01

Only 100 of 106 total aligned reads displayed.

AACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCA  >  NC_002947/4749914‑4750207
                                                                                                                                                   |                                                                                                                                                     
aaCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCCGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    <  1:2114625/150‑1 (MQ=255)
aaCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    >  2:2166907/1‑150 (MQ=255)
aaCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGAATGGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    <  1:2090658/150‑1 (MQ=255)
 aCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCa                                                                                                                                                   <  2:2262231/150‑1 (MQ=255)
     gcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAg                                                                                                                                               >  1:639212/1‑150 (MQ=255)
         aaCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                         <  2:2383086/150‑1 (MQ=255)
           cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                       >  1:619831/1‑150 (MQ=255)
           cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCGACCAAAATTCATCGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                       <  2:663486/150‑1 (MQ=255)
            tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCggggg                                                                                                                                      <  1:2029154/150‑1 (MQ=255)
                tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCttt                                                                                                                                  >  1:982560/1‑150 (MQ=255)
                    cTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgt                                                                                                                              <  1:479653/150‑1 (MQ=255)
                     tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgtt                                                                                                                             <  1:1630548/150‑1 (MQ=255)
                     tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGACCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAA‑‑CCGGGGGCTTTttgtt                                                                                                                             <  1:29668/150‑1 (MQ=38)
                     tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCGGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTttgtt                                                                                                                             <  2:1175627/150‑1 (MQ=38)
                      gCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTttgttg                                                                                                                            <  1:1760857/150‑1 (MQ=255)
                       cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTgg                                                                                                                           <  2:555923/150‑1 (MQ=255)
                          gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTg                                                                                                                        <  2:1191280/150‑1 (MQ=255)
                           aCCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGc                                                                                                                       <  1:984793/150‑1 (MQ=255)
                             cTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCa                                                                                                                     >  1:1931286/1‑150 (MQ=255)
                              tGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAt                                                                                                                    <  2:1931286/150‑1 (MQ=255)
                               gCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAtt                                                                                                                   <  2:338675/150‑1 (MQ=255)
                                 gagaTCAAGGACCCGGCCCAGAACATCCCCTGGGTCACTCCTGGGGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTAt                                                                                                                 <  2:982560/147‑1 (MQ=255)
                                 gACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGc                                                                                                                          <  2:1165840/141‑1 (MQ=255)
                                 gACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGc                                                                                                                          >  1:1165840/1‑141 (MQ=255)
                                    aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                              <  1:1079207/150‑1 (MQ=255)
                                       aaGGACCCGGCCCAGAGCATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                           <  2:1223472/150‑1 (MQ=255)
                                       aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                           >  1:605320/1‑150 (MQ=255)
                                       aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                           >  1:2195208/1‑150 (MQ=255)
                                        aGGCCCCGCCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAg                                                                                                          <  2:1641619/150‑1 (MQ=255)
                                         ggACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCa                                                                                                               <  2:47496/144‑1 (MQ=255)
                                         ggACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCa                                                                                                               >  1:47496/1‑144 (MQ=255)
                                           aCCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTAc                                                                                                       <  1:294717/150‑1 (MQ=255)
                                                  ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    >  1:858970/1‑100 (MQ=255)
                                                  ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    <  2:858970/100‑1 (MQ=255)
                                                  ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTc                                                                                                >  1:1185526/1‑150 (MQ=255)
                                                  ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTc                                                                                                <  2:2422898/150‑1 (MQ=255)
                                                  ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTc                                                                                                >  1:800889/1‑150 (MQ=255)
                                                    aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                         <  1:467344/107‑1 (MQ=255)
                                                    aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCgg                                                                                                                                         >  2:467344/1‑107 (MQ=255)
                                                     gAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTc                                                                                             <  2:639212/150‑1 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGATGGAGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCAGTATCAGGCAGTACGCACTTCCTCaaa                                                                                          <  2:1207410/150‑1 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                          <  2:2195208/150‑1 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                          <  1:1050860/150‑1 (MQ=255)
                                                          tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                        <  2:1185526/150‑1 (MQ=255)
                                                              ccTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                     <  2:467993/150‑1 (MQ=255)
                                                                tGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                   <  2:699316/150‑1 (MQ=255)
                                                                tGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGCTCGCTTCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTGTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCt                                                                                   >  1:1239488/1‑150 (MQ=255)
                                                                     cACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATg                                                                              <  1:2107915/150‑1 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAGCAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                       <  2:2415423/90‑1 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAGCAAAAAGCCCCCGGTTCTCATGAG‑‑CCgggg                                                                                                                                       >  1:2415423/1‑90 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGGTGGCg                                                                           >  1:1541585/1‑150 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                           >  2:546608/1‑150 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                           >  1:1059106/1‑150 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCGGTtctc                                                                                                                                                    <  2:1501371/78‑1 (MQ=255)
                                                                        cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCGGTtctc                                                                                                                                                    >  1:1501371/1‑78 (MQ=255)
                                                                         ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                          >  2:1897581/1‑150 (MQ=255)
                                                                           tGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATc                                                                        >  2:2065097/1‑150 (MQ=255)
                                                                                gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                   <  1:2065097/150‑1 (MQ=255)
                                                                                   cTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCt                                                                <  2:356373/150‑1 (MQ=255)
                                                                                    tGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                               >  1:831818/1‑150 (MQ=255)
                                                                                    tGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                               <  2:605320/150‑1 (MQ=255)
                                                                                       cGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAGGAGTTCTGGATGGCGATCCCGCCGCTTGAg                                                            >  2:786191/1‑150 (MQ=255)
                                                                                        gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    >  1:938765/1‑63 (MQ=255)
                                                                                        gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                    <  2:938765/63‑1 (MQ=255)
                                                                                        gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCAc                                                                                                   <  2:249999/110‑1 (MQ=255)
                                                                                        gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCAc                                                                                                   >  1:249999/1‑110 (MQ=255)
                                                                                              aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                          >  2:2277956/1‑129 (MQ=255)
                                                                                              aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                          <  1:2277956/129‑1 (MQ=255)
                                                                                                cTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCAc                                                   <  1:1678732/150‑1 (MQ=255)
                                                                                                     cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACg                                              >  2:2133917/1‑150 (MQ=255)
                                                                                                          aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACt                                         >  2:970355/1‑150 (MQ=255)
                                                                                                           tGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcag                                                                     >  1:1758547/1‑121 (MQ=255)
                                                                                                           tGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcag                                                                     <  2:1758547/121‑1 (MQ=255)
                                                                                                            gTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                             <  1:1262530/112‑1 (MQ=255)
                                                                                                            gTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTgg                                       >  1:742704/1‑150 (MQ=255)
                                                                                                            gTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCGTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAATAGTTCTGGATgg                                                                             >  2:1262530/1‑112 (MQ=255)
                                                                                                                    ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAg                                                            >  1:83729/1‑121 (MQ=255)
                                                                                                                    ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAg                                                            <  2:83729/121‑1 (MQ=255)
                                                                                                                    ccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtgg                               >  1:502855/1‑150 (MQ=255)
                                                                                                                     cTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCGGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtggt                              >  2:1493846/1‑150 (MQ=255)
                                                                                                                        cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                           >  1:371905/1‑150 (MQ=255)
                                                                                                                          cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                             <  2:2347456/98‑1 (MQ=255)
                                                                                                                          cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                             >  1:2347456/1‑98 (MQ=255)
                                                                                                                          cTTGCGGCAACAAAAAGCAACCGGGTGTGGTGTG‑‑CCGGGGGGTGTTTGGTGGAAACCAGTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                         >  1:1872716/1‑150 (MQ=255)
                                                                                                                                gCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                   <  2:433830/150‑1 (MQ=255)
                                                                                                                                ccaacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                     >  2:419706/2‑84 (MQ=18)
                                                                                                                                ccaacaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGtt                                                                                     <  1:419706/83‑1 (MQ=18)
                                                                                                                                   acaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                  <  2:1422886/148‑1 (MQ=38)
                                                                                                                                   acaaAAAGCCCCCGG‑‑CTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                  >  1:1422886/1‑148 (MQ=38)
                                                                                                                                    caaAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGTTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATa               >  1:2194441/1‑150 (MQ=255)
                                                                                                                                     aaaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                <  1:184405/132‑1 (MQ=255)
                                                                                                                                     aaaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                >  2:184405/1‑132 (MQ=255)
                                                                                                                                      aaaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAAc             <  2:800889/150‑1 (MQ=255)
                                                                                                                                       aaaGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACg            >  2:709875/1‑150 (MQ=255)
                                                                                                                                           cccccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGgag        >  1:1907591/1‑150 (MQ=255)
                                                                                                                                             cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg           <  1:983523/145‑1 (MQ=255)
                                                                                                                                             cccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa      >  1:1822364/1‑150 (MQ=255)
                                                                                                                                              ccGNTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTCGCCGTTTTCATAACGGGAGAAc     >  2:1504054/1‑150 (MQ=255)
                                                                                                                                              ccGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAAc     <  2:831818/150‑1 (MQ=255)
                                                                                                                                                 gTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCa  <  2:619831/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAG‑‑CCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCA  >  NC_002947/4749914‑4750207

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: