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

TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620973‑1621268                                                                                                                                                  |                                                                                                                                                        tttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcccc                                                                                                                                                     >  2:588553/1‑150 (MQ=255)  ttCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑At                                                                                                                                                  >  1:215469/1‑151 (MQ=255)   tCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATc                                                                                                                                                 >  2:527327/1‑151 (MQ=255)    cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCt                                                                                                                                                <  1:1066170/151‑1 (MQ=255)     ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCtt                                                                                                                                               <  1:537946/151‑1 (MQ=255)      gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCttt                                                                                                                                              >  2:384580/1‑151 (MQ=255)        cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTc                                                                                                                                            >  2:171159/1‑151 (MQ=255)                 aaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTc                                                                                                                                   <  2:783641/151‑1 (MQ=255)                  aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCt                                                                                                                                  >  2:1032403/1‑151 (MQ=255)                     aGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTAcc                                                                                                                               <  1:773049/151‑1 (MQ=255)                      gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCg                                                                                                                              >  1:960592/1‑151 (MQ=255)                          gCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCctc                                                                                                                          <  1:891734/151‑1 (MQ=255)                            cGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                        >  1:391453/1‑151 (MQ=255)                             ggAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                        >  1:435254/1‑150 (MQ=255)                                 aaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                   <  1:778782/151‑1 (MQ=255)                                 aaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAAttt                                                                                                                   <  1:150901/151‑1 (MQ=255)                                      aCCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                                <  1:62539/149‑1 (MQ=255)                                       ccAGCAAAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                               <  2:1037932/149‑1 (MQ=255)                                           cATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                         >  1:369457/1‑151 (MQ=255)                                            aaaaaGCGCAAACCGGCCAGCATCAGGGGGGGCATGTTATGAAGCCCCACTTGGATGACCACAAAATTTGGCCCCCATACGCCCACTACCAGTGGCGCCACCCCCCC‑ATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                           <  1:408610/146‑1 (MQ=255)                                            aTAAACCGAAACCCGCCAAGCATCAGCGGTGCCATGTTATAAAGCCCCACTTTGATGACCAAAAAATTTACCCCCCATACGACCCCTACCAGTAGCCCCAACCCCCC‑ATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCtttt                                                                                                           <  2:499310/148‑1 (MQ=255)                                              aaaCCGCAACCCGCCTACCACCAGCGGTGCCATGTTATGAACCCCCACTTTGATGACCACAAAATTTACCCCCAATACTCCCCCTACAAGTAGCCCCAACCCCCC‑ATCTTTGCGCCACTTTTTACCGCCTCTGAATTTCATCttt                                                                                                            <  2:379738/145‑1 (MQ=255)                                                     aaaCCGGCCAGCATCAGCGGGGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                                <  2:974684/150‑1 (MQ=255)                                                        cgggCAAGCAACGGCGGTGGCATGTTATGAACCCCCACTTTAATGACCAAAAAATTGCCCCCCCATCCGCCCCCT‑GCAGGTGCCCCAACCCCCCAATCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                                <  1:363966/145‑1 (MQ=255)                                                           gCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                         >  2:1026510/1‑151 (MQ=255)                                                            ccAGCATAAGCGGTGGCATTTTATGAAGCCCCACTTTGATAACCACAAAATTTACCCCCCATACGACCCCTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                         <  2:663803/150‑1 (MQ=255)                                                                 aTCAGCGTGGGCATTTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACGACCACTACCAGTAGCGCCAACCCCCC‑ATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                                   <  1:527072/151‑1 (MQ=255)                                                                 aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                                   >  1:777090/1‑151 (MQ=255)                                                                      cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                                 <  1:1024360/148‑1 (MQ=255)                                                                       ggTGGCATTTTATGAACCCCCACTTTGATGA‑CAAAAACTTTACCCCCCATACGACCTTTCCCAGTGGCGCCAACCCCCA‑ATCTTTTCGCGACATTCGACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                                 <  1:178165/146‑1 (MQ=255)                                                                       ggTGCCCTGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACACCCACCCCCAGTAGCCCCAACCCCCC‑ACCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                             <  2:810724/151‑1 (MQ=255)                                                                          ggCATGTTTTGAAGCCCACCTTTGATTACCACAAAATTTAGCCCCCATCCGCCCACTCCCAGTAGCCCCACCCCCCC‑ATCTTTGCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                              <  2:649138/147‑1 (MQ=255)                                                                            cATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                        >  1:957695/1‑151 (MQ=255)                                                                            cATGTTATGAAGCCCCACTTTGAAGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                        >  2:1034145/1‑151 (MQ=255)                                                                              tGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                        <  1:521581/149‑1 (MQ=255)                                                                              tGTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCTTACGACCCCTCCCAGTAGCCCCAACCCCCC‑ACTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                                         <  2:947788/148‑1 (MQ=255)                                                                              tGTTATGAACCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCtt                                                                         <  2:215298/148‑1 (MQ=255)                                                                               gTTATGAAGCCCCACTTTGATGACCACAAATTTTACCCCCCATACGCCCCCTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTaa                                                                      <  2:423281/150‑1 (MQ=255)                                                                                ttATGATCCCCCACTTTGATGCCCACAAAATTTCCCCCCCATACGCCCCCTCCCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTaa                                                                      <  2:510753/149‑1 (MQ=255)                                                                                ttATGAAGCCCCCCTTTAATGACCCCAAAATTTACCCCCCATACGACCACTCCCAGTGGCGCCAACACCCC‑ATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                     <  1:79229/150‑1 (MQ=255)                                                                                 tATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                                     >  2:195509/1‑149 (MQ=255)                                                                                    gAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                                >  1:508094/1‑151 (MQ=255)                                                                                     aacccccACTTTGATGACCACAAATTTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                                <  2:76337/147‑1 (MQ=255)                                                                                        ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                            >  1:877154/1‑151 (MQ=255)                                                                                         ccccCTTTGAGGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                           <  1:52934/151‑1 (MQ=255)                                                                                            aCTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                         >  2:994702/1‑150 (MQ=255)                                                                                                 gatgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGc                                                   <  1:384787/151‑1 (MQ=255)                                                                                                   tgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCg                                                 >  2:80001/1‑151 (MQ=255)                                                                                                      ccACAAAATTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                               <  1:371142/150‑1 (MQ=255)                                                                                                          aaatTTAAGCCCCAATACGACCACTACCAGTAGCCCAAACCCCCC‑ATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                            <  1:738364/149‑1 (MQ=255)                                                                                                           aaaTTTGCCCCCCATCCGCCCCCTCCCAGTAGCCCCAACACCCC‑TTTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                         <  1:490173/151‑1 (MQ=255)                                                                                                           aaaTTTACCCCCCATCCGCCCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                         <  2:645212/151‑1 (MQ=255)                                                                                                            aaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                        >  1:280473/1‑151 (MQ=255)                                                                                                              tttACCCCCCATCCGCCCCCTCCCAGTAGCCCCAACACCCC‑ATCTTTTCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                        <  2:50211/149‑1 (MQ=255)                                                                                                               ttACCCCCCAAACGCCCACTACAAGTAGCACCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                                      <  1:1037542/150‑1 (MQ=255)                                                                                                                 aGCCCCCATACGACCACTCCCAGTGGCCCCAACCCCCC‑ATCTTTTCGCGACATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCa                                     <  1:42352/149‑1 (MQ=255)                                                                                                                  gCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   >  2:576526/1‑150 (MQ=255)                                                                                                                   cccccATACACCCACTCCCATTAGCGCCAACACCCC‑ATTTTTTCCCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                   <  2:230447/149‑1 (MQ=255)                                                                                                                    ccccATACGACCACTCCCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                                <  1:274332/151‑1 (MQ=255)                                                                                                                    ccccATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                                >  2:241818/1‑151 (MQ=255)                                                                                                                      ccATCCCACCCCTACCATTAGCCCCAACCCCCC‑ATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                              <  2:150100/151‑1 (MQ=255)                                                                                                                      ccATACGACCACTACCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAGTTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                              <  1:1043162/151‑1 (MQ=255)                                                                                                                        aTACGACCACTCCCAGTAGCGCCAACCCCCC‑ATCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                             <  1:508872/150‑1 (MQ=255)                                                                                                                        aTACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                            >  1:53906/1‑151 (MQ=255)                                                                                                                          aCGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                           >  2:321841/1‑150 (MQ=255)                                                                                                                              ccACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTGCGACGGCATTAATGCGCAAATa                       <  1:484095/150‑1 (MQ=255)                                                                                                                               cACTCCCAGTAGCGCCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTATAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAg                     <  1:919778/151‑1 (MQ=255)                                                                                                                                  tACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACTGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg                  >  2:1011058/1‑151 (MQ=255)                                                                                                                                           cgcCAACCCCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc         <  1:956769/151‑1 (MQ=255)                                                                                                                                           cgcCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc         >  2:397664/1‑151 (MQ=255)                                                                                                                                                  acCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAtttt  >  2:437916/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                        TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCC‑ATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTTT  >  NC_000913/1620973‑1621268

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