Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,620,348 C→G 7.8% intergenic (‑114/‑87) opdH ← / → tctD tricarboxylate‑specific outer membrane porin/transcriptional regulator TctD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,620,3480CG 7.8% 197.3 / 6.6 64intergenic (‑114/‑87)opdH/tctDtricarboxylate‑specific outer membrane porin/transcriptional regulator TctD
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base G (2/3);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

GCGCGAACGCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTGGCCCTGGC  >  NC_002947/1620204‑1620490
                                                                                                                                                |                                                                                                                                              
gcgcGAACGCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTg                                                                                                                                             >  8:220706/1‑148 (MQ=255)
 cgcgAACGCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGaa                                                                                                                                           <  1:68362/149‑1 (MQ=255)
 cgcgAACGCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGa                                                                                                                                            >  8:128367/1‑148 (MQ=255)
      aCGCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGtt                                                                                                                                      <  8:165706/149‑1 (MQ=255)
        gCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCACGCTGAAAGGttag                                                                                                                                    <  5:35729/149‑3 (MQ=255)
        gCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCACGCTGAAAGGttag                                                                                                                                    <  1:219543/149‑3 (MQ=255)
         ccTGCGGTTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCACGCTGAAAGGttagc                                                                                                                                   <  2:276707/149‑4 (MQ=255)
          cTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGt                                                                                                                                   >  4:38927/1‑148 (MQ=255)
               gCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAg                                                                                                                             >  2:110842/1‑149 (MQ=255)
                  gCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAgggg                                                                                                                          >  2:153934/1‑149 (MQ=255)
                          aGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAg                                                                                                                  >  8:44966/1‑149 (MQ=255)
                            cATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCATGAAAGGTTTGTCACGAGGGGCTTCACAGGc                                                                                                                >  6:67532/1‑149 (MQ=255)
                              tcgtcgATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCgt                                                                                                              <  3:62825/149‑1 (MQ=255)
                              tcgtcgATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCgt                                                                                                              <  1:64311/149‑1 (MQ=255)
                                gtcgATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCgtgg                                                                                                            <  1:110842/149‑1 (MQ=255)
                                gtcgATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCgtg                                                                                                             <  5:142528/148‑1 (MQ=255)
                                 tcgATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCgtggt                                                                                                           >  5:255700/1‑149 (MQ=255)
                                 tcgATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCgtggt                                                                                                           <  4:134342/149‑1 (MQ=255)
                                        cTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCg                                                                                                    >  7:93505/1‑149 (MQ=255)
                                           cGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGATCCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTa                                                                                                 <  7:220706/149‑1 (MQ=255)
                                           cGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGt                                                                                                  <  7:44966/148‑1 (MQ=255)
                                             tttATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTaca                                                                                               <  2:52613/149‑1 (MQ=255)
                                                aTTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTc                                                                                            >  1:47877/1‑149 (MQ=255)
                                                    ttCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCcgcg                                                                                        >  4:28457/1‑149 (MQ=255)
                                                         aTTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCAcc                                                                                   >  5:13010/1‑149 (MQ=255)
                                                          ttCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCg                                                                                  >  3:310881/1‑149 (MQ=255)
                                                          ttCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCg                                                                                  <  6:255700/149‑1 (MQ=255)
                                                             gTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGccc                                                                               >  1:24028/1‑149 (MQ=255)
                                                                gCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCa                                                                             <  7:128367/148‑1 (MQ=255)
                                                                 cACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCAc                                                                            <  2:47877/148‑1 (MQ=255)
                                                                  aCCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCAccc                                                                          >  7:28555/1‑149 (MQ=255)
                                                                    cTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCAccc                                                                          >  1:262134/1‑147 (MQ=255)
                                                                    cTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCAccc                                                                          >  7:18182/1‑147 (MQ=255)
                                                                    cTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCAccc                                                                          >  7:2775/1‑147 (MQ=255)
                                                                           gTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAggg                                                                  <  3:38927/148‑1 (MQ=255)
                                                                           gTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGga                                                                 <  3:28457/149‑1 (MQ=255)
                                                                             gCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGgag                                                                >  7:176941/1‑148 (MQ=255)
                                                                               ttttttCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAAt                                                             >  6:31644/1‑149 (MQ=255)
                                                                                 ttttCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATc                                                            <  5:48560/148‑1 (MQ=255)
                                                                                    tCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGAt                                                        >  6:281949/1‑149 (MQ=255)
                                                                                     cGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGAt                                                        <  5:281949/148‑1 (MQ=255)
                                                                                      ggCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGc                                                      <  8:75977/149‑1 (MQ=255)
                                                                                        cgcgAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGTGAGAATCCGATGCgt                                                    >  6:229751/1‑149 (MQ=255)
                                                                                         gcgAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCgtg                                                   >  2:113946/1‑149 (MQ=255)
                                                                                            aaGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGc                                                >  6:70107/1‑149 (MQ=255)
                                                                                                  gACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTgg                                          <  1:113946/149‑1 (MQ=255)
                                                                                                       tggtgCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGTCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGaa                                     <  8:93505/149‑1 (MQ=255)
                                                                                                         gtgCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACTAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAg                                    <  6:301093/148‑1 (MQ=255)
                                                                                                         gtgCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAg                                    <  6:13010/148‑1 (MQ=255)
                                                                                                         gtgCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCAGAGCAGGGAGAATCAGATGCGTGTGCTACTGGTCGAAGa                                   >  5:65272/1‑149 (MQ=255)
                                                                                                         gtgCGATGCTAGGCAAAGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAg                                    <  6:305770/148‑1 (MQ=255)
                                                                                                            cGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAATGGTTTGTCACGAGGGGCTTCACAGGCGGGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGacc                                 <  5:309624/148‑1 (MQ=255)
                                                                                                             gATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGaccac                               <  8:190678/149‑1 (MQ=255)
                                                                                                                gCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTg                            >  8:1252/1‑149 (MQ=255)
                                                                                                                        aCGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGc                    <  8:18182/149‑1 (MQ=255)
                                                                                                                              cTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGc                                                <  8:75561/115‑1 (MQ=255)
                                                                                                                              cTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAg              >  4:232706/1‑149 (MQ=255)
                                                                                                                               tttCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGc             >  2:187822/1‑149 (MQ=255)
                                                                                                                                  gacaaaCCTTTCACGCTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTg          >  3:180268/5‑149 (MQ=255)
                                                                                                                                   gCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTgg         <  2:24028/149‑1 (MQ=255)
                                                                                                                                   gCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTgg         >  1:203925/1‑149 (MQ=255)
                                                                                                                                     aaaCCTTTCACGCTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTGGcc       >  5:214642/2‑149 (MQ=255)
                                                                                                                                          tttCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCAGCTGCAACTGGCCGAAAGCGTGGCCCTGGc  >  8:258859/1‑149 (MQ=255)
                                                                                                                                            tCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTGGCCCTGGc  <  8:176941/147‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
GCGCGAACGCCTGCGGCTGCGATGACAGCATCGTCGATGACTCCGTTTATTGTTCTTATTCGTTGCACCTTGCTGGTGCTTTTTTCGGCGCGAAGGGCGACCGTGGTGCGATGCTAGGCAACGAACCTTTCGCTAACCTTTCAGCGTGAAAGGTTTGTCACGAGGGGCTTCACAGGCGTGGTGACAGCGGTACACTCCGCGCCACCGCCCCACCCGAGCAGGGAGAATCCGATGCGTGTGCTACTGGTCGAAGACCACCTGCAACTGGCCGAAAGCGTGGCCCTGGC  >  NC_002947/1620204‑1620490

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: ‑