Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	43.3%	106.1 / 13.4	62	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (27/7);  new base C (0/26);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 8.34e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620973‑1621268                                                                                                                                                  |                                                                                                                                                       tttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                                    >  2:13489/1‑150 (MQ=255)    cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                                <  1:905831/151‑1 (MQ=255)      gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                              >  1:640082/1‑151 (MQ=255)         gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcg                                                                                                                                           <  1:8159/151‑1 (MQ=255)             gACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACGCCCCATCTTTTCGCGAc                                                                                                                                       <  1:14975/151‑1 (MQ=255)             gACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                       >  2:1106275/1‑151 (MQ=255)              aCAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                       >  1:963491/1‑150 (MQ=255)                  aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                                  <  1:380120/151‑1 (MQ=255)                   aaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                                  >  1:392249/1‑150 (MQ=255)                       aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                             >  1:252017/1‑151 (MQ=255)                       aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACCACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                             >  2:671806/1‑151 (MQ=255)                             gaaaaaGCGACCAGAATAAAGCGCAAACCGGCAACCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTTATGACCACAAATTTTAGCCCCCACACGACCACTCCCAGTAGCGCCAACACCCCATCTTTCCGCGACATTCTACCGCCTCTg                                                                                                                        <  1:141746/148‑1 (MQ=255)                               aaaaGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACCCCCAATCTTTTCGCAACATTCTACCGCCTCTGa                                                                                                                       <  2:220511/149‑1 (MQ=255)                                 aaGCCACCAGCATAAAGCCCAAACCGGCCAGCATCAGCGGTGGCATGTTGTGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGCCCACTACGATTAGCCCCAACCCCCCGTCTTTGCGCGATATTCTCCCGCCTCTGAAttt                                                                                                                   <  1:465724/151‑1 (MQ=255)                                       ccAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtt                                                                                                             >  1:521013/1‑151 (MQ=255)                                       ccAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtt                                                                                                             >  1:360699/1‑151 (MQ=255)                                         aGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                           >  2:159768/1‑151 (MQ=255)                                         aGCAAAAAGCGCAACCCGGCCAGCATCAGCGGGGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCttt                                                                                                            <  2:553223/150‑1 (MQ=255)                                            aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                        <  2:932185/151‑1 (MQ=255)                                              aaaGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTaa                                                                                                       <  2:381467/150‑1 (MQ=255)                                                 gcgcAACCCGGCCAGCATCAGCGGTGGCATGTTATGAGGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                    <  2:351745/150‑1 (MQ=255)                                                  cgcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                                  >  1:854410/1‑151 (MQ=255)                                                   gcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                                  >  1:926364/1‑150 (MQ=255)                                                     aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                                <  2:871895/147‑1 (MQ=255)                                                     aaacgagCCACCAACAGGGGTGGCTTGTTATGAAGCCCCTCTTTGATGACCACAAAAGTTAGCCCCTATACGCCCGCTCCCGTTAGCCCCAACCCCCCTTCTTTTCGGGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                               <  1:1116642/145‑1 (MQ=255)                                                        ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                            >  2:705496/1‑151 (MQ=255)                                                        ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                            >  2:130206/1‑151 (MQ=255)                                                              aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                       >  2:100938/1‑150 (MQ=255)                                                               gCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGa                                                                                      >  2:118892/1‑150 (MQ=255)                                                                cATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                    >  2:289420/1‑151 (MQ=255)                                                                 aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGCCCACAAAATTTACCCCCCATACGCCCCCTCCCATTAGCTCCAACCCCCCATCTTTTCGCACCATTCTCCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                    <  1:990272/150‑1 (MQ=255)                                                                 aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGCCCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                                   <  2:970778/151‑1 (MQ=255)                                                                       ggTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                             >  1:569712/1‑151 (MQ=255)                                                                         tGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                            <  2:755679/150‑1 (MQ=255)                                                                          ggCATGTTATGAAGCCCCTCTTTGATGACCACAAAATTTACCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                          <  2:786489/151‑1 (MQ=255)                                                                                 tATGAAGCCCCATTTGTATGACCAAAAAATTTCCCCCCCAAACCACCCCTACCAGTCGCGCCAACTCCCAATTTTTTCTCGATATGTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACa                                                                    <  1:137295/150‑1 (MQ=255)                                                                                  aTGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGt                                                                  >  1:465442/1‑151 (MQ=255)                                                                                     aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGa                                                               <  2:154947/151‑1 (MQ=255)                                                                                     aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                                >  2:310756/1‑150 (MQ=255)                                                                                        ccccACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                            <  2:683326/151‑1 (MQ=255)                                                                                         cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  2:387236/151‑1 (MQ=255)                                                                                               ttGATGCCCACAAAATTGAGCCCCCATACGCCCACTACCAGTCGCGCCAACCCCCCATCTTTGCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGt                                                      <  2:26877/150‑1 (MQ=255)                                                                                                  atgaCCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCg                                                    <  2:859530/149‑1 (MQ=255)                                                                                                      ccccaAAATTAAGCCCCCATACGCCCCCCCCCAGAACCCCCAACACCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                               <  1:22274/147‑1 (MQ=255)                                                                                                      ccACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGtt                                              >  2:563100/1‑151 (MQ=255)                                                                                                          aaaaTTTAGCCCCCATACACCCACTCCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacga                                          <  1:539989/151‑1 (MQ=255)                                                                                                              tttAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTCTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                      >  1:993898/1‑151 (MQ=255)                                                                                                                 aGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   >  1:442119/1‑151 (MQ=255)                                                                                                                   cccccATCCGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   <  2:628747/149‑1 (MQ=255)                                                                                                                    ccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                                >  1:1012228/1‑151 (MQ=255)                                                                                                                       cATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                             >  1:763724/1‑151 (MQ=255)                                                                                                                         tccGACCCCTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                            <  2:679265/148‑1 (MQ=255)                                                                                                                           cgcccACTCCCAGTAGCGCCAACCCCCCATTTTTCCGCGACATTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                           <  2:1034466/146‑1 (MQ=255)                                                                                                                              ccACTACCAGTAGCCCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                      <  2:496264/151‑1 (MQ=255)                                                                                                                                 ctcccgttggcgcGAACCCCCCATCTTTTCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAGTCGTTAGTCGCCGGTTACGGCGGCATGAATGCGCaa                          <  1:1115145/135‑1 (MQ=255)                                                                                                                                    ccAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt                <  2:369396/151‑1 (MQ=255)                                                                                                                                    caaGTAGCCCCAACCCCCCACCTTTCCGTGATATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt                <  2:653251/149‑1 (MQ=255)                                                                                                                                     cAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCta               <  1:118744/151‑1 (MQ=255)                                                                                                                                      aGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCtat              >  2:733096/1‑151 (MQ=255)                                                                                                                                         aGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATAc            <  2:293895/150‑1 (MQ=255)                                                                                                                                            gcCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg        >  2:596181/1‑151 (MQ=255)                                                                                                                                               aaCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACtt          <  1:996835/146‑1 (MQ=255)                                                                                                                                                 cacCCCATCTTTTCGGGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAttt   <  1:37289/151‑1 (MQ=255)                                                                                                                                                  acCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtttt  >  2:467930/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                       TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620973‑1621268

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