Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,639,576 G→T 13.6% intergenic (‑166/‑145) czcR‑III ← / → PP_1439 response regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,639,5760GT13.6% 148.9 / 14.9 59intergenic (‑166/‑145)czcR‑III/PP_1439response regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base T (3/5);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.28e-01

GAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGA  >  NC_002947/1639428‑1639721
                                                                                                                                                    |                                                                                                                                                 
gAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCcgcg                                                                                                                                                   >  3:124334/1‑149 (MQ=255)
   aGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCt                                                                                                                                                >  7:220883/1‑149 (MQ=255)
    ggggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTg                                                                                                                                               <  4:144336/149‑1 (MQ=255)
     gggACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTg                                                                                                                                               >  1:237853/1‑148 (MQ=255)
       gACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAg                                                                                                                                            <  3:149251/149‑1 (MQ=255)
               ttGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGt                                                                                                                                    <  1:190051/149‑1 (MQ=255)
                 gACCCTGTAGCGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAACGCCTGCGGTGCAGCTAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACggg                                                                                                                                     <  8:143922/146‑1 (MQ=255)
                     cTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTActc                                                                                                                              <  2:259693/149‑3 (MQ=37)
                        tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcgcg                                                                                                                           <  4:215897/149‑1 (MQ=37)
                        tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                            <  1:335722/148‑1 (MQ=255)
                        tAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                            <  1:340243/148‑1 (MQ=255)
                         gatctcATCAGGGTCATGCGCGCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcgcg                                                                                                                           <  6:62664/143‑1 (MQ=25)
                         aGTGACATCAGGGTCAAGCGCCACAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTcgcg                                                                                                                           <  8:125168/148‑1 (MQ=25)
                          gTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAc                                                                                                                                        <  1:154816/134‑1 (MQ=255)
                          gTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAc                                                                                                                                        >  2:160381/1‑134 (MQ=255)
                          gTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAc                                                                                                                                        >  2:154816/1‑134 (MQ=255)
                          gTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAAc                                                                                                                                        <  1:160381/134‑1 (MQ=255)
                               aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                             >  2:195606/1‑124 (MQ=255)
                               aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTa                                                                                                                                             <  1:195606/124‑1 (MQ=255)
                               aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                            <  8:190796/141‑1 (MQ=255)
                               aTCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                            >  7:190796/1‑141 (MQ=255)
                               aTCAGGGTCAAGCGCCCCAATGTTAAGGCATGAAAAAAGGCCTGCGGAGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGGGTAAACCACTTCATACAGGCGCGGGTTTACTCGCGAagatac                                                                                                                    >  5:62664/1‑146 (MQ=25)
                                  aGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGcac                                                                                                                  >  3:299861/1‑146 (MQ=255)
                                     gTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                            <  5:296916/135‑1 (MQ=255)
                                     gTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCcgc                                                                                                                            >  6:296916/1‑135 (MQ=255)
                                         aGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTCCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                          <  7:214883/149‑1 (MQ=255)
                                          gcgcCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCGTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                          >  3:342518/1‑148 (MQ=255)
                                          gcgcCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTc                                                                                                          >  1:197375/1‑148 (MQ=255)
                                             caccAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACTGATTTGTAATTCTGCTGTTATCTGGCTGGCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAATAGGCCCGAACAGTCcaca                                                                                                      <  6:170553/147‑1 (MQ=255)
                                                cAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCcacac                                                                                                     >  6:129559/1‑147 (MQ=255)
                                                     ttACGGCTTGAAAAACTGCCTGCGGTGCATATAAGGTTTTTTTATTTATGGTGTTATGTCGCTGTCTGTTGGGGCCTGTTCGCGAGTAAACCCGGGCGTGTAGGTACGGGTTAACCGGCTATGAGGCCCGAACAGTCCACACAAGGcc                                                                                               <  6:337894/148‑1 (MQ=255)
                                                      tACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTGca                                                                                                        >  3:154307/1‑138 (MQ=255)
                                                        cGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACAt                                                                                           <  5:129559/149‑1 (MQ=255)
                                                                aaaaaGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGAc                                                                                   <  1:90661/149‑1 (MQ=255)
                                                                      gCCTGAGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACACGCGCCTGTAGGAACGGGTTTACACGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAAt                                                                             >  5:281632/1‑149 (MQ=255)
                                                                                cAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCAGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGCCTGTAATCCCCGTCTCa                                                                   >  2:98789/1‑149 (MQ=255)
                                                                                 aGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCAc                                                                  >  6:97726/1‑149 (MQ=255)
                                                                                       cGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCCCCTTGcc                                                            >  3:18980/1‑149 (MQ=255)
                                                                                       cGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACACGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGcc                                                            >  7:328351/1‑149 (MQ=255)
                                                                                         gATTTCTAATTCTGCTGTTATCTGGCGGTCTGTTCGGGCCGCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCCAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGt                                                          <  4:124334/149‑1 (MQ=255)
                                                                                              gTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCca                                                                                                        >  8:350772/1‑98 (MQ=255)
                                                                                                  ttCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTagccagcc                                                 >  4:325714/1‑149 (MQ=255)
                                                                                                       ctgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCg                                             <  5:97726/148‑1 (MQ=255)
                                                                                                        tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTTCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCATCGcg                                           >  8:344740/1‑148 (MQ=255)
                                                                                                        tgTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGc                                            >  2:50485/1‑148 (MQ=255)
                                                                                                                     tctgtTCGGGCCTCTTCGCGCGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGCAGAGGCCCGAACAGTGCCCCCACGGCCACATTGCAGGACTGTAATCCCCGCCTGCCCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                              >  7:33534/1‑149 (MQ=255)
                                                                                                                       tgtTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATCGCCGTGTCt                            >  1:316935/1‑149 (MQ=255)
                                                                                                                         ttCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGa                          >  7:26629/1‑149 (MQ=255)
                                                                                                                               ccTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                              >  1:108774/1‑139 (MQ=255)
                                                                                                                               ccTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCgtgt                              <  2:108774/139‑1 (MQ=255)
                                                                                                                                   ttCGCGGGTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGttatta                >  5:292710/1‑149 (MQ=25)
                                                                                                                                     cgcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTAc               <  3:325714/148‑1 (MQ=255)
                                                                                                                                      gcgGGTAAACCAGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa             <  6:292710/149‑1 (MQ=37)
                                                                                                                                      gcgAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGa             <  2:197375/149‑1 (MQ=255)
                                                                                                                                        gggTAAACCCGTTCCTACAGGCGCGGGTTTACTCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAgg           >  4:168682/3‑149 (MQ=37)
                                                                                                                                         aGTAAACCCGCGCATGTAGGAACGGGTTTACCCGCGAAGAGGCCCTAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc          <  8:328351/149‑1 (MQ=255)
                                                                                                                                               cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc          >  5:265902/1‑143 (MQ=255)
                                                                                                                                               cccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGc          <  6:265902/143‑1 (MQ=255)
                                                                                                                                                ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac    >  3:251956/1‑148 (MQ=255)
                                                                                                                                                ccGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAacac    >  4:42680/1‑148 (MQ=255)
                                                                                                                                                 cgcgCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGAGCGGTTATTACGAGGCAAACACGa  <  6:281632/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GAAAGGGGACGATTCTTGACCCTGTAGTGACATCAGGGTCAAGCGCCCCAATGTTACGGCATGAAAAAAGGCCTGCGGTGCAGATAACGGATTTGTAATTCTGCTGTTATCTGGCTGTCTGTTCGGGCCTCTTCGCGAGTAAACCCGCGCCTGTAGGAACGGGTTTACCCGCGAAGAGGCCCGAACAGTCCACACAAGGCCACATTGCAGGACTGTAATCCCCGCCTCACCTTGCCGTTAGCCAGCCCTCGCTAGGATGGCCGTGTCTGATCGGTTATTACGAGGCAAACACGA  >  NC_002947/1639428‑1639721

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