Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,849 A→T 5.9% 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,8490AT 5.9% 603.5 / 20.2 186K120* (AAG→TAG) ‡PP_1703assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base A (69/106);  new base T (4/7);  total (73/113)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

Only 100 of 200 total aligned reads displayed.

CAGCCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTG  >  NC_002947/1899701‑1899998
                                                                                                                                                    |                                                                                                                                                     
cagccagcCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCcaacaa                                                                                                                                                      <  2:5239454/150‑1 (MQ=255)
cagccagcCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCcaacaa                                                                                                                                                      <  2:4856915/150‑1 (MQ=255)
 agccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAg                                                                                                                                                     <  1:5361600/150‑1 (MQ=255)
   ccagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGtt                                                                                                                                                   <  2:3250308/150‑1 (MQ=255)
    cagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                  <  1:5360970/150‑1 (MQ=255)
    cagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                  <  1:5251859/150‑1 (MQ=255)
    cagccGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                  <  2:4193375/150‑1 (MQ=255)
        cgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                  >  2:4533701/1‑146 (MQ=255)
        cgccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTg                                                                                                                                                  <  1:4533701/146‑1 (MQ=255)
         gccAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACTTGTTGGCCag                                                                                                                                             <  1:3225785/150‑2 (MQ=255)
            aGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCaa                                                                                                                                             >  2:933110/1‑147 (MQ=255)
            aGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAggg                                                                                                                                          <  2:4494507/150‑1 (MQ=255)
                cgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCttt                                                                                                                                      <  1:1892945/150‑1 (MQ=255)
                cgcACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCttt                                                                                                                                      <  2:675181/150‑1 (MQ=255)
                 gcACCAGCATGGACCAGGCCATCGCTGCCACCGCCTAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGAACCTGGCCAACAAGTTGGCCAAGGGCTTTa                                                                                                                                     >  1:4332485/1‑150 (MQ=255)
                   aCCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGcc                                                                                                                                               >  2:1649553/1‑138 (MQ=255)
                    ccAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATcc                                                                                                                                  <  1:4264381/150‑1 (MQ=255)
                    ccAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATcc                                                                                                                                  <  2:414292/150‑1 (MQ=255)
                          tGGACCCGGCCACGGGTGCCCCCGCCGTGCGGCTGCGCGGATTCTGGGCCCCGCCGGGCCCGGAGGCCGTGGCGCGGTAGGTGGCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTAcgc                                                                                                                            <  2:4285188/150‑1 (MQ=255)
                               cAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCcac                                                                                                                       >  2:5347198/1‑150 (MQ=255)
                                aGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCcaca                                                                                                                      >  1:5039834/1‑150 (MQ=255)
                                  gCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATc                                                                                                                    >  1:96298/1‑150 (MQ=255)
                                   ccATCGCTGCCACAGCCGAGCCTATGCGCGGAATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGACCGGGCAAATTTCGCTGGAGGCCCAGTACATGGCCAACAAGTTGGCCAAGGGATTTATCCGTACGCGCCACAACg                                                                                                                   >  2:790326/1‑150 (MQ=255)
                                     aTCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAg                                                                                                                 <  1:768430/150‑1 (MQ=255)
                                      tCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCca                                                                                                                        >  2:1854512/1‑142 (MQ=255)
                                       cGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                               <  1:3622870/150‑1 (MQ=255)
                                       cGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                               <  1:4972403/150‑1 (MQ=255)
                                            ccACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCTACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACt                                                                                                          >  2:4016296/1‑150 (MQ=255)
                                              aCCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTcg                                                                                                        >  2:2260890/1‑150 (MQ=255)
                                               ccgccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTcgc                                                                                                       >  2:4296471/1‑150 (MQ=255)
                                               ccgccgAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTcgc                                                                                                       >  2:3913285/1‑150 (MQ=255)
                                                    gAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCtg                                                                                                  >  2:1174664/1‑150 (MQ=255)
                                                     aGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCtgt                                                                                                 >  1:5340123/1‑150 (MQ=255)
                                                          cTGCGCGGCATCATCGACCAGCACGGCCCGGATGGCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATg                                                                                            <  1:2362411/150‑1 (MQ=255)
                                                          cTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATg                                                                                            <  2:2197517/150‑1 (MQ=255)
                                                           tGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATgg                                                                                           >  2:479204/1‑150 (MQ=255)
                                                            gcgcGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgc                                                                                          <  1:296128/150‑1 (MQ=255)
                                                              gcgGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgcca                                                                                        >  1:2420365/1‑150 (MQ=255)
                                                              gcgGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgcca                                                                                        >  1:3000861/1‑150 (MQ=255)
                                                              gcgGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgcca                                                                                        >  1:5148958/1‑150 (MQ=255)
                                                                ggCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagc                                                                                      <  1:2766938/150‑1 (MQ=255)
                                                                ggCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagc                                                                                      >  1:4983083/1‑150 (MQ=255)
                                                                  catcatCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgc                                                                                    <  1:2930822/150‑1 (MQ=255)
                                                                    tcatcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgcca                                                                                  <  1:1488676/150‑1 (MQ=255)
                                                                    tcatcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgcca                                                                                  <  2:1985774/150‑1 (MQ=255)
                                                                    tcatcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccagcgcca                                                                                  <  1:3928370/150‑1 (MQ=255)
                                                                       tcGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGgccggcgccaagc                                                                               >  2:4979930/1‑147 (MQ=255)
                                                                         gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             <  1:3713466/150‑1 (MQ=255)
                                                                         gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             <  2:2856265/150‑1 (MQ=255)
                                                                         gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             <  2:897191/150‑1 (MQ=255)
                                                                         gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             <  2:129802/150‑1 (MQ=255)
                                                                         gACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             <  1:2265139/150‑1 (MQ=255)
                                                                           ccAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                           <  2:2420365/150‑1 (MQ=255)
                                                                           ccAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTgg                                                                           <  2:5148958/150‑1 (MQ=255)
                                                                             aGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCt                                                                         >  1:173526/1‑150 (MQ=255)
                                                                             aGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCt                                                                         >  2:5396468/1‑150 (MQ=255)
                                                                              gCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTAt                                                                                                                                    <  2:4422159/90‑1 (MQ=255)
                                                                              gCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTAt                                                                                                                                    >  1:4422159/1‑90 (MQ=255)
                                                                              gCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATg                                                                                            <  1:1786561/130‑1 (MQ=255)
                                                                                aCGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACa                                                                      <  1:4092514/150‑1 (MQ=255)
                                                                                 cGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACaa                                                                     <  1:992273/150‑1 (MQ=255)
                                                                                    cccGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCt                                                                  >  1:2016633/1‑150 (MQ=255)
                                                                                        gATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCg                                                              >  1:3199830/1‑150 (MQ=255)
                                                                                        gATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCg                                                              <  2:1397424/150‑1 (MQ=255)
                                                                                         aTGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGc                                                             <  1:4708119/150‑1 (MQ=255)
                                                                                           gCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCtt                                                           <  2:1774247/150‑1 (MQ=255)
                                                                                              gTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGgc                                                        >  2:3807831/1‑150 (MQ=255)
                                                                                                 gcgcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCgcgc                                                                                                     >  1:4304475/1‑102 (MQ=255)
                                                                                                  cgcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCg                                                    <  2:3519951/150‑1 (MQ=255)
                                                                                                   gcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGa                                                   <  1:2845819/150‑1 (MQ=255)
                                                                                                   gcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGa                                                   <  1:76473/150‑1 (MQ=255)
                                                                                                   gcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGa                                                   >  2:4684937/1‑150 (MQ=255)
                                                                                                   gcTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGa                                                   >  1:1543327/1‑150 (MQ=255)
                                                                                                         cGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGccc                                             <  1:3908941/150‑1 (MQ=255)
                                                                                                                   cAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGt                                   >  1:4864894/1‑150 (MQ=255)
                                                                                                                       tGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTCCAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTAcc                               >  1:4377511/1‑150 (MQ=255)
                                                                                                                        gTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                              <  1:715259/150‑1 (MQ=255)
                                                                                                                        gTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                              <  1:1038261/150‑1 (MQ=255)
                                                                                                                           gCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGa                           <  2:432370/150‑1 (MQ=255)
                                                                                                                            cTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGAc                          <  2:5105125/150‑1 (MQ=255)
                                                                                                                               gAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTc                       <  1:3800860/150‑1 (MQ=255)
                                                                                                                               gAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTc                       <  2:2684749/150‑1 (MQ=255)
                                                                                                                                    ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                               >  1:1073619/1‑57 (MQ=255)
                                                                                                                                    ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTc                                                                                                               <  2:1073619/57‑1 (MQ=255)
                                                                                                                                    ccAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCa                  >  1:1670745/1‑150 (MQ=255)
                                                                                                                                    aaagccccTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCa                  >  2:3050916/7‑150 (MQ=255)
                                                                                                                                       gTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGc               >  1:1967500/1‑150 (MQ=255)
                                                                                                                                          ccTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGa            >  1:1241991/1‑150 (MQ=255)
                                                                                                                                           ttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTAc                                <  2:4929738/128‑1 (MQ=255)
                                                                                                                                           ttGGCCAACTTGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGAc           >  1:3411070/2‑150 (MQ=255)
                                                                                                                                            tGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACg          >  2:2114790/1‑150 (MQ=255)
                                                                                                                                              gCCAACTTGTTGGCCAGGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                              >  1:1513613/1‑128 (MQ=255)
                                                                                                                                              gCCAACTTGTTGGCCAGGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCa                              <  2:1513613/128‑1 (MQ=255)
                                                                                                                                                caacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAtt      >  1:571756/1‑150 (MQ=255)
                                                                                                                                                 aacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGTAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAttt     >  1:3777926/1‑150 (MQ=255)
                                                                                                                                                 aacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             <  1:2979554/78‑1 (MQ=255)
                                                                                                                                                 aacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGt                                                                             >  2:2979554/1‑78 (MQ=255)
                                                                                                                                                 aacaaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTAttt     >  2:624252/1‑150 (MQ=255)
                                                                                                                                                    aaGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTg  <  2:1549138/150‑1 (MQ=255)
                                                                                                                                                    aaGTTGGCCAAGGGCATTATCCGTTCGCGCCACATCGAGTCCAACTCGCGCCTGTTCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTg  <  2:4332485/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGCCAGCCGCCAGTGCGCACCAGCATGGACCAGGCCATCGCTGCCACCGCCGAGCGTCTGCGCGGCATCATCGACCAGCACGGCCCGGATGCCGTGGCGCTGTACGTGTCCGGGCAAATGTCGCTGGAGGCCCAGTACCTGGCCAACAAGTTGGCCAAGGGCTTTATCCGTACGCGCCACATCGAGTCCAACTCGCGCCTGTGCATGGCCAGCGCCAGCAGTGGCTACAAGCTGTCGCTTGGCGCCGATGGCCCACCTGGCAGTTACCAGGACTTCGAGCACGCCGACGTATTTCTG  >  NC_002947/1899701‑1899998

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: