Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,286,944 T→G 6.7% intergenic (+418/+122) narI → / ← rttR nitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,286,9440TG 6.7% 748.1 / 23.2 275intergenic (+418/+122)narI/rttRnitrate reductase 1, gamma (cytochrome b(NR)) subunit/rtT sRNA, processed from tyrT transcript
Reads supporting (aligned to +/- strand):  ref base T (66/181);  new base G (20/3);  total (91/184)
Fisher's exact test for biased strand distribution p-value = 1.60e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 277 total aligned reads displayed.

GTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCG  >  NC_000913/1286795‑1287092
                                                                                                                                                     |                                                                                                                                                    
 tCGAACGGTGGTCTAATCTTCTGCTCCTTCCCCGTTTGGGCAGATTATTTGATTTCGTATTCGTTTGAGTTTTGCGGTGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTg                                                                                                                                                     <  1:121363/150‑1 (MQ=14)
   gAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCg                                                                                                                                                   <  1:2057334/150‑1 (MQ=255)
    aaCCTTGGTCTAAGCTTCTCAACCTTCCCCGCTTGGGCAGAATATTTTATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCgg                                                                                                                                                  <  1:5572156/150‑1 (MQ=255)
      ccTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCa                                                                                                                                                >  1:1072652/1‑150 (MQ=255)
       cTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCaa                                                                                                                                               <  1:6238066/150‑1 (MQ=255)
        ttGGTCGAAGCTTCTCATCCTTCCCGGCTTGGGCATAATTTTTGATTGCGGTTTCGTTTGAGATTTCCGGGGCTTTTGAAAGTGTTGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGCTGGCAAc                                                                                                                                              <  1:4036267/150‑1 (MQ=14)
         tGGTCGAAGCTTCTCATCCTTCCCCGCTGGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCTTTTGAAAGTGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTGCGGGTCGTTGCCGGCGGCACCg                                                                                                                                             >  2:6172038/1‑150 (MQ=14)
            tCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCGGGGGGTTTTGAAAGTGATGGTGGTTGGGTATTGTTTTGTAAGAGATTCGATCTTCGGCGTTCTGTTCGTTGACTGCGGctaaggtc                                                                                                                                          >  2:4730436/1‑143 (MQ=9)
              gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGctcac                                                                                                                                        >  1:327989/1‑148 (MQ=255)
              gAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGCGTTATGGAGGGGAGGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCGGCGGGTCGTTGCCGGCGGCAACGcgcgc                                                                                                                                        >  2:3334538/1‑146 (MQ=14)
                aGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                                                                                                      <  2:1664648/150‑1 (MQ=255)
                aGCTTCTCAGCCTTCCCCGCTTGGGCAGAATACTTGAGTGCGTATTAGTTTGAGAATTCCTGGGCTTTTGAAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCtctctc                                                                                                                                      <  2:1625616/150‑1 (MQ=14)
                 gCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGGTTGAGAATTCCGGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCTCTCTCg                                                                                                                                     >  2:2019518/1‑150 (MQ=255)
                    tctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAGTTCAGGGGCTTTTGAAAGTGATGGTGGTCGGGGGAGCACTACTCAGAGCTCCGCTCTACGACCTACGGGTCGTTTGCAGGGGCATCGATCTCTCGcgg                                                                                                                                  >  1:4620524/1‑148 (MQ=11)
                     ctcATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGa                                                                                                                                 >  1:616094/1‑150 (MQ=255)
                        aTCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCTCTGCGCGCTTCGCCCTGCGGGTCGTTGCCGGCGGCAACGCACGCGAGCTGACGc                                                                                                                              >  2:5598059/1‑150 (MQ=14)
                            gtCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGa                                                                                                                          <  1:181110/149‑1 (MQ=255)
                             tCCCCGCTTGGGCAGAATATTTGAGTGCGGAGGCGTTTGAGAATTCCGGGTCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                         <  1:1586801/150‑1 (MQ=14)
                             tCCCCGCGTGGGCAGCTTTCTTGATTGGGTATTAGTTTGAGCATTCGGTGTCTTTTGTAAGGGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAg                                                                                                                         <  2:1552208/150‑1 (MQ=11)
                              ccccGCTTGGGCAGATTATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGc                                                                                                                        <  2:7042325/150‑1 (MQ=255)
                                  ggttGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGaa                                                                                                                    <  1:6340376/148‑1 (MQ=255)
                                    ttGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCTCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAAcc                                                                                                                  >  2:2850897/1‑150 (MQ=255)
                                       ggCAGAATCTTTGATTGCGGATTCGTTTGAGAATTCCGGGTCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTa                                                                                                               <  2:6523987/150‑1 (MQ=255)
                                        gCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCAGCGGCAACGCACTCACGCGGACGCTCGAGCCGAACCTTAg                                                                                                              >  1:3057075/1‑150 (MQ=255)
                                         cAGAATATTTGATTGCGGATTCGTTTGAGTATTCCGTGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGt                                                                                                             <  1:1115420/150‑1 (MQ=255)
                                          aGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTAAGAGCTGCGCGCTTCGCCCTTCGGGTCGTTGCGTGCGTCAAACCTCTCTATCTGGCGCTCGAGACTAACCTTAgac                                                                                                            >  2:4412343/1‑148 (MQ=11)
                                            aagatTTGAGTGCGGAGGCGTTTGAGAATTCCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGa                                                                                                          <  2:889696/147‑1 (MQ=14)
                                             atatTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGaa                                                                                                         <  1:4268785/150‑1 (MQ=255)
                                             atatTTGATTGCGGATTCGTTTGAGAATTACGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTAGCCCTTCGGGTCGTTGCCAGCGGCAACCCTCTATCGCGGACGCTCGAGCCGAACCATAGTCGaa                                                                                                         >  1:3816319/1‑150 (MQ=14)
                                              tatTTGATTGCGGATTCGTTTGAGAAGTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAg                                                                                                        <  2:5169853/150‑1 (MQ=255)
                                               atTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGc                                                                                                       <  1:4458071/150‑1 (MQ=255)
                                                  tGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTtc                                                                                                    <  1:3906956/150‑1 (MQ=255)
                                                   gATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGTTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTTGTCGAAGCTtct                                                                                                   <  1:2491641/150‑1 (MQ=255)
                                                      tGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCGAGTCGAACCTTGGTCGAAGCTTCTCAt                                                                                                >  2:2430284/1‑150 (MQ=12)
                                                       gCGGATTCGTTTGAGTATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATc                                                                                               <  2:1589305/150‑1 (MQ=255)
                                                              cGTTTGAGTATTGCGGGGCTTTTGCAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTcccc                                                                                        <  2:1242651/150‑1 (MQ=255)
                                                                tttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGGCGAACCTTCGTCGAAGCATCTCATCCTTCCCCGc                                                                                      >  2:4715211/1‑150 (MQ=255)
                                                                 ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGGCGAACCTTTGTCGAAGCTTCTCATCCTTACCCGCa                                                                                     >  2:136359/1‑150 (MQ=255)
                                                                 ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                     <  2:616094/150‑1 (MQ=255)
                                                                 ttGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCa                                                                                     <  1:5933019/150‑1 (MQ=255)
                                                                  tGAGTAGTCCGGGGCTTTTGTAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                    <  1:5205548/150‑1 (MQ=255)
                                                                  tGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCAt                                                                                    <  1:6963151/150‑1 (MQ=255)
                                                                    agcaTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATgg                                                                                  <  1:775445/147‑1 (MQ=255)
                                                                       aTTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCa                                                                               <  1:7116324/150‑1 (MQ=255)
                                                                           cGGGGCTTTTGAAGGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCCGCAACGCTCTCTCCCTGACGATCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCCGTGGCAGaag                                                                           >  2:4484393/1‑149 (MQ=14)
                                                                           cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGTATTACTCAGCTCTTCGAGCGAAGCACTTCCGGTCGTAGCCTGCGGCAACGCTCTCTCGCTGTCGGTCAAGCCGAACCTTCGTCGGAGCTTCTCATCCTTCCCCGCATGGGccggag                                                                           >  2:1094551/1‑145 (MQ=11)
                                                                           cGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAAACCACACACGCTGACGCACGAGCCCAACAAAAGTAGAATCTACTAAACCATACCCCCAAGGGCCGaaa                                                                           >  1:705411/1‑149 (MQ=2)
                                                                               gCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCTCGATCCGAACCTTAGTCGAAGCATCTCATACTTCCCCCCATGGGCAGAATAttt                                                                       >  2:2803264/1‑150 (MQ=11)
                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGa                                                                     <  2:5977920/150‑1 (MQ=38)
                                                                                 ttttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCACGCTGACGCGCGAGCCGGACCCTAGTCGAAGCTTCTCATCCATCCCCGCATGGGCAGAATATTTGa                                                                     >  2:1315462/1‑150 (MQ=11)
                                                                                   ttGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCATCGGGTCGTTGCATGCGGCAACGCTCTCACGCTGACGCTCGAGCCGAACATAAGTCGAAGCTACTCATCCTTCCCCGCATGGGCAGAATATTTGAtt                                                                   >  2:4185308/1‑150 (MQ=14)
                                                                                     gAAAGTGATGGTGGTGGGGGAAGGATTACTCAGAGCTGCGCGCTTCGCCCTTCGGTTCGTTGCCTGCGGCAACGCTCTCTCGCTGACCCACGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAATATTAGATTGc                                                                 >  2:471428/1‑150 (MQ=14)
                                                                                          tGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGTGCTTCGCCCTTCGGGTCGTTGCCGGCGGCAACGCACTCTCGCTGACGCTCGATCCGAACCTTAGGCGCATCATCGCAACCAAACCCGCCGGGGCAGTATATTTGATTGCCGAtt                                                            >  2:5489652/1‑150 (MQ=14)
                                                                                              ggtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCt                                                        <  1:2720382/150‑1 (MQ=40)
                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGGCGCTCTAGCCGAACCTTCGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                       <  2:3517252/150‑1 (MQ=255)
                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                       <  1:2528582/150‑1 (MQ=255)
                                                                                               gtggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCtt                                                       <  1:2193331/150‑1 (MQ=255)
                                                                                                 ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                     <  1:1952371/150‑1 (MQ=255)
                                                                                                 ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                     <  1:3645528/150‑1 (MQ=255)
                                                                                                 ggtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTga                                                     <  1:1855265/150‑1 (MQ=255)
                                                                                                  gtggGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTgag                                                    <  2:4963373/150‑1 (MQ=255)
                                                                                                     ggggAAGGAGGACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCACCGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCGTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGt                                                 <  1:6172038/150‑1 (MQ=255)
                                                                                                      gggAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGtt                                                <  1:4598598/150‑1 (MQ=255)
                                                                                                        gAAGGATTACTCAGCGCCGCGCGCTTCGCCCTTCGGGTCGGTGCCTGCGGCAACGCTCTCGCGCTGACGCTAGAGCCGAACCTTAGACGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                              <  1:1524987/150‑1 (MQ=14)
                                                                                                        gAAGGAGTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCa                                              <  1:2017739/150‑1 (MQ=255)
                                                                                                         aaGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAg                                             <  1:1883452/150‑1 (MQ=255)
                                                                                                             aTTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGAc                                         <  1:4961783/150‑1 (MQ=255)
                                                                                                              ttACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGGGCCGAACCTTAGTCGCAGCGTCTCATCCTTCCCCGCCTGGGCCGGATATTTGATTGCGGATTCGCTGGGGAGTTCAgggccg                                        >  2:3866021/1‑147 (MQ=14)
                                                                                                                    aGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCGTAGTCGAAGCTTCTCATCCTTCCCCCCATGGGCAGAATATTTGATTGCGGATTCGCTTGACAGTTCAGGGACATATGaa                                  >  2:1052931/1‑150 (MQ=14)
                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                             <  1:2165786/58‑1 (MQ=35)
                                                                                                                     gcgcTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCt                                                                                                                             >  2:2165786/1‑58 (MQ=35)
                                                                                                                     gcgcTGCGCGCTTCGCCCTGCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                 <  2:3767377/150‑1 (MQ=255)
                                                                                                                     ccgcTGCGCGCTTCGCCCTTCGGGTCGTTTCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTCAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGaaa                                 <  1:4730436/149‑1 (MQ=255)
                                                                                                                        cTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTg                              <  2:5733589/150‑1 (MQ=255)
                                                                                                                         tGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGTCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGA                             <  1:1827855‑M1/150‑1 (MQ=255)
                                                                                                                         tGCGCCCTTCGCCCTTCGGGTCGTTGCCTCCGGCAACGCTCTCTCGCTGACGCTCCACACGAACCTTAGTCGAAGCTTCTCAACCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGA                             <  2:5232257‑M1/150‑1 (MQ=255)
                                                                                                                           cgcgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATG                           <  2:4786456‑M1/150‑1 (MQ=255)
                                                                                                                             cgcTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGT                         <  1:6513814‑M1/150‑1 (MQ=255)
                                                                                                                                ttCACCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAAGATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGT                      >  2:5036487‑M1/1‑150 (MQ=255)
                                                                                                                                 tCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGATCCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGTTTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGTTGGTGGTG                     <  2:6947503‑M1/150‑1 (MQ=255)
                                                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGG                    <  1:3091419‑M1/150‑1 (MQ=255)
                                                                                                                                  cGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGG                    <  2:183823‑M1/150‑1 (MQ=255)
                                                                                                                                   gCCCTTCGGGTCGTTGCCGGCGGCAACGCACTCTCGCGGACGCGCGCGCCCAACCGTAGTCGAAGCTTCTCATCCGTCCCCCCATGGGCCGACTAGGTGATTGCGGCGTCGCTGGGGAGTTCAGGGACGGTTGAAAGGGAGGGTgggggg                   >  2:5606588/1‑150 (MQ=255)
                                                                                                                                    cccTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCGGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCAACCTTACCCGCAGGGGCAGAATATTTGATTGCGGATTCGCTTGTGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGG                  >  2:343235‑M1/1‑150 (MQ=255)
                                                                                                                                       ttCGGGTCGTTGCCTGGTGCAAAGCTCTGTGGCCGTCGGTCCTGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGTCTTTTGAAAGTGATGGTGGTGGGGGAA               <  1:6831606‑M1/150‑1 (MQ=255)
                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA               <  2:1542666‑M1/150‑1 (MQ=255)
                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGCCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA               <  1:7203702‑M1/150‑1 (MQ=255)
                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGGGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA               <  1:5830529‑M1/150‑1 (MQ=255)
                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA               <  2:3057075‑M1/150‑1 (MQ=255)
                                                                                                                                       ttCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAA               <  1:1642441‑M1/150‑1 (MQ=255)
                                                                                                                                           ggTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAT           <  2:6731378‑M1/150‑1 (MQ=255)
                                                                                                                                           gagcgttGCCGCAGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGAT           >  2:5769173‑M1/4‑150 (MQ=255)
                                                                                                                                            gTCGTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATT          <  2:662809‑M1/150‑1 (MQ=255)
                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa         <  1:3402461/150‑1 (MQ=255)
                                                                                                                                             tCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTa         <  1:2851074/150‑1 (MQ=255)
                                                                                                                                               gTTGCCTGCGGCAACGCTCTCTCGCTGTCGCTCTAGCCGAACCTTTGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACt       <  2:409871/150‑1 (MQ=255)
                                                                                                                                                ttGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTc      >  2:6358534/1‑150 (MQ=255)
                                                                                                                                                 tGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCa     <  1:2986627/150‑1 (MQ=255)
                                                                                                                                                   ccTGCGGCAACGCTCTCTCTCTGACGCTCGAGCCGAACCTTAGTCTAATCTTCTCATCCTTCCCCGCATGGGCAGAATATTTTATTGCGGATTCGCTTGAGAGTTCATGTAATTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc   <  1:6547066/150‑1 (MQ=17)
                                                                                                                                                   ccTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAgc   <  1:2800998/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTCGAACCTTGGTCGAAGCTTCTCATCCTTCCCCGCTTGGGCAGAATATTTGATTGCGGATTCGTTTGAGAATTCCGGGGCTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCGCTGCGCGCTTCGCCCTTCGGGTCGTTGCCTGCGGCAACGCTCTCTCGCTGACGCTCGAGCCGAACCTTAGTCGAAGCTTCTCATCCTTCCCCGCATGGGCAGAATATTTGATTGCGGATTCGCTTGAGAGTTCAGGGACTTTTGAAAGTGATGGTGGTGGGGGAAGGATTACTCAGCG  >  NC_000913/1286795‑1287092

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