Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,807 A→T 7.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 7.8% 173.3 / 3.1 61K120* (AAG→TAG) PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base T (3/2);  total (32/29)
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

CCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTT  >  NC_002947/1899662‑1899953
                                                                                                                                                 |                                                                                                                                                  
ccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGt                                                                                                                                                 <  8:172689/149‑1 (MQ=255)
ccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGt                                                                                                                                                 >  7:172689/1‑149 (MQ=255)
  agccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                               >  3:22521/1‑149 (MQ=255)
  agccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                               >  2:319171/1‑149 (MQ=255)
   gccgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTgg                                                                                                                                              <  5:106213/149‑1 (MQ=255)
    ccgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGc                                                                                                                                             <  3:117077/149‑1 (MQ=255)
      gccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGcc                                                                                                                                            >  2:361802/1‑148 (MQ=255)
                aCCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTAt                                                                                                                                 <  8:262630/149‑1 (MQ=255)
                aCCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTAt                                                                                                                                 >  4:131505/1‑149 (MQ=255)
                  cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATcc                                                                                                                               <  8:22437/149‑1 (MQ=255)
                  cAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGAGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATcc                                                                                                                               <  6:103020/149‑1 (MQ=255)
                     cATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTa                                                                                                                            <  7:286825/149‑1 (MQ=255)
                         gACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTAcgcg                                                                                                                        >  4:35720/1‑149 (MQ=255)
                             aGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                               <  5:282286/122‑1 (MQ=255)
                             aGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                               >  6:282286/1‑122 (MQ=255)
                                ccATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATc                                                                                                                 <  1:361802/149‑1 (MQ=255)
                                     gCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                            <  8:98577/149‑1 (MQ=255)
                                      cTGCCACCGCCGAGCGTCTGCGCGGCATCATCCACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                            >  4:207674/1‑148 (MQ=255)
                                         ccACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCaa                                                                                                         >  3:127942/1‑148 (MQ=255)
                                           aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTc                                                                                                      >  2:132691/1‑149 (MQ=255)
                                               ccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCgcgc                                                                                                  >  8:12043/1‑149 (MQ=255)
                                               ccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCgcgc                                                                                                  <  2:251792/149‑1 (MQ=255)
                                                 gAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACATGGCCAACAATTTGGACAAGGGGTTTAACCGTACGCGCCAAAGCGAGTCCAACACGCCCCt                                                                                                >  7:54244/1‑149 (MQ=255)
                                                     gTCTGCGCGGCATCATCGACCAGCAGGGCCGGGATGCCGTGGCGGTGTACGTGTCCGGGCAAATGTCGGTGGAGTCCCAGTACCTGGCCAACAAGTTGGCTAATGGCTTTATCCGTACGCGCCACATCGAGTCGAACTGGCGCCTGTGc                                                                                            <  4:127942/149‑1 (MQ=255)
                                                          cgcgGCATCATCCACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgc                                                                                       <  3:207674/149‑1 (MQ=255)
                                                               catcatCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcg                                                                                  >  6:271719/1‑149 (MQ=255)
                                                                   atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCtgt                                                                                              <  1:147169/133‑1 (MQ=255)
                                                                   atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCtgt                                                                                              >  2:147169/1‑133 (MQ=255)
                                                                   atcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgccag                                                                              >  1:29530/1‑149 (MQ=255)
                                                                        ccAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTg                                                                         >  6:180455/1‑149 (MQ=255)
                                                                         cAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                        <  2:29530/149‑1 (MQ=255)
                                                                             aCGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTAc                                                                    <  4:291685/149‑1 (MQ=255)
                                                                                gCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAg                                                                 >  2:150613/1‑149 (MQ=255)
                                                                                    ggATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCACCTCGGGCCGGTGCATGGCCCGCGCCAGCAGTGGCGACAAGgtgt                                                             >  8:151128/1‑149 (MQ=255)
                                                                                       tGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGc                                                          <  1:150613/149‑1 (MQ=255)
                                                                                            tGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGgc                                                     <  3:35720/149‑1 (MQ=255)
                                                                                                   gTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATg                                              >  1:323386/1‑149 (MQ=255)
                                                                                                      cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGcc                                           <  2:323386/149‑1 (MQ=255)
                                                                                                            cGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCt                                      <  3:131505/148‑1 (MQ=255)
                                                                                                            cGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTg                                     >  6:267612/1‑149 (MQ=255)
                                                                                                             gggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTgg                                    >  1:319687/1‑149 (MQ=255)
                                                                                                              ggCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGc                                   <  4:22521/149‑1 (MQ=255)
                                                                                                                   aTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                        >  7:146337/1‑107 (MQ=255)
                                                                                                                   aTGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                        <  8:146337/107‑1 (MQ=255)
                                                                                                                      tCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                           <  1:319171/149‑1 (MQ=255)
                                                                                                                           ggAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACt                      >  7:317652/1‑149 (MQ=255)
                                                                                                                           ggAGGCCCAGTACCAGGCCAACAAGTTGGCCAAGGGCTTTAAACGTACGCGCCACATCGAGTCCAACTCGCGCGTTTGCATGGCCAGCGCCAGCAGTGGATACAAGAAGACGCTTGCCGCGGATGGCCCAACTGCCAGATACCAGGACt                      >  6:56850/1‑149 (MQ=255)
                                                                                                                             cggCCCACTACCTGGCCAACAACTTGGCCAATGGCTTTATCCGTACGCGCCAAATCGATTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTc                    <  8:54244/148‑1 (MQ=255)
                                                                                                                                 ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGtt                               >  7:238545/1‑134 (MQ=255)
                                                                                                                                 ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGtt                               <  8:238545/134‑1 (MQ=255)
                                                                                                                                 ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGc                <  6:235767/149‑1 (MQ=255)
                                                                                                                                 ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGc                >  3:241627/1‑149 (MQ=255)
                                                                                                                                 ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGCCAGTTACCAGCACTACGac                 >  8:70867/1‑147 (MQ=255)
                                                                                                                                  cAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCa               >  4:106163/1‑149 (MQ=255)
                                                                                                                                     tACCTGGCCAACTAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGc            <  7:70867/149‑1 (MQ=255)
                                                                                                                                     cccttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGGTACCAGGACTTCGAGCACGc            >  3:57561/5‑149 (MQ=255)
                                                                                                                                       cttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGc                                   >  7:282078/3‑124 (MQ=255)
                                                                                                                                       cttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGc                                   <  8:282078/122‑1 (MQ=255)
                                                                                                                                       cttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGcc           >  8:279160/3‑148 (MQ=255)
                                                                                                                                           gCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCt                                      >  5:325115/1‑117 (MQ=255)
                                                                                                                                           gCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCt                                      <  6:325115/117‑1 (MQ=255)
                                                                                                                                             caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAt    >  4:11061/1‑149 (MQ=255)
                                                                                                                                               acaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAttt  <  5:180455/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTT  >  NC_002947/1899662‑1899953

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