Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,906,189 T→C 6.2% F815S (TTT→TCT)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,906,1890TC 6.2% 184.9 / 2.0 64F815S (TTT→TCT) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base C (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

GACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCAATGCTGCTGTGAACACGCAA  >  NC_002947/1906044‑1906330
                                                                                                                                                 |                                                                                                                                             
gACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGt                                                                                                                                            >  2:20900/1‑149 (MQ=255)
gACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGt                                                                                                                                            >  6:87533/1‑149 (MQ=255)
      cTGAAGCAGGTGATAATCGACGACAGCCTGGGCCTGTCCAGCGAGCTGGAAAGACAGAAGCAGCATGTGGTCGACCAGTATGAGTGCGAATGGGAAAATGCCCGCAACGACCCGGAGAAGATCAAGCTCTGCCGTACCTTTGTCAACg                                                                                                                                       >  6:146103/1‑148 (MQ=255)
            cAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGCCCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATcgccg                                                                                                                                >  6:265263/1‑149 (MQ=255)
             aGGTGATCATCGACGACAGCCTGGGCATGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATcgccgc                                                                                                                               <  1:20900/149‑1 (MQ=255)
                       cgacgaCAGCCTGGGCCTGGCCAGCGAGCTAGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCg                                                                                                                      >  1:198238/1‑148 (MQ=255)
                        gacgacAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGa                                                                                                                    >  2:179186/1‑149 (MQ=255)
                         acgacAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGAc                                                                                                                   >  5:183443/1‑149 (MQ=255)
                             cAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACTATCGCCGCGCCGACCCGGACGTGc                                                                                                               <  6:183443/149‑1 (MQ=255)
                              aGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGc                                                                                                               >  5:63048/1‑148 (MQ=255)
                              aGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCa                                                                                                              >  2:13451/1‑149 (MQ=255)
                                cctgggcctggCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACt                                                                                                            >  7:218761/1‑149 (MQ=255)
                                     gcctggCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTg                                                                                                       >  1:40580/1‑149 (MQ=255)
                                      cctggCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTgc                                                                                                      <  5:311162/149‑1 (MQ=255)
                                      cctggCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTgc                                                                                                      <  6:255447/149‑1 (MQ=255)
                                          gCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTgcgc                                                                                                    <  8:242057/147‑1 (MQ=255)
                                                 aGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGcgcgaacgcgaa                                                                                           <  6:301673/149‑1 (MQ=255)
                                                  gCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCgcgaacgcgaac                                                                                          <  5:87533/149‑1 (MQ=255)
                                                      gAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACg                                                                                      >  2:46406/1‑149 (MQ=255)
                                                           ccAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAGGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGc                                                                                 >  8:175329/1‑149 (MQ=255)
                                                           ccAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGCACGTGAACTTCGAGCGCGAACGCCAACAACGCCGGc                                                                                 >  6:202986/1‑149 (MQ=255)
                                                                   agcaTGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCCACCCGGAAGTGCACTTCCTGCGCGAACGCGAACAACGACGGCccgctgca                                                                         >  7:250209/1‑149 (MQ=255)
                                                                        gtgGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCccgctgcaccgct                                                                    >  2:269453/1‑149 (MQ=255)
                                                                           gTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCccg                                                                              <  6:133681/136‑1 (MQ=255)
                                                                           gTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCccg                                                                              >  5:133681/1‑136 (MQ=255)
                                                                            tCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCcgctgcaccgctgcac                                                                <  7:238180/149‑1 (MQ=255)
                                                                            tCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCcgctgcaccgctgcac                                                                >  8:238180/1‑149 (MQ=255)
                                                                              gtccATTATGAGTGCGAATGGGCCATTGCCCTCAACGACCCGGAGAAGCTCTAGCGCTTCCGTACGTTTGTCATCGATCGCCGCGCCGTCCCGGTCGTGCACTTCGTGCGCGATCGCGAACTACGCCGGCCCgctgcaccgctgcacc                                                               <  5:146103/146‑1 (MQ=255)
                                                                              gACCAGTATGAATGGGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGTACAACGCCGGCCCGCTGCACCGCTGCACCt                                                              <  3:282089/149‑1 (MQ=255)
                                                                                 cAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGAt                                                           <  8:218761/149‑1 (MQ=255)
                                                                                     aTGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTTCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATccc                                                        <  6:63048/148‑1 (MQ=255)
                                                                                      tGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAGGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCaa                                                      <  7:175329/149‑1 (MQ=255)
                                                                                             gAATGGGCCAATGCCCTCAACGTCCCGGAGACGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATACCAACTGTCga                                               <  1:135686/149‑1 (MQ=255)
                                                                                             gAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCaa                                                      >  4:163246/1‑142 (MQ=255)
                                                                                             gAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCaa                                                      <  3:163246/142‑1 (MQ=255)
                                                                                                  ggCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCATCGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGaggagg                                          <  1:269453/149‑1 (MQ=255)
                                                                                                   gCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCAAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGc                                         >  2:61151/1‑149 (MQ=255)
                                                                                                    ccAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATccc                                                        >  8:223067/1‑133 (MQ=255)
                                                                                                    ccAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATccc                                                        <  7:223067/133‑1 (MQ=255)
                                                                                                      aaTGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGt                                      <  1:13451/149‑1 (MQ=255)
                                                                                                      aaTGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCAGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCCGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCAGt                                      >  5:158620/1‑149 (MQ=255)
                                                                                                         gCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCa                                                                                                              <  1:87495/74‑1 (MQ=255)
                                                                                                         gCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCa                                                                                                              >  2:87495/1‑74 (MQ=255)
                                                                                                          cccTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTtg                                   <  1:179186/148‑1 (MQ=255)
                                                                                                           ccTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTCTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCAGCACCGCTGCACCTGATACCAACTGTCGAAGAGCCTGTTtgct                                 >  6:38435/1‑147 (MQ=255)
                                                                                                           ccTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTCTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAAAGCGAACAAAGCCGGCCCGCTGCAACGCTGCACATGATCCAAACTGTCGATGAGGCTGTTtgat                                 >  6:33509/1‑149 (MQ=255)
                                                                                                            cTCAACGACCCGGAGATGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACGCGGACGTGCACTTCGTGCGCGAACGCGAACAACGGCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTtgatg                                <  5:265263/149‑1 (MQ=255)
                                                                                                            cTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTc                                                 <  1:246871/132‑1 (MQ=255)
                                                                                                            cTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTc                                                 >  2:246871/1‑132 (MQ=255)
                                                                                                             tCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTtgatg                                >  2:97614/1‑148 (MQ=255)
                                                                                                               aaCGACCCGGATAAGCTCAAGCGCTTCCGTACCTCTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAAc                             <  5:33509/149‑1 (MQ=255)
                                                                                                               aaCGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAAc                             <  1:61151/149‑1 (MQ=255)
                                                                                                               aaCGACCCGGAGAAGCTCAAGCGCTTCCGTACCTCTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAAc                             <  5:38435/149‑1 (MQ=255)
                                                                                                                 cGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTtg                                   <  8:210815/141‑1 (MQ=255)
                                                                                                                 cGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTtg                                   >  7:210815/1‑141 (MQ=255)
                                                                                                                    cccGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGc                        >  3:107860/1‑149 (MQ=255)
                                                                                                                       ggAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCaa                     >  2:53734/1‑149 (MQ=255)
                                                                                                                         agaAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCTCGAACGCGAACAACTCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCAAtg                   <  5:202986/149‑1 (MQ=255)
                                                                                                                             gCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCCCCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCAAtgctgc               >  4:174825/1‑149 (MQ=255)
                                                                                                                                 aaGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACCCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTTAAGAACCTGGCCCATGCTGCtgt            >  6:207312/1‑148 (MQ=255)
                                                                                                                                     gcTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGt                                      <  4:59448/118‑1 (MQ=255)
                                                                                                                                     gcTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGt                                      >  3:59448/1‑118 (MQ=255)
                                                                                                                                          cGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCAATGCTGCTGTGAACACGCaa  <  1:46406/149‑1 (MQ=255)
                                                                                                                                          cGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCAATGCTGCTGTGAACACGCaa  <  5:207312/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
GACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGAATGCGAATGGGCCAATGCCCTCAACGACCCGGAGAAGCTCAAGCGCTTCCGTACCTTTGTCAACGATCGCCGCGCCGACCCGGACGTGCACTTCGTGCGCGAACGCGAACAACGCCGGCCCGCTGCACCGCTGCACCTGATCCCAACTGTCGAGGAGGCTGTTTGATGAACCTGGCCAATGCTGCTGTGAACACGCAA  >  NC_002947/1906044‑1906330

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: