Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,750,020 T→A 6.3% 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,750,0200TA 6.3% 218.4 / 2.7 74intergenic (+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 (32/37);  new base A (3/2);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAC  >  NC_002947/4749874‑4750167
                                                                                                                                                     |                                                                                                                                                   
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                                                                                                                                                   <  1:84243/128‑1 (MQ=255)
                     gACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGa                                                                                                                                                   >  2:84243/1‑128 (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                                                                                               >  1:215272/1‑137 (MQ=255)
                                                                aCCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCaaa                                                                                               <  2:215272/137‑1 (MQ=255)
                                                                  cccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGc                                                                                 >  7:259118/1‑149 (MQ=255)
                                                                   ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                  <  8:275981/115‑1 (MQ=255)
                                                                   ccTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAgg                                                                                                                  >  7:275981/1‑115 (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                                                    >  5:302169/1‑149 (MQ=255)
                                                                                               cGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAAc                                                    >  2:296340/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                <  5:189650/149‑1 (MQ=255)
                                                                                                                                   aGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACgg                <  1:209606/149‑1 (MQ=255)
                                                                                                                                               ttCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCa       >  3:86648/1‑149 (MQ=255)
                                                                                                                                                 ctcATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCAt      <  5:372056/148‑1 (MQ=255)
                                                                                                                                                    aTGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAc  >  5:24151/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCAGAACTGCGTGCACTGCAAGACCTGCGACATCAAGGACCCGGCCCAGAACATCACCTGGGTCACCCCTGAAGGCGCTGGCGGGCCGAACTATCCGAACATGTAAGGCCCCTGCCCTTGCGGCAACAAAAAGC‑‑CCCCGGTTCTCATGAGCCGGGGGCTTTTTGTTGGCTGCCATTATCAGGCAGTACGCACTTCCTCAAAGAGTTCTGGATGGCGATCCAGCAGCTTGAGCAGTTTCACCAACGCCACTGGCGGGGTGGTCTTGCCGTTTTCATAACGGGAGAACGCATTGAC  >  NC_002947/4749874‑4750167

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: