Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,017 T→A 5.4% 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,750,0170TA 5.4% 234.1 / 2.7 74intergenic (+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 (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

ACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAT  >  NC_002947/4749870‑4750163
                                                                                                                                                      |                                                                                                                                                  
aCGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTtct                                                                                                                                                    >  2:119276/1‑148 (MQ=255)
    gcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCAt                                                                                                                                                 >  3:297756/1‑147 (MQ=255)
    gcAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCAt                                                                                                                                                 >  4:325675/1‑147 (MQ=255)
     cAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAg                                                                                                                                              <  1:42630/149‑1 (MQ=255)
          cTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCggg                                                                                                                                         <  4:297756/149‑1 (MQ=255)
           tGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCgggg                                                                                                                                        >  6:367377/1‑149 (MQ=255)
              gTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTCGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATCAGCCGGGGGCt                                                                                                                                     >  6:357368/1‑149 (MQ=255)
                  aCTGCAAGACGTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCttttt                                                                                                                                 <  2:104782/149‑1 (MQ=255)
                    tGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTttgt                                                                                                                               >  2:16084/1‑149 (MQ=255)
                      cAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTttgtt                                                                                                                              <  3:77343/148‑1 (MQ=255)
                         gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGa                                                                                                                                               >  2:84243/1‑128 (MQ=255)
                         gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGa                                                                                                                                               <  1:84243/128‑1 (MQ=255)
                         gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCt                                                                                                                          <  5:367377/149‑1 (MQ=255)
                            cTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGCTCATGAGACCCGGGGGCTTTttg                                                                                                                                >  6:252229/1‑140 (MQ=38)
                            cTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGCTCATGAGAACCGGGGGCTTTttg                                                                                                                                <  5:252229/140‑1 (MQ=38)
                             tGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCa                                                                                                                      >  6:320226/1‑149 (MQ=255)
                               cGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCAt                                                                                                                     >  6:106226/1‑148 (MQ=255)
                               cGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAACC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCa                                                                                                                      <  5:247751/147‑1 (MQ=255)
                                gACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGACCATGTAACGCCCCTGCCCTCGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTa                                                                                                                   <  5:266951/149‑1 (MQ=255)
                                    tCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAg                                                                                                               >  2:209606/1‑149 (MQ=255)
                                     cAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                              >  5:296096/1‑149 (MQ=255)
                                      aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCggg                                                                                                                                         <  2:290818/121‑1 (MQ=255)
                                      aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCggg                                                                                                                                         >  1:290818/1‑121 (MQ=255)
                                      aaGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGc                                                                                                             <  7:22760/149‑1 (MQ=255)
                                        ggACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAg                                                                                                           >  2:205863/1‑149 (MQ=255)
                                          aCCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCTGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGt                                                                                                          >  4:41420/1‑148 (MQ=255)
                                            ccGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCAAGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACg                                                                                                       >  6:76163/1‑149 (MQ=255)
                                              ggCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCa                                                                                                     >  8:218975/1‑149 (MQ=255)
                                                 ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                            <  3:51014/123‑1 (MQ=255)
                                                 ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTgg                                                                                                                            >  4:51014/1‑123 (MQ=255)
                                                 ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACtt                                                                                                  >  8:208934/1‑149 (MQ=255)
                                                 ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGCTCATGAGAACCGGGGGCttttt                                                                                                                                 <  6:222/118‑1 (MQ=37)
                                                 ccAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGCTCATGAGAACCGGGGGCttttt                                                                                                                                 >  5:222/1‑118 (MQ=37)
                                                      aCATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCTAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCa                                                                                             <  5:357368/149‑1 (MQ=255)
                                                        aTCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCTGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  3:41420/149‑1 (MQ=255)
                                                        aTCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  6:296096/149‑1 (MQ=255)
                                                         tCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAag                                                                                          <  1:16084/149‑1 (MQ=255)
                                                           aCCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAagag                                                                                        >  8:22877/1‑149 (MQ=255)
                                                             cTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAG‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGc                                                                                                                           >  6:270055/1‑111 (MQ=255)
                                                             cTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAG‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGc                                                                                                                           <  5:270055/111‑1 (MQ=255)
                                                                ggTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTc                                                                                     <  5:320226/147‑1 (MQ=255)
                                                                  tCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGa                                                                                 <  3:325675/149‑1 (MQ=255)
                                                                    aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           <  2:215272/137‑1 (MQ=255)
                                                                    aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                           >  1:215272/1‑137 (MQ=255)
                                                                      cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGc                                                                             >  7:259118/1‑149 (MQ=255)
                                                                       ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                              >  7:275981/1‑115 (MQ=255)
                                                                       ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                              <  8:275981/115‑1 (MQ=255)
                                                                       ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGc                                                                             >  2:342006/1‑148 (MQ=255)
                                                                              gcgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagc                                                                     <  1:342006/149‑1 (MQ=255)
                                                                               cgcTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCcagca                                                                    <  7:218975/149‑1 (MQ=255)
                                                                                     cGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                              <  7:22877/149‑1 (MQ=255)
                                                                                        gCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAACCTTCTCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGa                                                              <  1:119276/149‑1 (MQ=255)
                                                                                         ccGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAg                                                          >  7:341128/1‑149 (MQ=255)
                                                                                             aCTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTc                                                      >  6:372056/1‑149 (MQ=255)
                                                                                                 tCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCAcc                                                   >  8:165165/1‑148 (MQ=255)
                                                                                                   cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAAc                                                >  2:296340/1‑149 (MQ=255)
                                                                                                   cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAAc                                                >  5:302169/1‑149 (MQ=255)
                                                                                                    gAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCAccccgc                                               >  3:51301/1‑145 (MQ=255)
                                                                                                             aGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGc                                       <  5:76163/148‑1 (MQ=255)
                                                                                                               gCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCggg                                    <  1:205863/149‑1 (MQ=255)
                                                                                                               gCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCCCGCCAGTGGCgttg                                   <  4:51301/149‑4 (MQ=255)
                                                                                                                ccccTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCgggg                                   >  7:248371/1‑149 (MQ=255)
                                                                                                                    tGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGggtggt                               >  4:196493/1‑149 (MQ=255)
                                                                                                                            cGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCg                        <  8:341128/148‑1 (MQ=255)
                                                                                                                             ggCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtt                      <  7:165165/149‑1 (MQ=255)
                                                                                                                              gCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGttt                     <  3:196493/149‑1 (MQ=255)
                                                                                                                               caacaaAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGtttt                    <  8:248371/149‑1 (MQ=255)
                                                                                                                                 acaaAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                              >  7:61158/1‑89 (MQ=255)
                                                                                                                                 acaaAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATgg                                                                              <  8:61158/89‑1 (MQ=255)
                                                                                                                                 acaaAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCa                  <  6:302169/149‑1 (MQ=255)
                                                                                                                                   aaaaaGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCAt                 >  6:189650/1‑148 (MQ=255)
                                                                                                                                      aaGC‑‑CCCCGGGTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACg             <  8:259118/149‑1 (MQ=255)
                                                                                                                                       aGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg            <  1:209606/149‑1 (MQ=255)
                                                                                                                                       aGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg            <  5:189650/149‑1 (MQ=255)
                                                                                                                                                   ttCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCa   >  3:86648/1‑149 (MQ=255)
                                                                                                                                                     ctcATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAt  <  5:372056/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
ACGCGCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAT  >  NC_002947/4749870‑4750163

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