Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,022,306 G→C 6.9% intergenic (‑1/+667) PP_3547 ← / ← emrB NAD(P)‑binding oxidoreductase/multidrug transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,022,3060GC 6.9% 159.6 / 3.8 58intergenic (‑1/+667)PP_3547/emrBNAD(P)‑binding oxidoreductase/multidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (2/2);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.49e-01

CGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGA  >  NC_002947/4022160‑4022453
                                                                                                                                                  |                                                                                                                                                    
cgcTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGa                                                                                                                                                    <  2:82309/149‑1 (MQ=255)
cgcTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGAGCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGa                                                                                                                                                    <  5:42931/149‑1 (MQ=255)
     cGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                                               >  6:81842/1‑149 (MQ=255)
     cGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGccc                                                                                                                                               >  2:64200/1‑149 (MQ=255)
     cGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGc                                                                                                                                                 <  2:212169/147‑1 (MQ=255)
       ttGGAATGGCGAACATCAGCACAGCCAACCGTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGATGACGGGCATCCATGCCCGt                                                                                                                                             >  1:38552/1‑149 (MQ=21)
         ggAATGGCGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGATGACGGGCATCCATGCCCGTc                                                                                                                                            <  2:38552/148‑1 (MQ=255)
           aaTGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                           >  4:21851/1‑147 (MQ=255)
           aaTGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATCCATGCCCGTCa                                                                                                                                           <  3:21851/147‑1 (MQ=255)
             tGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATGTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                                                       >  7:100589/1‑149 (MQ=255)
             tGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                                                       >  2:142866/1‑149 (MQ=255)
             tGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAgcg                                                                                                                                       <  1:224160/149‑1 (MQ=255)
               aCGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCt                                                                                                                                     <  5:67385/149‑1 (MQ=255)
                 gAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATGTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                   <  8:100589/149‑1 (MQ=255)
                  aaCATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGAAGCCCGTCAAGCGCTGAg                                                                                                                                  <  3:207769/149‑1 (MQ=255)
                          cacaGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAg                                                                                                                          >  1:431276/1‑148 (MQ=255)
                                 aaCCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGc                                                                                                                  >  2:345055/1‑149 (MQ=255)
                                      tctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                             >  4:275111/1‑149 (MQ=255)
                                      tctGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCg                                                                                                             >  4:81766/1‑149 (MQ=255)
                                          cttttgatctggcttttgatctggCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGgcgc                                                                                                         >  7:140094/1‑149 (MQ=255)
                                           ttttgatctggcttttgatctggcTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                   <  7:154951/123‑1 (MQ=255)
                                           ttttgatctggcttttgatctggcTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGa                                                                                                                                   >  8:154951/1‑123 (MQ=255)
                                             ttgatctggcttttgatctggcttTTGGTCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGcgcg                                                                                                        <  2:118198/147‑1 (MQ=255)
                                             ttgatctggcttttgatctggcttTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGt                                                                                                      >  5:34824/1‑149 (MQ=255)
                                               gatctggcttttgatctggcttttGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGGAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTc                                                                                                     <  1:142866/148‑1 (MQ=255)
                                               gatctggcttttgatctggcttttGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTc                                                                                                     <  6:202185/148‑1 (MQ=255)
                                                 tctggctattgatctggcatttgaTCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTg                                                                                                  >  6:339923/1‑149 (MQ=255)
                                                     gcttttgatctgggttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCgg                                                                                              <  1:357085/149‑1 (MQ=38)
                                                         ttgatctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGCGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGa                                                                                            <  3:392300/147‑1 (MQ=37)
                                                             tctggcttttggtctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAgg                                                                                      >  8:284417/1‑149 (MQ=25)
                                                             tctggcttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAgg                                                                                      >  2:96968/1‑149 (MQ=35)
                                                                  cttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCTAACGCTGGCCTTGCGGAGGCAGGTGACGGGCATGGAGGCCCGTCATGCGCTGAGGCCCCTGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGa                                                                                 <  8:140218/149‑1 (MQ=11)
                                                                  cttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGCGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCgg                                                                                  <  1:96968/149‑1 (MQ=17)
                                                                  cttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGa                                                                                 <  3:134848/149‑1 (MQ=32)
                                                                   ttttgatctTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCg                                                                                   >  4:134848/1‑146 (MQ=31)
                                                                    tttgatctTGATCTTGATCTTGATCTACGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTCACGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGt                                                                               >  1:70127/1‑149 (MQ=2)
                                                                           ttgatcttgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGaa                                                                        >  6:124902/1‑149 (MQ=31)
                                                                                  tgatcttgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCg                                                                   >  2:51486/1‑147 (MQ=32)
                                                                                        tgatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcag                                                           >  3:204406/1‑149 (MQ=32)
                                                                                         gatcttCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagc                                                               <  7:284417/144‑1 (MQ=32)
                                                                                               cgcgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGagcgcagcgc                                                        <  7:5933/145‑1 (MQ=32)
                                                                                                gcgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAgggg                                                   >  3:155611/1‑149 (MQ=32)
                                                                                                gcgACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAggg                                                    <  4:338394/148‑1 (MQ=32)
                                                                                                       cAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGAGGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGgat                                             <  5:81842/148‑1 (MQ=17)
                                                                                                                gAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGa                                      >  1:209399/1‑146 (MQ=32)
                                                                                                                 aaCGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAgcgc                                  <  3:275111/149‑1 (MQ=32)
                                                                                                                        gCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGttt                           <  1:345055/149‑1 (MQ=32)
                                                                                                                        gCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGttt                           >  6:377611/1‑149 (MQ=32)
                                                                                                                         gcTTGAGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCGGGGAGCAAGGCCGGATTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGtttt                          <  5:179774/148‑1 (MQ=12)
                                                                                                                          cTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGttttt                         >  6:23299/1‑149 (MQ=32)
                                                                                                                            tGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                       >  2:114884/1‑149 (MQ=32)
                                                                                                                             ggggAGGCAGGTTCCGGGACTGGAGGCCCGTCAAGCGCTGAGGCCCCCGGTTGGGGCCTGCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                       <  5:339923/146‑1 (MQ=11)
                                                                                                                             gCGGTGGCAGGTGACGGGCATGTATGCCCGTAAAGGGGTGAGGCCGCTGGTTGGGGCCTTCGGGGGGGTCCTGCCGGTAGCAAGGCGGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                       <  2:70127/148‑1 (MQ=11)
                                                                                                                             gCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTgg                      >  3:291588/1‑149 (MQ=32)
                                                                                                                             gCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTg                       <  5:23299/148‑1 (MQ=32)
                                                                                                                                     aGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACtttt              >  3:255845/1‑149 (MQ=32)
                                                                                                                                              gCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTTGCGtt     >  6:114132/1‑149 (MQ=17)
                                                                                                                                              gCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGtt     <  4:291588/149‑1 (MQ=32)
                                                                                                                                                  ggATGCCGGTCATGCGCTGAGGCCCCTGGAGGGGGCGTTCGGCGGGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGa  <  4:155611/148‑1 (MQ=11)
                                                                                                                                                  |                                                                                                                                                    
CGCTGCGTTGGAATGACGAACATCAGCACAGCCAACCCTCTGCTTTTGATCTGGCTTTTGATCTGGCTTTTGATCTTGATCTTGATCTTGATCTTCGCGACTTCAGGAGGCCGAACGCAGGCCTTGCGGAGGCAGGTGACGGGCATGGATGCCCGTCAAGCGCTGAGGCCCCAGGATGGGGCCTTCGGCGCGGTCCTGCCGGGAGCAAGGCCGGAGTGAGGGAACCCCGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCGCAGCGTTTTTTGGTTACTTTTTGTCGCGTTTGA  >  NC_002947/4022160‑4022453

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