Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	40.7%	161.1 / 25.2	91	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (45/9);  new base C (0/37);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 3.10e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAT  >  NC_000913/1620975‑1621265                                                                                                                                                |                                                                                                                                                    tCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATc                                                                                                                                              >  2:1083314/1‑151 (MQ=255)  cGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCAtat                                                                                                                                             <  2:751439/151‑3 (MQ=255)      cGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTc                                                                                                                                         >  2:495876/1‑151 (MQ=255)       gTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcg                                                                                                                                        <  1:633269/151‑1 (MQ=255)         gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgcg                                                                                                                                      <  1:1315975/151‑1 (MQ=255)              aaaaaaGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAtt                                                                                                                                 <  1:229592/151‑1 (MQ=255)                   aGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTAcc                                                                                                                            >  2:1363320/1‑151 (MQ=255)                    gATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTAcc                                                                                                                            >  2:206702/1‑150 (MQ=255)                     aTAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGc                                                                                                                          <  1:937266/151‑1 (MQ=255)                      tAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGcc                                                                                                                         >  1:931758/1‑151 (MQ=255)                          cGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTg                                                                                                                     >  1:64150/1‑151 (MQ=255)                            gAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGa                                                                                                                    >  1:759046/1‑150 (MQ=255)                             aaaaGCGACCAGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATAAAGCCCCCCTTTGATGCCCACAAAATTAACCCCCCATCCGACCACTCCCAGTAGCCCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAAt                                                                                                                  <  1:1146946/151‑1 (MQ=255)                              aaaGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAAtt                                                                                                                 >  1:658885/1‑151 (MQ=255)                                  cGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                             >  1:530702/1‑151 (MQ=255)                                   gACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                            >  1:496932/1‑151 (MQ=255)                                    aCCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                           >  2:1059772/1‑151 (MQ=255)                                    aCCAGAATAAAGCGCAAACCGGCCAGCACCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCACTACCATTAGCCCCACCCCCCCATTTTTCCGCGACATTCTACCGCCTCTGAATTTCATc                                                                                                            <  1:760149/150‑1 (MQ=255)                                     ccAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtt                                                                                                          >  1:1037055/1‑151 (MQ=255)                                      cAGCAAAAAGCGAAACCCGGCCAGCATGAGCGGTGGCATTTTATGAAGCCCCTTTTTGATGACCACAAAATTTGCCCCCCATACGCCCCCTCCAAGTAGCCCCAACCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTc                                                                                                               <  2:926028/145‑1 (MQ=255)                                          aaaaaGCCAAACCCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCACTTTGATGCCCAAAAAATTTAGCCCCCAACCGCCCCCTCCCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCtt                                                                                                          <  2:796987/144‑1 (MQ=255)                                          aTAAAGCGCAACCCGCCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTGATGACCAAAAAATTTACCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                     <  1:1029844/151‑1 (MQ=255)                                          aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCATTTTGATAACCACAAATTTTACCCCCCACCCTCCCACTCCCAGTAGCGCCAACCCCCCATTTTTTCGCGAAATTCTACCGCCTCTGAATTTCATCtttt                                                                                                        <  1:181674/148‑1 (MQ=255)                                            aaaGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAACCCCCACTTTGATAACCACAAAATTTAGCCCCCATCCGCCCACTACCAGTAGCCCCACCCCCCCATTTTTTCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTaa                                                                                                    <  1:464332/150‑1 (MQ=255)                                            aaaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTa                                                                                                     >  2:130889/1‑149 (MQ=255)                                             aaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                  >  1:569743/1‑151 (MQ=255)                                                 gcaccccGCCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTACCCCCCTTACCCCCACTACCAGAAGCGCCAACACCCCATTTTTGCCCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCa                                                                                                 <  2:548781/143‑1 (MQ=255)                                                   aacccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                            <  1:138436/148‑1 (MQ=255)                                                   aaaCGGCCCCGCATGAGTGGTGGAATCTTAGGAAGCCCATTTTTGAGAACCAAAAATTTTTGCCCCCATACGACCCCCCCAATTAGCGCAAGCACCCCACCTTTTCGGGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                             <  1:1298345/150‑1 (MQ=255)                                                   aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                             >  2:1061245/1‑150 (MQ=255)                                                   aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCa                                                                                             >  2:826024/1‑150 (MQ=255)                                                   aaaCCGGCCAGCATCAGCGGTGGCATGATAAGAAGACCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACAGCCTCTGAATTTCATCTTTTGTAAGCAATc                                                                                              >  2:411872/1‑149 (MQ=255)                                                      cggcccAGCATCAGGGGTGGCATGTTATGAACCCCCACTTTGACAACCACAAAATTTACCCCCCATACGCCCACTACAAAAAGCCCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                          <  2:710496/146‑1 (MQ=255)                                                      ccGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATACGACCACTCCCAGTAGCCCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                          <  1:361126/150‑1 (MQ=255)                                                      ccGGCCAGCATCAGCGGTGGCATGTTATAAAGCCCCATTTTGATGACCACAAAATTTACCCCCCATACGCCCCCTACCAGTAGCCCCAACCCCCCATCTTTTCGCGAAATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                          <  1:987732/150‑1 (MQ=255)                                                          ccAGCACCAGCGGTGGCATTTTATGAAGCCCCGCTTTGATGACCAAAAAATTTACCCCCCATACCCCCACTCCCAGTAGCGCCAACCCCCCCTCTTTCCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGc                                                                                      <  2:1040769/150‑1 (MQ=255)                                                           cAGCACCGGCGGTGGCTGGTTATGACCCCCCCTTTTGATGACCCCAAAATTTACCCCCAATACGCCCACCCCCAGTAGCGCCAACCCCCTATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACtt                                                                                         <  2:1354760/146‑1 (MQ=255)                                                            aGCGTCAGCGGTAGCATGTTACGAAACCCCTCTTTATTAACCAAAAAATTTAGCCCCCACTCCCCCACTCCTAGGAGCGCCAACCCCCCATTTATAGGGGACATTAAACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                    <  2:1207404/150‑1 (MQ=255)                                                               aTCAGCGGGGGCATGTTATAAAGCCCCATTTTGATGCCAAAAAAATTTACCCCCCATACGCCCACTCCCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                 <  1:936902/150‑1 (MQ=255)                                                                  aGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                             >  1:863522/1‑151 (MQ=255)                                                                   ggggTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTAGCCCCCATCCACCCACTCCCAGTAGCGCCAACCCCCCATTTTTTCGCGACATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTAc                                                                             <  1:1083366/148‑1 (MQ=255)                                                                    cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                           >  1:350331/1‑151 (MQ=255)                                                                    cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                           >  1:926318/1‑151 (MQ=255)                                                                     tgTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCCCCAACACCCCATCTTTTCGCGACATTCCACCGCCTCTGCATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                          >  2:639836/2‑151 (MQ=255)                                                                      gTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                         >  2:784045/1‑151 (MQ=255)                                                                        ggaattttATGACGCCCCATTTTAATGCCCAAAAATTTTAGCCCCCATCCGCCCACTACCAGTAGCCCCAACCCCCCATTTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCt                                                                       <  2:1331446/145‑1 (MQ=255)                                                                        ggCATTTTATGAACCCCCCCTTTGATGACCACAAAATTTACCCCCCATCCGCCCACTACCAGTAGCCCCAACCCCCCTTCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACttt                                                                          <  2:935038/148‑1 (MQ=255)                                                                        ggCATGTTATGAAGCCCCTCTTTGATGCCCACAAAATTTAGCCCCCACACGACCCCTCCCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTc                                                                        <  1:558266/150‑1 (MQ=255)                                                                                   aaccccccttttGATGACCAAAAAATTTACCCCCCATACGACCCCTCCCAGTAGCGCCAACCCCCCATTTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGa                                                            <  2:332414/142‑1 (MQ=255)                                                                                   aacccccATTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATTTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                              <  1:1428393/146‑1 (MQ=255)                                                                                   aaGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTg                                                             >  1:430301/1‑150 (MQ=255)                                                                                    aGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGAt                                                           >  2:60291/1‑151 (MQ=255)                                                                                      ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                         <  1:439268/151‑1 (MQ=255)                                                                                      ccccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTc                                                         >  1:1446310/1‑151 (MQ=255)                                                                                       cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                        >  1:769463/1‑151 (MQ=255)                                                                                       cccACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                        >  1:887865/1‑151 (MQ=255)                                                                                       cccACTTTGATGACCAAAAAATTTACCCCCCATACGACCACTACCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCg                                                        <  1:1341900/151‑1 (MQ=255)                                                                                          aCTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGtt                                                      >  1:746256/1‑150 (MQ=255)                                                                                                atgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                               >  2:1178310/1‑151 (MQ=255)                                                                                                  gaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCgg                                             >  2:187375/1‑151 (MQ=255)                                                                                                   aCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGt                                            >  2:300035/1‑151 (MQ=255)                                                                                                      acaAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTac                                         >  1:165424/1‑151 (MQ=255)                                                                                                       caAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                        >  1:669764/1‑151 (MQ=255)                                                                                                       caAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                        >  1:1448203/1‑151 (MQ=255)                                                                                                        aaaaTTAAGCCCCCATACGCCCCCTCCCATTAGCGCCACCCCCCCTTCTTTCCCCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacg                                        <  1:844386/150‑1 (MQ=255)                                                                                                         aaaTTTGCCCCCCATCCGCCCACTCCCAGTAGCGCCACCCCCCCATTTTTTCGCAAATTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                      <  1:401587/151‑1 (MQ=255)                                                                                                            tttACCCCCCATACGCCCCCTCCCAGTGGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACgg                                    <  2:178227/150‑1 (MQ=255)                                                                                                              tAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAt                                 >  1:294552/1‑151 (MQ=255)                                                                                                                 cccccaaccccccACTACCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                <  1:1441017/138‑1 (MQ=255)                                                                                                                 cccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTa                               >  2:661038/1‑150 (MQ=255)                                                                                                                    ccATACGACAACTACCAGTGGCGCCAACACCCCATCTTTTCGCGACATTCTACTGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGCCGGCATTaa                              >  2:1057956/1‑148 (MQ=255)                                                                                                                       tACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                         >  1:24184/1‑150 (MQ=255)                                                                                                                         cGACCACTCCCAGTAGCCCCACCCCCCCATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCa                        <  2:179687/149‑1 (MQ=255)                                                                                                                         cGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCaa                       >  1:1434642/1‑150 (MQ=255)                                                                                                                            ccccccccAGTAGCGCAAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                    <  1:1043943/144‑1 (MQ=255)                                                                                                                            cccccacaaGTAGCACCACCCCCCCATCTTTGCGCGTCATTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                    <  2:1459510/142‑1 (MQ=255)                                                                                                                            ccccTACCAGTGGCCCCAACACCCCACCTTTCCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                    <  1:536680/147‑1 (MQ=255)                                                                                                                            ccACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                   <  1:734866/151‑1 (MQ=255)                                                                                                                            ccACTCCCAGTAGCGCCAACCCCCAATCTTTCCGCGACTTTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                   <  1:1048306/151‑1 (MQ=255)                                                                                                                            ccACTCCCAGTAGCCCCACCCCCCCATCTTTTCGCGACTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                    <  2:703302/150‑1 (MQ=255)                                                                                                                            ccACTACCAGTAGCCCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATaa                   <  2:249852/151‑1 (MQ=255)                                                                                                                             cACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATa                    >  1:1256494/1‑149 (MQ=255)                                                                                                                               cTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc                >  1:534119/1‑151 (MQ=255)                                                                                                                                 aCCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTc                >  2:1257853/1‑149 (MQ=255)                                                                                                                                  ccggaggcgcCACCCCCCCTTTTTTTCGCGACATCCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGc              <  1:328990/144‑1 (MQ=255)                                                                                                                                  ccAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt                 <  1:661052/147‑1 (MQ=255)                                                                                                                                           ccAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCg     <  1:201943/150‑1 (MQ=255)                                                                                                                                           ccAACCCCCCATCTTTCCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg    <  1:71154/151‑1 (MQ=255)                                                                                                                                            cAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCgg    >  2:140530/1‑150 (MQ=255)                                                                                                                                             aaCACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAt  >  2:447966/1‑151 (MQ=255)                                                                                                                                             aaCACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAt  >  1:576812/1‑151 (MQ=255)                                                                                                                                                |                                                                                                                                                    TCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAT  >  NC_000913/1620975‑1621265

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: ‑