Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,905,984 C→A 5.2% R747S (CGC→AGC)  nirB → nitrite reductase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,905,9840CA 5.2% 247.1 / 2.4 77R747S (CGC→AGC) nirBnitrite reductase large subunit
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base A (2/2);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.03e-01

CATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGA  >  NC_002947/1905836‑1906132
                                                                                                                                                    |                                                                                                                                                     
cATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATcc                                                                                                                                                       <  7:327042/149‑1 (MQ=255)
    ggCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGAGCCCGCGGCATGCCGAGCTGTTCGCCACCTCCATGGACGACGAATCCATGGTACGCCTGATCGACAGCGTGCTGATGTTCTACATCCGCAc                                                                                                                                                   >  8:123544/1‑149 (MQ=255)
          aTCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGa                                                                                                                                             >  8:276718/1‑149 (MQ=255)
             gCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACa                                                                                                                                           <  3:360024/148‑1 (MQ=255)
              ccACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAg                                                                                                                                         >  4:103638/1‑149 (MQ=255)
              ccACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCGGATGTAGAACATCAGCACgcggt                                                                                                                                              <  3:104836/144‑5 (MQ=255)
              ccACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCGGATGTAGAACATCAGCACgcggt                                                                                                                                              >  4:104836/1‑140 (MQ=255)
               cACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGc                                                                                                                                        >  7:38710/1‑149 (MQ=255)
                aCCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCt                                                                                                                                       >  3:248422/1‑149 (MQ=255)
                 ccGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTg                                                                                                                                      <  7:276718/149‑1 (MQ=255)
                    aCAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAg                                                                                                                                   >  5:388221/1‑149 (MQ=255)
                       aGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAgcgc                                                                                                                                <  1:307058/149‑1 (MQ=255)
                         ggCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTTGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCAc                                                                                                                              >  4:12606/1‑149 (MQ=255)
                         ggCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCAc                                                                                                                              >  2:425849/1‑149 (MQ=255)
                         ggCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCAc                                                                                                                              >  7:282209/1‑149 (MQ=255)
                                 ccTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGgtgt                                                                                                                      <  5:222930/149‑1 (MQ=255)
                                  cTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGtgtg                                                                                                                     >  5:61968/1‑149 (MQ=255)
                                  cTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGtgtg                                                                                                                     <  3:103638/149‑1 (MQ=255)
                                     tACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCATGGCCGACAAGCTGCAGCGCACCTCGGTGTgg                                                                                                                    <  1:265294/147‑1 (MQ=255)
                                     tACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGGCGACAAGCTGCAGCGCACCTCGGTGTgg                                                                                                                    <  2:338389/147‑1 (MQ=255)
                                     tACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGgc                                                                                                                   <  5:382988/148‑1 (MQ=255)
                                       cGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTTGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGcgcg                                                                                                                <  3:12606/149‑1 (MQ=255)
                                          gtgCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGaa                                                                                                              >  6:241780/1‑148 (MQ=255)
                                             cGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAAc                                                                                                           <  6:61968/148‑1 (MQ=255)
                                               gTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGgcg                                                                                                                  >  4:352885/1‑139 (MQ=255)
                                               gTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGgcg                                                                                                                  <  3:352885/139‑1 (MQ=255)
                                                tAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTg                                                                                                        >  5:124753/1‑148 (MQ=255)
                                                tAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTg                                                                                                        >  8:15786/1‑148 (MQ=255)
                                                          aTGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCt                                                                                             >  3:342143/1‑149 (MQ=255)
                                                           tGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTg                                                                                            <  1:244808/149‑1 (MQ=255)
                                                             cgcCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCt                                                                                                         >  7:359461/1‑134 (MQ=255)
                                                             cgcCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCt                                                                                                         <  8:359461/134‑1 (MQ=255)
                                                                 cgcgGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTAc                                                                                      <  8:38710/149‑1 (MQ=255)
                                                                  ggcgCATGCCGCGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTAcc                                                                                     <  1:15330/146‑1 (MQ=255)
                                                                  gcgGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTAcc                                                                                     <  8:98262/149‑1 (MQ=255)
                                                                          ccGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAg                                                                             >  5:255314/1‑149 (MQ=255)
                                                                                  ttCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGatc                                                                       <  5:122349/147‑1 (MQ=255)
                                                                                    cGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGatcatc                                                                    >  2:251178/1‑148 (MQ=255)
                                                                                     gCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcga                                                                  >  6:16548/1‑149 (MQ=255)
                                                                                      ccACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgac                                                                 <  4:342143/149‑1 (MQ=255)
                                                                                       cACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgac                                                                 >  2:367635/1‑148 (MQ=255)
                                                                                       cACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgac                                                                 >  6:11610/1‑148 (MQ=255)
                                                                                       cACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACATGGAAGGCGGGATGCACTACCTGAAGCAGGTGATCAGcgcc                                                                 >  4:307328/1‑146 (MQ=255)
                                                                                                ggACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgac                                                                 >  1:193781/1‑139 (MQ=255)
                                                                                                ggACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgac                                                                 <  2:193781/139‑1 (MQ=255)
                                                                                                ggACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAgcctgg                                                      >  2:407395/1‑149 (MQ=255)
                                                                                                 gacgacGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTgg                                                                                                                    >  8:253087/1‑87 (MQ=255)
                                                                                                 gacgacGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAgcctg                                                       <  1:251178/147‑1 (MQ=255)
                                                                                                  acgacgAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAgcctgggc                                                    <  8:166531/149‑1 (MQ=255)
                                                                                                   cgacgaAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAgcctgggcc                                                   >  5:422177/1‑149 (MQ=255)
                                                                                                      cgaAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAgcctgggcctgg                                                <  7:123544/149‑1 (MQ=255)
                                                                                                         aGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCa                                             <  5:11610/149‑1 (MQ=255)
                                                                                                          gCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAg                                            <  7:15786/149‑1 (MQ=255)
                                                                                                           cccTGGTACGTCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGc                                           <  1:367635/149‑1 (MQ=255)
                                                                                                              tGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAg                                        >  5:388183/1‑149 (MQ=255)
                                                                                                                gTACGCCTGATCGCCCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCt                                      <  4:427215/149‑1 (MQ=255)
                                                                                                                gTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGc                                       <  6:388183/148‑1 (MQ=255)
                                                                                                                    gCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGa                                   >  3:287126/1‑148 (MQ=255)
                                                                                                                     ccTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGa                                   >  4:190612/1‑147 (MQ=255)
                                                                                                                     ccTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGa                                   >  7:90371/1‑147 (MQ=255)
                                                                                                                     ccTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGa                                   >  2:367143/1‑147 (MQ=255)
                                                                                                                            gACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGa                           <  1:259862/149‑1 (MQ=255)
                                                                                                                            gACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGa                           <  1:425849/148‑1 (MQ=255)
                                                                                                                              ccGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATgc                        <  6:388221/149‑1 (MQ=255)
                                                                                                                               cgcgTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATgca                       <  5:241780/149‑1 (MQ=255)
                                                                                                                               cgcgTGCTGATGTTCTACATCCGCACGGCCGAAAATCTGAAGCGCAACTCGGTGTGGCGAGAATACCTTGAATGCGGGCTGGACTACCTGAAGCAGGTTATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATgca                       <  6:422177/149‑1 (MQ=255)
                                                                                                                                gcgTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATgcag                      >  1:433310/1‑149 (MQ=255)
                                                                                                                                  gTGCGGATGTAGAACATCAGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgacga                                                               <  2:32220/107‑1 (MQ=255)
                                                                                                                                  gTGCGGATGTAGAACATCAGCACGGCCGACAAGCTGCAGAGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATcgacga                                                               >  1:32220/1‑107 (MQ=255)
                                                                                                                                   ggatGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCAt                   <  4:387281/146‑1 (MQ=255)
                                                                                                                                      tGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCAtgtg                >  4:49622/1‑149 (MQ=255)
                                                                                                                                      tGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCAtgtg                >  6:85231/1‑149 (MQ=255)
                                                                                                                                         tGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTc             <  2:302117/149‑1 (MQ=255)
                                                                                                                                           ttCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGa           >  4:72792/1‑149 (MQ=255)
                                                                                                                                             cTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCCGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGAc          >  8:239836/1‑148 (MQ=255)
                                                                                                                                             cTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGAc          <  1:407395/148‑1 (MQ=255)
                                                                                                                                             cTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGAc          >  4:84088/1‑148 (MQ=255)
                                                                                                                                                cATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAg       >  7:187697/1‑148 (MQ=255)
                                                                                                                                                    cGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGa  <  8:90371/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATCGGCGTTATCGCCACCGACAAGGGCTGGAACCTGTACGTGTGCGGTAACGGCGGCATGCGCCCGCGGCATGCCGAGCTGTTCGCCACCGACCTGGACGACGAAGCCCTGGTACGCCTGATCGACCGCGTGCTGATGTTCTACATCCGCACGGCCGACAAGCTGCAGCGCACCTCGGTGTGGCGCGAAAACCTGGAAGGCGGGCTGGACTACCTGAAGCAGGTGATCATCGACGACAGCCTGGGCCTGGCCAGCGAGCTGGAAAGCCAGATGCAGCATGTGGTCGACCAGTATGA  >  NC_002947/1905836‑1906132

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: