Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,766 G→A 6.1% R674H (CGC→CAC)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,7660GA 6.1% 197.8 / 2.3 65R674H (CGC→CAC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base A (2/2);  total (34/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.98e-01

GCCGGTTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCC  >  NC_002947/1905627‑1905911
                                                                                                                                           |                                                                                                                                                 
gCCGGTTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGccccc                                                                                                                                           <  8:19240/148‑1 (MQ=255)
  cGGTTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCcac                                                                                                                                         >  7:10155/1‑148 (MQ=255)
     tttCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCACAGCCCCCACAAGc                                                                                                                                     >  8:241541/1‑149 (MQ=255)
        cGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTc                                                                                                                                   >  5:261121/1‑148 (MQ=255)
         gAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTc                                                                                                                                   >  7:260954/1‑147 (MQ=255)
              aggCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGTGCAGCACCGGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTcg                                                                                                                            >  2:220700/2‑149 (MQ=255)
                gCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTTGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAATATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTcgcg                                                                                                                          <  3:138034/149‑1 (MQ=255)
                gCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTcgcg                                                                                                                          <  6:39649/149‑1 (MQ=255)
                          cGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCggg                                                                                                                  <  4:38260/147‑1 (MQ=255)
                          cGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGAt                                                                                                                >  2:289211/1‑149 (MQ=255)
                                gTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGc                                                                                                                                     >  2:203799/1‑122 (MQ=255)
                                gTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGc                                                                                                                                     <  1:203799/122‑1 (MQ=255)
                                gTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCACAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCAccc                                                                                                          >  7:57456/1‑149 (MQ=255)
                                   gACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGt                                                                                                        <  1:10562/148‑1 (MQ=255)
                                      ccGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAAt                                                                                                    <  7:30913/149‑1 (MQ=255)
                                           cGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCAccc                                                                                                          >  4:22660/1‑138 (MQ=255)
                                           cGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCAccc                                                                                                          <  3:22660/138‑1 (MQ=255)
                                            ggTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCg                                                                                              >  1:154856/1‑149 (MQ=255)
                                                 aGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGgcg                                                                                         >  6:77353/1‑149 (MQ=255)
                                                 aGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGgcg                                                                                         >  2:54474/1‑149 (MQ=255)
                                                  gTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGt                                                                                                        <  5:171035/133‑1 (MQ=255)
                                                  gTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGt                                                                                                        >  6:171035/1‑133 (MQ=255)
                                                  gTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGgcg                                                                                         >  7:79774/1‑148 (MQ=255)
                                                   tCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGgcg                                                                                         >  5:270572/1‑147 (MQ=255)
                                                          tgGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAg                                                                                                                                                <  5:330667/85‑1 (MQ=255)
                                                          tgGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAg                                                                                                                                                >  6:330667/1‑85 (MQ=255)
                                                          tgGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAg                                                                                <  8:10155/149‑1 (MQ=255)
                                                          tgGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCACAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAg                                                                                <  8:57456/149‑1 (MQ=255)
                                                             gcagcaCCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGAc                                                                             <  1:54474/149‑1 (MQ=255)
                                                                gcaCCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCa                                                                                       >  3:230587/1‑136 (MQ=255)
                                                                gcaCCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCa                                                                                       <  4:230587/136‑1 (MQ=255)
                                                                   ccTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGc                                                                       >  8:233308/1‑149 (MQ=255)
                                                                     tGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGt                                                                     >  1:109180/1‑149 (MQ=255)
                                                                        tGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTAt                                                                  <  7:233308/149‑1 (MQ=255)
                                                                                       cAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAggg                                                   >  8:2186/1‑149 (MQ=255)
                                                                                       cAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAggg                                                   >  2:197445/1‑149 (MQ=255)
                                                                                           aCAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTgg                                               <  8:260954/149‑1 (MQ=255)
                                                                                            cAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTTAATGCGCCGAGGCGGAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGa                                              <  1:220700/149‑1 (MQ=255)
                                                                                             aGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAggg                                                   >  5:210596/1‑143 (MQ=255)
                                                                                             aGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAggg                                                   <  6:210596/143‑1 (MQ=255)
                                                                                                 tCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTg                                         <  1:289211/149‑1 (MQ=255)
                                                                                                   gCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAgg                                                                               <  1:115310/109‑1 (MQ=255)
                                                                                                   gCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAgg                                                                               >  2:115310/1‑109 (MQ=255)
                                                                                                   gCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGt                                        <  6:270572/148‑1 (MQ=255)
                                                                                                   gCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTa                                       >  7:206359/1‑149 (MQ=255)
                                                                                                       tGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAAc                                            >  8:299901/1‑140 (MQ=255)
                                                                                                       tGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAAc                                            <  7:299901/140‑1 (MQ=255)
                                                                                                         gCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCACAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGtgtg                                 <  7:241541/149‑1 (MQ=255)
                                                                                                          cccTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAg                                                                                <  1:41958/101‑1 (MQ=255)
                                                                                                          cccTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAg                                                                                >  2:41958/1‑101 (MQ=255)
                                                                                                           ccTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCg                               <  8:206359/149‑1 (MQ=255)
                                                                                                                  cTGGAAGATCGCTACAAGGGCCTGCGCACCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAAcgg                        >  1:183692/1‑149 (MQ=255)
                                                                                                                     gAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCaga                                                                                     <  4:159957/85‑1 (MQ=255)
                                                                                                                     gAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCaga                                                                                     >  3:159957/1‑85 (MQ=255)
                                                                                                                        gATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCAt                  <  2:12886/149‑1 (MQ=255)
                                                                                                                        gATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCAt                  >  2:191183/1‑149 (MQ=255)
                                                                                                                          tCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATgc                >  7:69586/1‑149 (MQ=255)
                                                                                                                            gCTACAAGGGCCTGCGCAGCACCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCAGTGAATGCGCAGAGGCGCAGAGCAAGGCAATAGGCGTTATAGCCACCGACAAGGGCTGGAACCTGTACGTGTGCTGTAACTGCGGCATgcgc              >  5:200975/1‑149 (MQ=255)
                                                                                                                             cTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGcc             >  3:236208/1‑149 (MQ=255)
                                                                                                                             cTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGACTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGcc             >  5:223639/1‑149 (MQ=255)
                                                                                                                                cAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCcgc          <  8:69586/149‑1 (MQ=255)
                                                                                                                                cAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCcgc          <  8:79774/149‑1 (MQ=255)
                                                                                                                                  aGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCgg        >  3:61933/1‑149 (MQ=255)
                                                                                                                                       cTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGc   <  4:61933/149‑1 (MQ=255)
                                                                                                                                       cTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGc   <  2:109180/149‑1 (MQ=255)
                                                                                                                                          cgcAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGcc  <  1:197445/147‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                 
GCCGGTTTCGAAACCGGCCACGCCTACGGCAAGTCGACCCGCACGGTAAAGTCGTGCGTGGGCAGCACCTGGTGCCGTTACGGCGTGCAGGACAGTGTCGCCATGGCCCTGCGCCTGGAAGATCGCTACAAGGGCCTGCGCAGCCCCCACAAGCTCAAGTTCGCGGTATCGGGATGCACCCGTGAATGCGCCGAGGCGCAGAGCAAGGACATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCC  >  NC_002947/1905627‑1905911

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