Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,807 A→T 6.8% K120* (AAG→TAG)  PP_1703 → assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,8070AT 6.8% 211.2 / 2.7 68K120* (AAG→TAG) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base T (3/2);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCT  >  NC_002947/1899662‑1899955
                                                                                                                                                 |                                                                                                                                                      
ccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGt                                                                                                                                                     <  3:177192/149‑1 (MQ=255)
        cAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAg                                                                                                                                            <  8:280036/149‑1 (MQ=255)
             cgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCt                                                                                                                                        <  5:56208/148‑1 (MQ=255)
                 ccAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATc                                                                                                                                   <  5:181985/149‑1 (MQ=255)
                  cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCATCAAGTTGGCC‑AAGGGCTTTATcc                                                                                                                                  <  6:353913/149‑1 (MQ=255)
                  cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATcc                                                                                                                                  >  6:263869/1‑149 (MQ=255)
                  cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATcc                                                                                                                                  <  7:169959/149‑1 (MQ=255)
                    gCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGt                                                                                                                                <  1:104272/149‑1 (MQ=255)
                        ggACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTAcgc                                                                                                                            >  3:307303/1‑149 (MQ=255)
                        ggACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTAcgc                                                                                                                            >  2:231821/1‑149 (MQ=255)
                                       tGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTcc                                                                                                              >  3:290352/1‑148 (MQ=255)
                                          cACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTgg                                                                                                                                                  >  8:131335/1‑110 (MQ=255)
                                          cACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTgg                                                                                                                                                  <  7:131335/110‑1 (MQ=255)
                                          cACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACt                                                                                                          >  1:219642/1‑149 (MQ=255)
                                          cACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACt                                                                                                          >  2:187720/1‑149 (MQ=255)
                                           aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTc                                                                                                         >  3:104745/1‑149 (MQ=255)
                                           aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTc                                                                                                         >  3:114544/1‑149 (MQ=255)
                                           aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTc                                                                                                         >  8:12617/1‑149 (MQ=255)
                                           aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTc                                                                                                         >  5:351766/1‑149 (MQ=255)
                                             cgccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTcgc                                                                                                       >  3:280682/1‑149 (MQ=255)
                                                cgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGcc                                                                                                    >  5:186347/1‑149 (MQ=255)
                                                   gCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCtgt                                                                                                 <  4:280682/149‑1 (MQ=255)
                                                         gcgcGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                           >  5:197442/1‑149 (MQ=255)
                                                          cgcgGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGg‑c                                                                                         >  5:201502/1‑149 (MQ=255)
                                                              gCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGg‑ccagc                                                                                     <  6:201502/149‑1 (MQ=255)
                                                                 tcatcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGg‑ccagcgc                                                                                   >  7:233338/1‑148 (MQ=255)
                                                                   atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGg‑ccagcgccag                                                                                <  5:263869/149‑1 (MQ=255)
                                                                      gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCcagcag                                                                             >  7:279624/1‑149 (MQ=255)
                                                                        ccAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTg                                                                           <  4:290352/149‑1 (MQ=255)
                                                                         cAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AATGGCTTTATCCGTACGCGCCACATCGAGTCCTACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTgg                                                                          <  8:1044/149‑1 (MQ=255)
                                                                          aGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGc                                                                         <  6:186347/149‑1 (MQ=255)
                                                                          aGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGc                                                                         <  7:107537/149‑1 (MQ=255)
                                                                             aCGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCTCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTAc                                                                      <  6:351766/149‑1 (MQ=255)
                                                                              cGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGg‑ccagcgccag                                                                                >  6:33017/1‑138 (MQ=255)
                                                                              cGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGg‑ccagcgccag                                                                                <  5:33017/138‑1 (MQ=255)
                                                                                 cccGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGc                                                                  >  2:279953/1‑149 (MQ=255)
                                                                                     gATGCCGTGGCGCTGTACGTGTCAGGGCAAATGTCGCTGGAGGGCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAAATCGCGCGTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTc                                                              >  5:14781/1‑149 (MQ=255)
                                                                                          cGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCCAGCGCCAGCAGTGGCTACAAGCTGTCGCtt                                                          <  1:231821/149‑1 (MQ=255)
                                                                                          cGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                         <  8:233338/149‑1 (MQ=255)
                                                                                           gTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTgg                                                        >  2:34951/1‑149 (MQ=255)
                                                                                               cgcTGTACGTGGCCGGGCAAATGTAGCTGGCGGCCCAGGCCGTGGCCAACTAGTTGGCC‑CAGGGCTTTTTGCGTTCGCGCGACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTTCAAGCTGTCGCTTGGCGcc                                                    <  6:293873/149‑1 (MQ=255)
                                                                                                  tGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGAt                                                 <  4:104745/149‑1 (MQ=255)
                                                                                                    tACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATg                                                <  2:219642/148‑1 (MQ=255)
                                                                                                      cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGcc                                             <  4:307303/149‑1 (MQ=255)
                                                                                                       gtgtCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGccc                                            >  8:323921/1‑149 (MQ=255)
                                                                                                             gggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCGAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTg                                       >  1:76124/1‑149 (MQ=255)
                                                                                                             gggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATg                                                <  7:241287/139‑1 (MQ=255)
                                                                                                             gggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATg                                                >  8:241287/1‑139 (MQ=255)
                                                                                                             gggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTgg                                      <  1:34951/149‑1 (MQ=255)
                                                                                                              ggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGc                                     <  2:76124/149‑1 (MQ=255)
                                                                                                               gCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCa                                    <  6:197442/149‑1 (MQ=255)
                                                                                                               gCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCa                                    <  4:114544/149‑1 (MQ=255)
                                                                                                                cAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAg                                   <  7:12617/149‑1 (MQ=255)
                                                                                                                 aaaTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGt                                  >  6:215984/1‑149 (MQ=255)
                                                                                                                   aTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTa                                >  5:256909/1‑149 (MQ=255)
                                                                                                                      tCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                             <  8:279624/149‑1 (MQ=255)
                                                                                                                        gCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAg                            >  2:369227/1‑148 (MQ=255)
                                                                                                                        gCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAg                            <  5:215984/148‑1 (MQ=255)
                                                                                                                            gAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                           >  4:61564/1‑82 (MQ=255)
                                                                                                                            gAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                           <  3:61564/82‑1 (MQ=255)
                                                                                                                                 ccAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAg                   >  2:342192/1‑148 (MQ=255)
                                                                                                                                  cAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                             >  6:18639/1‑137 (MQ=255)
                                                                                                                                      ccttGGCCAACTTGTTGGCC‑AGGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                         >  5:237619/4‑105 (MQ=255)
                                                                                                                                      ccttGGCCAACTTGTTGGCC‑AGGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTg                                                         <  6:237619/102‑1 (MQ=255)
                                                                                                                                      ccttGGCCAACTTGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGcc             >  7:303201/4‑149 (MQ=255)
                                                                                                                                       cttGGCCAACTTGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGcc             >  5:329677/3‑148 (MQ=255)
                                                                                                                                        cTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGa           >  2:251980/1‑149 (MQ=255)
                                                                                                                                         tGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGAc          >  5:293157/1‑149 (MQ=255)
                                                                                                                                              aacaaGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt     >  4:61204/1‑149 (MQ=255)
                                                                                                                                              aacaaGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt     >  6:250262/1‑149 (MQ=255)
                                                                                                                                               acaaGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAttt    <  1:369227/149‑1 (MQ=255)
                                                                                                                                                 aaGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCt  <  6:293157/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCC‑AAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGG‑CCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCT  >  NC_002947/1899662‑1899955

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