Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,751,986 A→G 5.8% intergenic (+40/+32) 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,9860AG 5.8% 489.8 / 12.2 158intergenic (+40/+32)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 A (71/76);  new base G (2/7);  total (74/84)
Fisher's exact test for biased strand distribution p-value = 1.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 165 total aligned reads displayed.

ACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGACC  >  NC_002947/4751837‑4752135
                                                                                                                                                     |                                                                                                                                                     
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)
   cgcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGa                                                                                                                                                    <  1:225993/150‑1 (MQ=255)
         aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                            <  2:2577320/150‑8 (MQ=255)
         aCTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGCCAACAAAAAGCCCCCGG‑‑CTCATGAgaaccggg                                                                                                                                            <  1:4703790/150‑8 (MQ=255)
          cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCgggg                                                                                                                                             >  1:2671403/1‑150 (MQ=255)
          cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCgggg                                                                                                                                             <  2:4984499/150‑1 (MQ=255)
           tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCggggg                                                                                                                                            <  1:2039954/150‑1 (MQ=255)
            gCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGc                                                                                                                                           <  1:1126863/150‑1 (MQ=255)
               tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCttt                                                                                                                                        >  2:1377484/1‑150 (MQ=255)
               tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCttt                                                                                                                                        <  1:2128962/150‑1 (MQ=39)
               tGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCttt                                                                                                                                        <  1:2480529/150‑1 (MQ=39)
                gCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCtttt                                                                                                                                       <  2:3566405/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:5378284/150‑1 (MQ=255)
                    tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTttgtt                                                                                                                                   <  1:3887883/150‑1 (MQ=255)
                    tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGCTCATGAGAACCGGGGGCtttt                                                                                                                                       >  1:1974225/1‑146 (MQ=39)
                      cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                                 >  2:403301/1‑150 (MQ=255)
                      cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                                 <  1:5435979/150‑1 (MQ=255)
                        aGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCt                                                                                                                               <  2:3187512/150‑1 (MQ=255)
                                   aTCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                                    <  2:683007/150‑1 (MQ=255)
                                    tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   >  1:945520/1‑150 (MQ=255)
                                    tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                   >  1:1192337/1‑150 (MQ=255)
                                      aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCa                                                                                                                 >  2:2527157/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)
                                            ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTGTTGTTGGCTGCCATTCTCAGGCAGTACGc                                                                                                           >  2:2945764/1‑150 (MQ=255)
                                                   aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCggggg                                                                                                                                            <  1:4113204/110‑1 (MQ=255)
                                                   aGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCggggg                                                                                                                                            >  2:4113204/1‑110 (MQ=255)
                                                    gAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTc                                                                                                   >  2:3431760/1‑150 (MQ=255)
                                                     aaCATCACCTGCGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAATTAGCCCCCGGTTCTGATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTGCTCa                                                                                                  >  2:145752/1‑150 (MQ=255)
                                                       cATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                                <  2:4204175/150‑1 (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)
                                                                 gTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGa                                                                                      <  1:3465145/150‑1 (MQ=255)
                                                                  tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGAt                                                                                     >  1:5225996/1‑150 (MQ=255)
                                                                  tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTAAGCACTTCCTCAAAGAGTGCTGGAt                                                                                     >  1:5225881/1‑150 (MQ=255)
                                                                   cACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATg                                                                                    <  1:243091/150‑1 (MQ=255)
                                                                    aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                                   <  1:1089235/150‑1 (MQ=255)
                                                                      cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCg                                                                                 >  1:5253353/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                                                                            >  2:4742272/1‑150 (MQ=255)
                                                                           aaGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCca                                                                            >  1:4250465/1‑150 (MQ=255)
                                                                              gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                         <  1:340796/150‑1 (MQ=255)
                                                                                  tGGCGGGCCGAACTATCCGACCATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCTTGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCtt                                                                     <  2:1165254/150‑1 (MQ=255)
                                                                                    gCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                                   <  1:1048902/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)
                                                                                                 tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCaa                                                      >  2:4806898/1‑150 (MQ=255)
                                                                                                 tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCaa                                                      <  1:2704088/150‑1 (MQ=255)
                                                                                                 tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCaa                                                      >  2:2069837/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)
                                                                                                                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)
                                                                                                                    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:5402161/150‑1 (MQ=255)
                                                                                                                        cTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGc                               >  1:1925211/1‑150 (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:3197533/1‑150 (MQ=255)
                                                                                                                              gCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTGTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAATGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                         >  2:4029895/1‑150 (MQ=255)
                                                                                                                                 acaaAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAAc                                                     <  1:2819497/119‑1 (MQ=255)
                                                                                                                                        gCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAACGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGga               <  2:4217695/150‑1 (MQ=255)
                                                                                                                                        gCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGga               >  2:1921197/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:2082881/1‑150 (MQ=255)
                                                                                                                                           cccGGTTCTCATGAGCCGGGGGCTTGTTGGTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGaa            >  2:3725410/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)
                                                                                                                                                     |                                                                                                                                                     
ACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGCCCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGACC  >  NC_002947/4751837‑4752135

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: