Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,200	0	G	A	40.3%	81.8 / 28.1	68	S217F (TCC→TTC)	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (34/6);  new base A (0/27);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 2.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTT  >  NC_000913/3556050‑3556346                                                                                                                                                       |                                                                                                                                                     aCCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAggg                                                                                                                                                      >  2:88796/1‑151 (MQ=255)           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:175814/1‑151 (MQ=255)                                            cgacgaCAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATg                                                                                                         >  1:787487/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)                                                                                                                                                       |                                                                                                                                                     ACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAA‑A‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTT  >  NC_000913/3556050‑3556346

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