Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	37.1%	114.7 / 13.5	71	S215G (AGT→GGT)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (36/8);  new base C (0/26);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556060‑3556357                                                                                                                                                      |                                                                                                                                                       ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTccc                                                                                                                                                     <  1:1139225/151‑1 (MQ=255)     agccaccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTcccc                                                                                                                                                    <  2:1107976/146‑1 (MQ=255)        caccTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAg                                                                                                                                               >  2:1194719/1‑150 (MQ=255)           ctcgcccgggcagcgcTTTTCCCGGCAGAAAAAAAAGCTTCGGGGGATATTGTCCCTTCGCACTACAAGCCGGGCGGTATTCCCCGGCCCCTGGTCGCGCGCAGGGTTATGAATATTTGGTTCAGGAAGGGAAAA‑A‑CTCCCCGCCAgtgt                                                                                                                                            <  2:795055/138‑1 (MQ=255)                   tGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGCCGCTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGAttt                                                                                                                                    >  2:744563/1‑150 (MQ=255)                   tGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGAttt                                                                                                                                    >  2:44283/1‑150 (MQ=255)                   tGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTc                                                                                                                                   >  1:1281567/1‑151 (MQ=255)                           ttttCACCGACAACAAAAAAGCTTTCGGTGATATTTTCCCTTCCGCCTTCAAGCGAGACGTTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCa                                                                                                                            <  1:827617/150‑1 (MQ=255)                            tttCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                          >  1:1180824/1‑151 (MQ=255)                              cccccaaggaaaaaaaaGCTTGGGGTTTTATTTTCCCTTCCTCCTTAAAGCGAGCCGTTATTCCCCGGCCCCTGTCCGCGGGCCAGGTTATAATTTTTTTTTCCAGCCAGGAAAAA‑C‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                                        <  2:658337/140‑1 (MQ=255)                                  cgacgaCAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                    >  1:787487/1‑151 (MQ=255)                                  cgacgaCAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                                    >  1:175814/1‑151 (MQ=255)                                    acgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGCATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                                  >  1:123733/1‑151 (MQ=255)                                    acaaaaaaaaaGCGTTGGGTTATATTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                                  <  1:473222/146‑1 (MQ=255)                                    acaaaaaaaaaGCGTTCGGTAATTTTTTCCCTTCCGACTTCAAGCGAGACGGTATCCCCCGGCCCCTGGTCGCGCGGAAGTTTATGAATTTTCTGTTCATGGGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCaa                                                                                                                        <  1:118117/140‑1 (MQ=255)                                       acAAAAAAGCGTTCGGTGATTTTTTCCCTTTCGACTTCAAGCGAGACGGTATTCCCGGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCGGGCGGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                                  <  1:332565/148‑1 (MQ=255)                                       aaaaaaaaGCGTTCGGTGATTTTTTCCTTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTTTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgc                                                                                                                  <  1:704966/146‑1 (MQ=255)                                         aaaaagGCGTTGGGGGAAATTTTCACTTCGGCCTCAAAGGGATACGTTATTCCCCGGCCCCGGGCCGGGGGGGAGGTTATAATTATCCGTTTCA‑GGAGGGAAAAGC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGgcg                                                                                                                 <  2:157070/147‑1 (MQ=255)                                         aaaaaaGCTTGGGGTGTTTTTTTCACTTCCCCCTCAAAGCGAGACGTTTTTCCCCGGCCCCTGGCCGCGGGGCAGGTTATAAATATTCTGTCCATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                              <  2:260485/150‑1 (MQ=255)                                         aaaaaaGCGTTGGGTGATATTTCCACTTTCGACTACAAGCGAGACGGTATCCCCCGCCCCTTGGTCGCGCGGCAGGTTAAAAATATTCTGTTCATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             <  1:751021/151‑1 (MQ=255)                                         aaaaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGGTCATGCAGGGAAAA‑A‑CTCCCCGCCCGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                             >  2:562053/1‑151 (MQ=255)                                           aaaaGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGcac                                                                                                           >  1:417201/1‑151 (MQ=255)                                              aGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGCGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACAcc                                                                                                        >  1:284313/1‑151 (MQ=255)                                                     ggTGATATTTTCACTTTCGACTTCAAGCGAGACGGTCTTACCCGGCCCCTGGTCGCGCGGCAGGTTATGATTATTCTGATCATGCAGGGAACA‑A‑CTCCCCGCCAGTGCAGCGATTTCACGCTCAGCCACATGACGCCGCACACCAGCTAAt                                                                                                 >  2:622077/1‑151 (MQ=255)                                                        gATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAAt                                                                                                 >  1:657565/1‑148 (MQ=255)                                                            ttttCACTTTGGTCTCAAGGGGGGACTGTATTACCCACCCTCGGGGCGCGGGGCAGGTTTGTAATTTTCGTTTCAGGCGGAGAAAA‑C‑CGCCCCGCCAGTGTAGCGACTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                            <  2:1029994/149‑1 (MQ=255)                                                                  cTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTACTAGAACTTCTCCAc                                                                                    >  2:562353/1‑151 (MQ=255)                                                                   tttCCACTAAAAGCTGGACGTTACTCCCCGCCCCCGGGCCGCGCGGCAGTTTATAAATATCTGGTTCATGAAGAAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAGCTtc                                                                                         <  2:946246/145‑1 (MQ=255)                                                                    tccGACTTCAGGCGAAGCGGTATTCCCCGGCCCCGGGTCGCGCGGCGGTTTATGAATTTTCTTTCCATGCAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCAc                                                                                    <  1:139096/147‑1 (MQ=255)                                                                     tCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGc                                                                                  >  2:130381/1‑150 (MQ=255)                                                                        aCTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATTAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                              >  1:501273/1‑151 (MQ=255)                                                                         cTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                             >  1:621428/1‑151 (MQ=255)                                                                          ttaaaGCGAGACGGTATTCCCGGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCGGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCt                                                                              <  1:616037/146‑1 (MQ=255)                                                                            cAAGGGACACGGTATGGCCCGCCCCCGGGCCGCGCGCCAGGTTAGAATTCTTTTTTTCATGCAGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGc                                                                            <  2:1177970/149‑1 (MQ=255)                                                                             aaGCGAGACGGTATTCCCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTTTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                                         <  1:348273/151‑1 (MQ=255)                                                                                 gagaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                                       >  2:1125765/1‑149 (MQ=255)                                                                                      gtttttccccGCCCCCGGGTCGCGCGGCAGTTTATAAATATTCTTTTCATGAAGAAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                 <  1:954551/143‑1 (MQ=255)                                                                                      gttATTCCCCGCCCCCGGGCCGCGCGACGGTTTATAAATTTTTTTTACAGGCGGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGt                                                                     <  1:178298/144‑1 (MQ=255)                                                                                      ggTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                 >  2:802538/1‑150 (MQ=255)                                                                                          ttACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCgcgc                                                             >  2:797773/1‑150 (MQ=255)                                                                                             cccGCCCCCTGGTCGCGCGGCGGTTTATAAATATTCTGTTCATGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGc                                                         <  1:1157053/151‑1 (MQ=255)                                                                                              ccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                        >  2:1059310/1‑151 (MQ=255)                                                                                               cGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                        >  2:260967/1‑150 (MQ=255)                                                                                                   cccTGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                   <  1:838409/151‑1 (MQ=255)                                                                                                       ggTCGCGCGGCAGGTTATGAATCTTCTGATCATGCAGGGAAAC‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACACTGTTCCCGCGCTCGCCAAGTGGCa                                                >  2:152168/1‑150 (MQ=255)                                                                                                       ggTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                                <  1:336941/150‑1 (MQ=255)                                                                                                         tcggtcggcggGTTATAAATATTCTGTGCTGACAGGTAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCACGCGCTCGCCAAGtt                                                   <  2:766608/140‑1 (MQ=255)                                                                                                         tCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAgg                                             >  1:259999/1‑151 (MQ=255)                                                                                                            cgcgGCAGTTTATAATTATTCTGTT‑AAGGAGGGAAAA‑ACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgt                                            <  2:437087/149‑1 (MQ=255)                                                                                                               ggCAGTTTATGAATTTTCTTTTCATGCAGGAAAAA‑C‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                        <  1:989687/150‑1 (MQ=255)                                                                                                               ggCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                        >  1:307701/1‑150 (MQ=255)                                                                                                                     ttAAGAATATTCTTTCAATGAAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                                     <  2:3713/147‑1 (MQ=255)                                                                                                                     gcaTGAATAATCTGTTCATCCAGGGAAAA‑A‑GTCCCCGCCAGTGTAGCGATTTCACGCTCACCAACATCGCGCGGCACCCCCGATAATAAAACTTCACGACGCATCTGACCAATGTTCCCGCCCTCGCCAGGTTGCAAGGTGTTAAACGAt                                  >  2:726735/3‑150 (MQ=255)                                                                                                                      tataAATATTCTGTTGAGGAAGGAAAAA‑C‑CTCCCCGCCAGTGTAGTGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATg                                 <  2:814224/150‑1 (MQ=255)                                                                                                                        tGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGcc                               >  1:1023119/1‑150 (MQ=255)                                                                                                                          aaTATTCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                            <  1:342680/151‑1 (MQ=255)                                                                                                                          aaTATTCTGTTCAGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                              <  1:552342/149‑1 (MQ=255)                                                                                                                          aaTATTCTGTTCAGGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                            <  1:763931/151‑1 (MQ=255)                                                                                                                              ttttGTTCATGCGGGAAAAA‑A‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc                         <  1:559692/147‑1 (MQ=255)                                                                                                                              ttCTTTTCGTGCAGGAAAAA‑C‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACg                            <  2:962189/147‑1 (MQ=255)                                                                                                                              ttCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  1:563941/1‑151 (MQ=255)                                                                                                                              ttCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTACCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  2:3213/1‑151 (MQ=255)                                                                                                                               tCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  1:149289/1‑150 (MQ=255)                                                                                                                               tCTGTTCATGCAGGAAAAA‑C‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc                         <  1:1122335/149‑1 (MQ=255)                                                                                                                                 tGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCg                        >  1:1120588/1‑148 (MQ=255)                                                                                                                                     cATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTc                 >  1:1198696/1‑151 (MQ=255)                                                                                                                                          aGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGtt             >  2:28298/1‑150 (MQ=255)                                                                                                                                          aGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt              >  2:2954/1‑149 (MQ=255)                                                                                                                                              aaaa‑g‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGTCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGGGTTAAACGATGCCACGGGCGAGACTTCCGTTGccac        <  2:1125717/151‑1 (MQ=255)                                                                                                                                               aaa‑a‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGccacc       >  1:595128/1‑151 (MQ=255)                                                                                                                                                   a‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCcaccac     >  1:556966/1‑150 (MQ=255)                                                                                                                                                      tCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACAcc  >  1:506154/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACC  >  NC_000913/3556060‑3556357

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