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.7%	94.4 / 25.6	72	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (31/9);  new base C (0/31);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 2.51e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTC  >  NC_000913/1620973‑1621261                                                                                                                                                  |                                                                                                                                                attCGGTCGTGCGACAAAAAAGATAGCCGGAAACGCGACCAGCATAAAGCTCAAACCGCCCAGCATCAGCGGCGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACCACCAGTATCGCCAACAcccc                                                                                                                                             >  1:1285654/2‑150 (MQ=255)     ggTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCt                                                                                                                                         <  1:543888/150‑1 (MQ=255)        cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                     <  1:650355/151‑1 (MQ=255)             gACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                <  2:582483/151‑1 (MQ=255)              aCAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGAc                                                                                                                                <  1:470473/150‑1 (MQ=255)                aaaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAtt                                                                                                                             >  2:1135607/1‑151 (MQ=255)                 aaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTc                                                                                                                            >  2:392162/1‑151 (MQ=255)                  aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                           >  1:822876/1‑151 (MQ=255)                  aaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCt                                                                                                                           >  2:144065/1‑151 (MQ=255)                    aaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGGGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTAc                                                                                                                         >  1:1193372/1‑151 (MQ=255)                      gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCg                                                                                                                       >  2:999545/1‑151 (MQ=255)                       aTAGCCGGAAAAGAGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCCTGTTATGAAGCCCCACTTTGATGACCCCAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                      >  1:181477/1‑151 (MQ=255)                          gCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                    <  1:833864/150‑1 (MQ=255)                               aaaaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAt                                                                                                              >  2:854224/1‑151 (MQ=255)                                 aaGCGACCAGCAAAAACCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCTTTTGATGCCCACAAAATTTAGCCCCCATACGACCCCTACCAGTAGCGCCAACACCCCATCTTTCCGCGACTTTCTACCGCCTCTGAAttt                                                                                                            <  1:293449/151‑1 (MQ=255)                                     ggccAGAAAAAAGCCAAAACCGGCCAGCACCAGGGGTGGCATGTTATGAACCCCCACTTAGATAACCACAAAATTTGGCCCCCATACGCCCTCTACCAGTAGCCCCAACCCCCCACTTTTCCGCGAAATTCTACCGCCTCTGAATTTc                                                                                                           <  1:402622/146‑1 (MQ=255)                                     gACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                        >  2:195523/1‑151 (MQ=255)                                            aaaaaGCCCAAGCCGCCCACCACCGGCGTTGGCATTTTATAAACCCCCATTTTAATGACCACAAAATTTTCCCCCCATCCGACCACTCCCATTAGCGCCACCCCCCCTCCTTTTCGCGATTTTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                  <  1:372990/148‑1 (MQ=255)                                            aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                  >  2:1052236/1‑150 (MQ=255)                                             tAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                 >  1:11880/1‑150 (MQ=255)                                               aaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                              >  1:1090270/1‑151 (MQ=255)                                                   gcAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                          >  1:966197/1‑151 (MQ=255)                                                     aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCaa                                                                                            <  2:1295141/147‑1 (MQ=255)                                                       aCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                      >  2:314122/1‑151 (MQ=255)                                                        ccGGCCAGCACCGGGGGGGGCATGTTAGGAAACCCCATTTTGATGACCAAAAATTTAAGCCCCAACACGCCCCCTCCCAGTCGCCCCAACCCCCCATCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                     <  2:95017/151‑1 (MQ=255)                                                          ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                     >  2:788596/1‑149 (MQ=255)                                                          ggCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                     >  1:1058092/1‑149 (MQ=255)                                                            ccAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAATTTTACCCCCCATACGCCCCCTACCAGTAGCGCCAACCCCCCATTTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                  <  2:582437/150‑1 (MQ=255)                                                            ccAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCct                                                                                                                    <  2:1078968/116‑1 (MQ=255)                                                                cATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                             >  2:761671/1‑151 (MQ=255)                                                                   cAGCGGTGGCATGTTATGAACCCCCCCTTTGATGACCACAAAATTTAGCCCCCATCCGCCCACCCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                          <  2:428796/151‑1 (MQ=255)                                                                    aGCGGTGTGATGTTATGAGCCCCCCCTTTAAACACCAAAAATTTTACCCCCAATACAACCCCCCCCATTAGCCCCACCGCCCAATCTTTTCGAGAAATTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                         <  1:416604/151‑1 (MQ=255)                                                                     ggggTGGCATGTTATGAACCCCCCCTTTGATGACCAAAAAATTTCCCCCCCATACGACCCCTCCCATTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                        <  1:1086063/149‑1 (MQ=255)                                                                      cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                       >  2:714952/1‑151 (MQ=255)                                                                       ggTGCCATTTTATGAACCCCCACTTTGATCACCACAAATTTTGCCCCCCATCCGCCCACTCCCAGTAGCGCAACCCCCCCATTTTTCCGCAACTTTCAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                            <  1:1139408/145‑1 (MQ=255)                                                                          ggaatattagaaaGCCCCCTTTTGATCACCAAAAATTTTAGCCCCCATATAACCACTCCCAACAGCGCAACCACCCTTTTTTTTCCCTAAATTCTTCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                             <  1:487556/130‑1 (MQ=255)                                                                          ggCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCGATCAACTTAGCTGAATTTACtttt                                                                     <  2:969613/149‑1 (MQ=255)                                                                             aTGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                 >  1:299505/1‑150 (MQ=255)                                                                                ttATGAAGCCCTACTTTGATCACCAAAAAATTTAGCCCCAATACAACCACTCCCAGTAGCGCCAACCCCCCTTCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                 <  2:756227/147‑1 (MQ=255)                                                                                 tATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAAc                                                              >  2:286885/1‑149 (MQ=255)                                                                                  aTGAAGCCCCACTTTTATGACCACAAAATTTAGCCCCCATACGACCTCTACCAGTAGCCCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTTACa                                                             <  2:196916/149‑1 (MQ=255)                                                                                  aTGAACCCCCACTTGGTTGACCAAAAAATTAAGCCCCCATACGCCCCCTCCCAGAGGCCCCACCCCCCCATTTTTTCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACa                                                             <  2:196918/149‑1 (MQ=255)                                                                                   tGAACCCCCCTTTTGATGCCCAAAAAATTTACCCCCCATCCGCCCACTCCCAGTAGCGCAACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAg                                                            <  2:1143163/149‑1 (MQ=255)                                                                                     aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                         >  1:812703/1‑150 (MQ=255)                                                                                       gCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                      >  1:1173520/1‑151 (MQ=255)                                                                                        ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAt                                                       <  2:782444/149‑1 (MQ=255)                                                                                        ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                     <  2:86504/151‑1 (MQ=255)                                                                                        ccccACTTTGATCACCAAAAAATTTACCCCCCATCCGCCCCCTACCAGTAGCCCCAACCCCCCATCTTTCCGCGATTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAtt                                                      <  1:779796/150‑1 (MQ=255)                                                                                        ccccACTTTGATACCAACAAAATTTACCCCCCATACGACCCCTCCCAGTAGCACCACCCCCCCTTTTTTCCGCGATATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAATAGTTGAtt                                                      <  1:285025/150‑1 (MQ=255)                                                                                         cccACTTTGATCACCACAAATTTTAGCCCCCATCCGCCCAGTCCCAGTAGCCCAACCCCCCCATCTTTCCCCCACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                    <  1:1059596/151‑1 (MQ=255)                                                                                            aCTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                  >  1:936597/1‑150 (MQ=255)                                                                                               ttGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCACCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAg                                                <  2:81141/149‑1 (MQ=255)                                                                                                  atgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                           >  2:558994/1‑151 (MQ=255)                                                                                                        acaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                     >  2:425091/1‑151 (MQ=255)                                                                                                        aaaaaaTTTAGCCCCCATCCACCCCCTCCAGGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                     <  1:254793/149‑1 (MQ=255)                                                                                                          aaatTTTACCCCCCATACGCCCACTACCAGTAGCGCAAACCCCCATTATTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                    <  1:231815/150‑1 (MQ=255)                                                                                                            aaTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                 >  2:123426/1‑151 (MQ=255)                                                                                                             aTTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                >  1:597606/1‑151 (MQ=255)                                                                                                              tttAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                <  2:877841/150‑1 (MQ=255)                                                                                                                   cccccATACGCCCCCTCCCAGTACCCCCAACCCCCCTTCTTTCCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGc                               <  2:660087/146‑1 (MQ=255)                                                                                                                     cccATACGCCCACTCCCAGAGGCCCAACCCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATg                        <  1:386431/151‑1 (MQ=255)                                                                                                                           cgccctctcccAGTAGCCCCAACCCCCCATCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                   <  2:454915/141‑1 (MQ=255)                                                                                                                           ccaCCACTCCCAGTAGCGCCAACCCCCAATTTTTTCGCGACATTCTTCCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaaa                  <  2:544766/149‑1 (MQ=255)                                                                                                                           cGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                   >  1:795556/1‑150 (MQ=255)                                                                                                                              cccctcccAGTAGCGCCAACACCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaaa                  <  2:963891/142‑1 (MQ=255)                                                                                                                              ccACTACCAGTAGCCCCACCCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                    <  2:709333/146‑1 (MQ=255)                                                                                                                                  tACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc            >  1:85012/1‑150 (MQ=255)                                                                                                                                  tACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg           >  1:153394/1‑151 (MQ=255)                                                                                                                                    ccAGTAGCGCCAACCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg           <  1:251082/149‑1 (MQ=255)                                                                                                                                    ccAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCg           <  2:487694/149‑1 (MQ=255)                                                                                                                                     cAGTAGCGCCACCCCCCCATCTTTCCGCCACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc            <  2:111687/147‑1 (MQ=255)                                                                                                                                            gcCAACCCCCCTTTTTTCCCCGACGTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTc  <  1:1250751/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                TTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTC  >  NC_000913/1620973‑1621261

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