Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,751,984 T→A 5.6% intergenic (+38/+34) 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,751,9840TA 5.6% 548.9 / 14.5 161intergenic (+38/+34)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 (76/76);  new base A (2/7);  total (78/83)
Fisher's exact test for biased strand distribution p-value = 1.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 164 total aligned reads displayed.

AACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGA  >  NC_002947/4751836‑4752133
                                                                                                                                                    |                                                                                                                                                     
aaCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTtctc                                                                                                                                                      <  1:214715/150‑1 (MQ=255)
 aCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAATGCCCCCGGTTCTCa                                                                                                                                                     >  2:1487024/1‑150 (MQ=255)
 aCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCa                                                                                                                                                     >  2:1487201/1‑150 (MQ=255)
 aCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCa                                                                                                                                                     >  2:1487062/1‑150 (MQ=255)
  cgcgCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCAt                                                                                                                                                    <  2:2473629/150‑1 (MQ=255)
          aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                          <  2:2577320/150‑8 (MQ=255)
          aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGCCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                          <  1:4703790/150‑8 (MQ=255)
           cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCgggg                                                                                                                                           <  2:4984499/150‑1 (MQ=255)
             gCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGc                                                                                                                                         <  1:1126863/150‑1 (MQ=255)
               gTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCtt                                                                                                                                       <  1:274809/150‑1 (MQ=255)
                tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttt                                                                                                                                      >  2:1377484/1‑150 (MQ=255)
                tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCttt                                                                                                                                      <  1:2480529/150‑1 (MQ=39)
                 gCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCtttt                                                                                                                                     <  2:3566405/150‑1 (MQ=255)
                  cACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttttt                                                                                                                                    <  1:3483420/150‑1 (MQ=255)
                  cACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttttt                                                                                                                                    <  1:1360877/150‑1 (MQ=255)
                  cACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttttt                                                                                                                                    <  1:2318640/150‑1 (MQ=255)
                    cTGCAAGGCCTGCGACTTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgt                                                                                                                                  <  1:5378231/150‑1 (MQ=255)
                    cTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgt                                                                                                                                  <  1:5378288/150‑1 (MQ=255)
                     tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgtt                                                                                                                                 <  1:3887883/150‑1 (MQ=255)
                     tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCtttt                                                                                                                                     >  1:1974225/1‑146 (MQ=39)
                     tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCTTTttgtt                                                                                                                                 <  1:2088628/150‑1 (MQ=38)
                       cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                               <  1:5435979/150‑1 (MQ=255)
                       cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                               >  2:403301/1‑150 (MQ=255)
                                    aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                                  <  2:683007/150‑1 (MQ=255)
                                     tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                 >  1:1192337/1‑150 (MQ=255)
                                     tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                 >  1:945520/1‑150 (MQ=255)
                                     tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGACCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                 >  1:2436270/1‑150 (MQ=255)
                                       aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                               >  1:2036260/1‑150 (MQ=255)
                                       aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                               >  2:2527157/1‑150 (MQ=255)
                                             ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGc                                                                                                         >  1:3143520/1‑150 (MQ=255)
                                             ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGc                                                                                                         >  2:2327552/1‑150 (MQ=255)
                                                gCCCAGAACAACACCTGGGTCACCCCTGAAGGCGCGGGCGGGCCGAACTATCCCAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTATTGTTGGCTGCCATTATCAGGCAGTACGCACt                                                                                                      <  2:3885465/150‑1 (MQ=255)
                                                    aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCggggg                                                                                                                                          <  1:4113204/110‑1 (MQ=255)
                                                      aaCATCACCTGCGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAATTAGCCCCCGGTTCTGATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTGCTCa                                                                                                >  2:145752/1‑150 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  1:733286/150‑1 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  1:733270/150‑1 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  1:118711/150‑1 (MQ=255)
                                                        cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                              <  1:2041625/150‑1 (MQ=255)
                                                          tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                            <  1:475244/150‑1 (MQ=255)
                                                          tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                            >  1:1390624/1‑150 (MQ=255)
                                                            aCCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCGGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCTTGAGCCGGGGGGTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCGTCAAAGAGt                                                                                          >  2:3299624/1‑150 (MQ=255)
                                                                   tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTAAGCACTTCCTCAAAGAGTGCTGGAt                                                                                   >  1:5225881/1‑150 (MQ=255)
                                                                     aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                                 <  1:1089235/150‑1 (MQ=255)
                                                                       cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                               >  1:3338103/1‑150 (MQ=255)
                                                                       cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                               >  1:1130646/1‑150 (MQ=255)
                                                                       cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                               >  1:5253353/1‑150 (MQ=255)
                                                                        ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                              >  2:5402161/1‑150 (MQ=255)
                                                                            aaGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCca                                                                          >  2:4742272/1‑150 (MQ=255)
                                                                            aaGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCca                                                                          >  1:5246279/1‑150 (MQ=255)
                                                                            aaGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCca                                                                          >  1:4250465/1‑150 (MQ=255)
                                                                              ggCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagc                                                                        <  1:3262732/150‑1 (MQ=255)
                                                                               gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                       <  1:340796/150‑1 (MQ=255)
                                                                                   tGGCGGGCCGAACTATCCGACCATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCTTGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                   <  2:1165254/150‑1 (MQ=255)
                                                                                     gCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                 <  1:1048902/150‑1 (MQ=255)
                                                                                     gCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                 <  1:2684347/150‑1 (MQ=255)
                                                                                            gAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGttt                                                          >  1:4669074/1‑150 (MQ=255)
                                                                                             aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTc                                                         <  2:4199713/150‑1 (MQ=255)
                                                                                              aCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCa                                                        <  1:3487802/150‑1 (MQ=255)
                                                                                                tATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGGTGCGATTATCAGGCAGTACGCAGTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCTcc                                                      >  2:4111821/1‑150 (MQ=255)
                                                                                                  tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCaa                                                    >  2:2069837/1‑150 (MQ=255)
                                                                                                  tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCaa                                                    >  2:4806898/1‑150 (MQ=255)
                                                                                                     gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATCGCCATCcagaa                                                                       <  2:3021572/128‑3 (MQ=255)
                                                                                                       aCATGTAAGGCCCCTGCCCTTGCTTCAACAAAAAGCCCCCGGTTATCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCa                                               <  1:4299791/150‑1 (MQ=255)
                                                                                                         aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTGTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACt                                             >  2:2532126/1‑150 (MQ=255)
                                                                                                             aaGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCg                                         <  1:902626/150‑1 (MQ=255)
                                                                                                               ggCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCggg                                       <  2:4236168/150‑1 (MQ=255)
                                                                                                                 ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggt                                     >  2:4950780/1‑150 (MQ=255)
                                                                                                                  cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                    <  1:762905/150‑1 (MQ=255)
                                                                                                                  cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                    >  1:5117851/1‑150 (MQ=255)
                                                                                                                  cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATGATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                    <  1:4340975/150‑1 (MQ=255)
                                                                                                                  cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTATTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                    <  2:606562/150‑1 (MQ=255)
                                                                                                                   ccTGCCCTTGCGGCAACAAACAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTCTCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCCCGCCAGTGGCgttggtga                                  <  2:4156113/150‑8 (MQ=255)
                                                                                                                     tGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttt                                                                                                                                      >  2:1799484/1‑49 (MQ=255)
                                                                                                                     tGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttt                                                                                                                                      <  1:1799484/49‑1 (MQ=255)
                                                                                                                       cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                               <  1:131954/150‑1 (MQ=255)
                                                                                                                        ccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTg                              >  2:1603951/1‑150 (MQ=255)
                                                                                                                         cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                             >  1:1925211/1‑150 (MQ=255)
                                                                                                                         cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                             <  1:5402161/150‑1 (MQ=255)
                                                                                                                          ttGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGcc                            <  2:3079639/150‑1 (MQ=255)
                                                                                                                            gCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGt                          <  2:3609588/150‑1 (MQ=255)
                                                                                                                             cGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtt                         <  1:403301/150‑1 (MQ=255)
                                                                                                                             agccaacaaAAAGCCCCCGGCTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtt                         >  1:2815176/4‑150 (MQ=38)
                                                                                                                               gCAACAACAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACGGGCGGGGTGGTCTTGCCGtttt                       <  2:771182/150‑1 (MQ=255)
                                                                                                                               gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                               >  2:107161/1‑142 (MQ=255)
                                                                                                                               gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                               <  1:107161/142‑1 (MQ=255)
                                                                                                                               gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                       >  1:2269458/1‑150 (MQ=255)
                                                                                                                                  acaaAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAAc                                                   >  2:2819497/1‑119 (MQ=255)
                                                                                                                                     aaaaGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAAc                 >  1:286866/1‑150 (MQ=255)
                                                                                                                                      aaaGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACg                >  1:5012159/1‑150 (MQ=255)
                                                                                                                                         gCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCACAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGga             <  2:4217647/150‑1 (MQ=255)
                                                                                                                                         gCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGga             >  2:1921197/1‑150 (MQ=255)
                                                                                                                                            cccGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                      <  2:2982826/138‑1 (MQ=255)
                                                                                                                                            cccGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTc                      >  1:2982826/1‑138 (MQ=255)
                                                                                                                                            cccGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa          <  2:1201559/150‑1 (MQ=255)
                                                                                                                                            cccGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa          >  1:2082881/1‑150 (MQ=255)
                                                                                                                                                gTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCa      <  2:2205906/150‑1 (MQ=255)
                                                                                                                                                 ttCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAt     >  2:154434/1‑150 (MQ=255)
                                                                                                                                                 ttCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAt     <  1:2359671/150‑1 (MQ=255)
                                                                                                                                                   ctcATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTg   <  1:1688473/150‑1 (MQ=255)
                                                                                                                                                    tcATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGa  >  2:1697725/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGA  >  NC_002947/4751836‑4752133

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