Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,751,987 T→A 5.7% intergenic (+41/+31) 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,9870TA 5.7% 528.7 / 15.4 159intergenic (+41/+31)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 (74/76);  new base A (2/7);  total (76/83)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 160 total aligned reads displayed.

CGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGACC  >  NC_002947/4751838‑4752135
                                                                                                                                                     |                                                                                                                                                    
cgcgCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCAt                                                                                                                                                      <  2:2473629/150‑1 (MQ=255)
        aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCggg                                                                                                                                              <  2:3003964/150‑1 (MQ=255)
        aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                            <  2:2577320/150‑8 (MQ=255)
        aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGCCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                            <  1:4703790/150‑8 (MQ=255)
           gCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGc                                                                                                                                           <  1:1126863/150‑1 (MQ=255)
             gTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCtt                                                                                                                                         <  1:274809/150‑1 (MQ=255)
              tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCttt                                                                                                                                        <  2:1696771/150‑1 (MQ=39)
              tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCttt                                                                                                                                        <  1:2480529/150‑1 (MQ=39)
              tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCttt                                                                                                                                        <  1:2128962/150‑1 (MQ=39)
               gCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCtttt                                                                                                                                       <  2:3566405/150‑1 (MQ=255)
                  cTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgt                                                                                                                                    <  1:5378284/150‑1 (MQ=255)
                  cTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgt                                                                                                                                    <  1:5378288/150‑1 (MQ=255)
                   tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgtt                                                                                                                                   >  1:5364234/1‑150 (MQ=255)
                   tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgtt                                                                                                                                   <  1:3887883/150‑1 (MQ=255)
                   tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCTTTttgtt                                                                                                                                   <  1:2088628/150‑1 (MQ=38)
                    gCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgttg                                                                                                                                  >  1:1332123/1‑150 (MQ=255)
                     cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                                 <  1:5435979/150‑1 (MQ=255)
                     cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                                 >  2:403301/1‑150 (MQ=255)
                                  aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                                    <  2:683007/150‑1 (MQ=255)
                                   tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   >  1:945520/1‑150 (MQ=255)
                                   tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   >  1:1192337/1‑150 (MQ=255)
                                   tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGACCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   >  1:2436270/1‑150 (MQ=255)
                                     aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGcc                                                                                                                 >  1:196878/1‑149 (MQ=255)
                                     aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                                 >  1:2036260/1‑150 (MQ=255)
                                     aaGGACCCGGCCCAGAACATCACCAGGGTCACCACTCAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCAGGGGGCTTTTTGTTGGCTGCCATTATAAGGCa                                                                                                                 >  1:4825689/1‑150 (MQ=255)
                                           ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGc                                                                                                           >  2:2327552/1‑150 (MQ=255)
                                           ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGc                                                                                                           >  1:3143520/1‑150 (MQ=255)
                                              gCCCAGAACAACACCTGGGTCACCCCTGAAGGCGCGGGCGGGCCGAACTATCCCAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTATTGTTGGCTGCCATTATCAGGCAGTACGCACt                                                                                                        <  2:3885465/150‑1 (MQ=255)
                                                 cAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTcc                                                                                                     >  1:5150099/1‑150 (MQ=255)
                                                  aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCggggg                                                                                                                                            <  1:4113204/110‑1 (MQ=255)
                                                    aaCATCACCTGCGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAATTAGCCCCCGGTTCTGATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTGCTCa                                                                                                  >  2:145752/1‑150 (MQ=255)
                                                      cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                <  2:1457573/150‑1 (MQ=255)
                                                      cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                <  1:733282/150‑1 (MQ=255)
                                                      cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                <  1:118711/150‑1 (MQ=255)
                                                      cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                <  1:2041625/150‑1 (MQ=255)
                                                        tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAaga                                                                                              >  1:1390624/1‑150 (MQ=255)
                                                         cACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAagag                                                                                             <  1:1766157/150‑1 (MQ=255)
                                                          aCCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCGGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCTTGAGCCGGGGGGTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCGTCAAAGAGt                                                                                            >  2:3299624/1‑150 (MQ=255)
                                                                gTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGa                                                                                      <  1:3465145/150‑1 (MQ=255)
                                                                 tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGAt                                                                                     >  1:5225996/1‑150 (MQ=255)
                                                                 tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGAt                                                                                     >  1:4403106/1‑150 (MQ=255)
                                                                   aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                                   <  1:1089235/150‑1 (MQ=255)
                                                                     cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                                 >  1:743603/1‑150 (MQ=255)
                                                                     cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                                 >  2:5358085/1‑150 (MQ=255)
                                                                     cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                                 >  1:1130646/1‑150 (MQ=255)
                                                                     cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                                 >  1:3338103/1‑150 (MQ=255)
                                                                      ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGa                                                                                >  2:5402161/1‑150 (MQ=255)
                                                                          aaGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCca                                                                            >  1:5246279/1‑150 (MQ=255)
                                                                            ggCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagc                                                                          <  1:3262732/150‑1 (MQ=255)
                                                                             gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                         <  1:340796/150‑1 (MQ=255)
                                                                                   gCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                   <  1:1048902/150‑1 (MQ=255)
                                                                                   gCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                   <  1:2684347/150‑1 (MQ=255)
                                                                                     gggCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGGTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCACCGAGTTCTGGATGGCGATCCAGCAGCTTGAGc                                                                 <  1:942831/150‑1 (MQ=255)
                                                                                           aaCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTc                                                           <  1:2327552/150‑1 (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)
                                                                                                     aCATGTAAGGCCCCTGCCCTTGCTTCAACAAAAAGCCCCCGGTTATCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCa                                                 <  1:4299791/150‑1 (MQ=255)
                                                                                                       aTGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTGTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACt                                               >  2:2532126/1‑150 (MQ=255)
                                                                                                             ggCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCggg                                         <  2:4236168/150‑1 (MQ=255)
                                                                                                               ccccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggt                                       >  2:4950780/1‑150 (MQ=255)
                                                                                                                cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                      >  1:5117851/1‑150 (MQ=255)
                                                                                                                cccTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtg                                      <  1:762905/150‑1 (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)
                                                                                                                     cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                                 <  1:131954/150‑1 (MQ=255)
                                                                                                                     cccTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                                 <  2:2475306/150‑1 (MQ=255)
                                                                                                                       cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                               <  1:5402196/150‑1 (MQ=255)
                                                                                                                       cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                               >  1:1925211/1‑150 (MQ=255)
                                                                                                                       cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                               >  1:4207469/1‑150 (MQ=255)
                                                                                                                        ttGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGcc                              <  2:3079639/150‑1 (MQ=255)
                                                                                                                          gCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGt                            <  2:3609588/150‑1 (MQ=255)
                                                                                                                           agccaacaaAAAGCCCCCGGCTCATGAGAACCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtt                           >  1:2815176/4‑150 (MQ=38)
                                                                                                                             gCAACAACAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACGGGCGGGGTGGTCTTGCCGtttt                         <  2:771182/150‑1 (MQ=255)
                                                                                                                             gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                                 <  1:107161/142‑1 (MQ=255)
                                                                                                                             gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCtt                                 >  2:107161/1‑142 (MQ=255)
                                                                                                                             gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                         >  1:2269458/1‑150 (MQ=255)
                                                                                                                                acaaAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAAc                                                     >  2:2819497/1‑119 (MQ=255)
                                                                                                                                   aaaaGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAAc                   >  1:286866/1‑150 (MQ=255)
                                                                                                                                    aaaGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATACCg                  >  2:4769408/1‑150 (MQ=255)
                                                                                                                                     aaGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg                 <  2:1192337/150‑1 (MQ=255)
                                                                                                                                       gCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCACAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGga               <  2:4217647/150‑1 (MQ=255)
                                                                                                                                       gCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGga               >  2:1921197/1‑150 (MQ=255)
                                                                                                                                         ccccGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGgaga             >  1:4722319/1‑150 (MQ=255)
                                                                                                                                         ccccGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGgaga             >  1:3914756/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:2082739/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)
                                                                                                                                                    aTGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAcc  >  2:4340292/1‑150 (MQ=255)
                                                                                                                                                    aTGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAcc  >  2:3304347/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGACC  >  NC_002947/4751838‑4752135

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