Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,899,025 T→C 6.0% intergenic (+53/‑467) PP_t18 → / → PP_1703 tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,899,0250TC 6.0% 291.0 / 4.7 99intergenic (+53/‑467)PP_t18/PP_1703tRNA‑Asp/assimilatory nitrate reductase/sulfite reductase
Reads supporting (aligned to +/- strand):  ref base T (43/50);  new base C (3/3);  total (46/53)
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.40e-01

Only 100 of 104 total aligned reads displayed.

CCACCACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCA  >  NC_002947/1898877‑1899173
                                                                                                                                                    |                                                                                                                                                    
ccaccaCTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCtt                                                                                                                                                     >  2:1089361/1‑150 (MQ=32)
      cTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTc                                                                                                                                               <  1:779699/150‑1 (MQ=35)
       tCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTc                                                                                                                                               <  1:2056148/149‑1 (MQ=35)
       tCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTc                                                                                                                                               >  2:2056148/1‑149 (MQ=35)
       tCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCg                                                                                                                                              >  1:131167/1‑150 (MQ=35)
         gATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCt                                                                                                                                            <  1:1089361/150‑1 (MQ=35)
             tataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCttttt                                                                                                                                        <  2:2118699/150‑1 (MQ=35)
             tataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAAAAAGCGACCCTAAGGTCGCtttct                                                                                                                                        <  1:1237742/150‑3 (MQ=255)
              ataGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACATTAGGGTCGCtttttt                                                                                                                                       <  2:1252951/150‑1 (MQ=25)
                         gTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCCCCACTTGGGTCGCTTTTTTTATTTCTGcat                                                                                                                            >  1:252/1‑150 (MQ=14)
                          tAGTTCAGTCGGTTAGACAGACGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGcatc                                                                                                                           <  1:120965/150‑1 (MQ=255)
                           aGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGcatca                                                                                                                          <  2:2138139/150‑1 (MQ=255)
                             ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGTGTCGCtttt                                                                                                                                         >  2:1489485/1‑133 (MQ=18)
                             ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGTGTCGCtttt                                                                                                                                         <  1:1489485/133‑1 (MQ=18)
                             ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATg                                                                                                                        <  2:131167/150‑1 (MQ=255)
                             ttCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATg                                                                                                                        <  2:1017456/150‑1 (MQ=255)
                               cAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCa                                                                                                                      >  1:629262/1‑150 (MQ=255)
                                  tCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATc                                                                                                                   <  1:369028/150‑1 (MQ=255)
                                    ggTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGt                                                                                                                 <  2:2234813/150‑1 (MQ=255)
                                     gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGtt                                                                                                                >  2:549628/1‑150 (MQ=255)
                                     gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCACGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCCTAAGGTc                                                                                                                                               <  1:2223832/119‑1 (MQ=255)
                                     gTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCACGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCCTAAGGTc                                                                                                                                               >  2:2223832/1‑119 (MQ=255)
                                          aaTACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCAAGGGACCACTGAACACCCGGAAGTAAAAAAGTAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGc                                                                                                           <  1:2034169/150‑1 (MQ=255)
                                            tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTa                                                                                                                                                    >  1:1674215/1‑107 (MQ=33)
                                            tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTa                                                                                                                                                    <  2:1674215/107‑1 (MQ=33)
                                            tACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGccc                                                                                                         >  1:1153059/1‑150 (MQ=255)
                                             aCCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCt                                                                                                        <  2:1590479/150‑1 (MQ=255)
                                               cGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGc                                                                                                      >  2:922909/1‑150 (MQ=255)
                                                 gCCTGTCAGGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTACAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaa                                                                                                    >  2:515967/1‑150 (MQ=255)
                                                 gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaa                                                                                                    >  1:1720981/1‑150 (MQ=255)
                                                 gCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaa                                                                                                    >  2:189125/1‑150 (MQ=255)
                                                  ccTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaac                                                                                                   >  1:1932053/1‑150 (MQ=255)
                                                    tGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGcaacaa                                                                                                 >  1:2193980/1‑150 (MQ=255)
                                                     gTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAg                                                                                                >  1:393025/1‑150 (MQ=255)
                                                      tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgg                                                                                               <  2:90954/150‑1 (MQ=255)
                                                      tCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATAATGCAATCGTTTACGCCCTGCAACAAgg                                                                                               <  2:1818606/150‑1 (MQ=255)
                                                       cACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAggg                                                                                              <  1:189125/150‑1 (MQ=255)
                                                       cACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAggg                                                                                              <  1:1860792/150‑1 (MQ=255)
                                                          gCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                  <  1:205266/111‑1 (MQ=255)
                                                          gCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTAttt                                                                                                                                  >  2:205266/1‑111 (MQ=255)
                                                                  tCGCGGGTTCGAGTCCCCTCCGCTGCGCCATATTCCACCTCCAGGTCCCACTGAACACCCGCAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCTTCATGCAATCGTTTCCGCCCTGCAACAATGGGCTAGCTgcgc                                                                                   <  1:849039/150‑1 (MQ=21)
                                                                   cgcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGcc                                                                                  <  1:2225854/150‑1 (MQ=255)
                                                                    gcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCtt                                                                                                                                                     <  1:1716917/82‑1 (MQ=35)
                                                                    gcgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCtt                                                                                                                                                     >  2:1716917/1‑82 (MQ=35)
                                                                     cgGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACCTACCGGCCTCGCAGCGTTCTCGGAATTAAAAAAGAACGAGATCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCAc                                                                                <  1:523413/150‑1 (MQ=25)
                                                                      gggTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCCTAAGGTCGctttcttttt                                                                                                                                    >  2:1111088/1‑95 (MQ=14)
                                                                      gggTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCCTAAGGTCGctttcttttt                                                                                                                                    <  1:1111088/97‑3 (MQ=14)
                                                                          tCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtcttc                                                                           <  2:324832/150‑1 (MQ=255)
                                                                             aGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCttcttc                                                                        <  2:393025/150‑1 (MQ=255)
                                                                               tCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGACGTAAAAAAGAACGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTccc                                                                      <  1:2173222/150‑1 (MQ=255)
                                                                                 ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGGAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAg                                                                                                >  1:1696801/1‑122 (MQ=255)
                                                                                 ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGGAGTAAAAAAGAAAGCGACCTTAGGGTCGCTCTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAg                                                                                                <  2:1696801/122‑1 (MQ=255)
                                                                                 ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTccc                                                                      <  2:750750/148‑1 (MQ=255)
                                                                                 ccGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTccc                                                                      >  1:750750/1‑148 (MQ=255)
                                                                                  cGTCGGCTGCGCAATATTCCAATTCAAGGGCCGACTGAACAGCCGGAAGTAAAAAAGAAAACGACCTTGGGGTCGCTTATTTTATTGCTGCATCATGCAATCGTTTACACACTGCAGCAAGGGGCTAGCTGCGTCACTCTTCTTcccctg                                                                   >  2:745674/1‑147 (MQ=255)
                                                                                       gCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTc                                                              <  1:1633570/150‑1 (MQ=255)
                                                                                       gCTGCGCCATATTCCACTTCCACGGCCCACTGACCACCCGCAAGTAAAAAAGAATGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTc                                                              <  2:1360439/150‑1 (MQ=255)
                                                                                            gcCCTATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCGGCAACAAGGGGCTAGCTGCGCCACTCCTCTTCCCGCAGCGTCAGccc                                                         <  2:629262/150‑3 (MQ=255)
                                                                                             ccATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCAcc                                                        >  1:576103/1‑150 (MQ=255)
                                                                                              cATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACATCAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGAGACGCTCTTCTTCCCGCCGCGTCTGCActt                                                       >  1:1071701/1‑148 (MQ=25)
                                                                                                   tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAACTCGGGGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                  <  1:613172/150‑1 (MQ=38)
                                                                                                   tCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAAc                                                  >  1:40318/1‑150 (MQ=255)
                                                                                                    ccACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAggg                                                                                              <  2:2083268/105‑1 (MQ=255)
                                                                                                    ccACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAggg                                                                                              >  1:2083268/1‑105 (MQ=255)
                                                                                                     cACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATc                                                                                                                   <  2:1781967/83‑1 (MQ=255)
                                                                                                     cACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATc                                                                                                                   >  1:1781967/1‑83 (MQ=255)
                                                                                                      aCTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCt                                                                                                        <  1:688596/93‑1 (MQ=255)
                                                                                                      aCTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCg                                               >  2:2183100/1‑150 (MQ=255)
                                                                                                      aCTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCg                                               >  1:1090344/1‑150 (MQ=255)
                                                                                                      aCTTCCAGGGACCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCt                                                                                                        >  2:688596/1‑93 (MQ=255)
                                                                                                        ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCGCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTc                                             >  1:1251644/1‑150 (MQ=255)
                                                                                                        ttCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTc                                             >  1:1019660/1‑150 (MQ=255)
                                                                                                           cAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCCACGGGGGTCGACGTGTACTACTCTTCTTCCTGCATCGTCAGCACCTCAAACCCGTCCTc                                          >  2:1814673/1‑150 (MQ=25)
                                                                                                              ggCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGAt                                       <  2:2193980/150‑1 (MQ=255)
                                                                                                                  cACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCAcc                                   <  1:922909/150‑1 (MQ=255)
                                                                                                                   aCTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTGATTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCg                                  >  1:1009096/1‑150 (MQ=255)
                                                                                                                    cTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaaa                                                   >  1:709693/1‑132 (MQ=255)
                                                                                                                    cTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCaaa                                                   <  2:709693/132‑1 (MQ=255)
                                                                                                                    cTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCgg                                 <  1:549628/150‑1 (MQ=255)
                                                                                                                    cTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCgg                                 <  2:821454/150‑1 (MQ=255)
                                                                                                                      gAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGccc                                                                                                         >  1:1396669/1‑76 (MQ=255)
                                                                                                                      gAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGccc                                                                                                         <  2:1396669/76‑1 (MQ=255)
                                                                                                                       aaCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtctt                                                                            >  1:1050768/1‑104 (MQ=255)
                                                                                                                       aaCACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACtctt                                                                            <  2:1050768/104‑1 (MQ=255)
                                                                                                                            ccGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGGTTACGCCCTGCAACAAGGGGCTAGCTGCGCAACTCTTCTTCCCGCAGCGTCAGAAACTCAAACCCGTCCTAGCTCACCGGCATCTCCt                         >  1:247230/1‑150 (MQ=255)
                                                                                                                             cGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTg                        <  2:1932053/150‑1 (MQ=255)
                                                                                                                              ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgggg                                                                                             >  1:1739542/1‑80 (MQ=255)
                                                                                                                              ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAgggg                                                                                             <  2:1739542/80‑1 (MQ=255)
                                                                                                                              ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCGAGTCTTCTTCCCGCAGCGTCAGCACCTCAAACCGGTACTCGATCACGGGCATCTCATgg                       >  1:947097/1‑150 (MQ=255)
                                                                                                                              ggAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTgg                       <  1:2183100/150‑1 (MQ=255)
                                                                                                                                     aaaaaGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTTCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATgg                >  2:2093014/1‑150 (MQ=255)
                                                                                                                                          gaaagCGACCTTAGGGTCGCTTTTTTTTTTTTTTCTTCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg           >  2:1413450/1‑150 (MQ=255)
                                                                                                                                          gaaagCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAggg           <  2:576103/150‑1 (MQ=255)
                                                                                                                                             agCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAg        >  2:1990371/1‑150 (MQ=255)
                                                                                                                                             agCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAg        >  1:505040/1‑150 (MQ=255)
                                                                                                                                                 aCCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAg                                                                  >  2:1507050/1‑88 (MQ=255)
                                                                                                                                                 aCCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAg                                                                  <  1:1507050/88‑1 (MQ=255)
                                                                                                                                                   cTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCa  <  2:1090344/150‑1 (MQ=255)
                                                                                                                                                    ttAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCttcttc                                                                        <  2:492735/79‑1 (MQ=255)
                                                                                                                                                    ttAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCttcttc                                                                        >  1:492735/1‑79 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCACCACTCGATGTATAGCGCAGCGGTAGTTCAGTCGGTTAGAATACCGGCCTGTCACGCCGGGGGTCGCGGGTTCGAGTCCCGTCCGCTGCGCCATATTCCACTTCCAGGGCCCACTGAACACCCGGAAGTAAAAAAGAAAGCGACCTTAGGGTCGCTTTTTTTATTTCTGCATCATGCAATCGTTTACGCCCTGCAACAAGGGGCTAGCTGCGCCACTCTTCTTCCCGCAGCGTCAGCACCTCAAACCCGTCCTCGATCACCGGCATCTCCTGGGACATGGGAGGGCAGTGCTCA  >  NC_002947/1898877‑1899173

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