Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,610	0	A	C	40.5%	128.1 / 19.9	75	W186G (TGG→GGG)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (32/12);  new base C (0/30);  total (32/43)
Fisher's exact test for biased strand distribution p-value = 2.86e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATC  >  NC_000913/1174463‑1174755                                                                                                                                                      |                                                                                                                                                  aGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCCCCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCcca                                                                                                                                                   <  2:299865/148‑1 (MQ=255)  ggCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCcaccac                                                                                                                                               <  1:382666/151‑1 (MQ=255)  gcccccgaccgCCAGACCGGAAGGAAAAGAGCCCACAAGCAGAGGGGGGCG‑AACGGCCCCCCCTTAACACCGGTTAATTTTGGTTGC‑CCAAACCATTGGCCAAAGCTGGAAATAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACg                                                                                                                                              <  1:580589/147‑1 (MQ=255)   gCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACg                                                                                                                                              >  2:485752/1‑151 (MQ=255)     accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGc                                                                                                                                            <  2:568674/151‑1 (MQ=255)         accgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTgg                                                                                                                                        >  1:910593/1‑151 (MQ=255)              cAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACg                                                                                                                                   >  2:1031648/1‑151 (MQ=255)                 aaCGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGaaa                                                                                                                                <  1:28974/151‑1 (MQ=255)                   cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                              >  2:50574/1‑151 (MQ=255)                     gTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTaa                                                                                                                            <  2:275512/151‑1 (MQ=255)                       aaGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAAc                                                                                                                           <  1:42051/150‑1 (MQ=255)                          cAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                       >  2:312661/1‑151 (MQ=255)                           aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                       >  1:288028/1‑150 (MQ=255)                           aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGcc                                                                                                                       <  2:157765/150‑1 (MQ=255)                                     tCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTcagc                                                                                                              >  1:1088690/1‑149 (MQ=255)                                         cagATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCAGGCCCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTa                                                                                                          <  2:321493/149‑1 (MQ=255)                                           gttGGCGGCGCACCGCCACCTCTTTAACACAGG‑TCATTATTTTTGCGCCAACCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACCCTG‑TTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTac                                                                                                         <  2:714745/146‑1 (MQ=255)                                           gATGGCGGCGCAACGGCCCCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACCCTG‑TTCATTCCCCCCCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTaca                                                                                                        <  2:582078/149‑1 (MQ=255)                                              ggcggcGCAACGGCACCTCTTTAACACAGG‑TCATTTTTGTGGCGCCAAACCATGCGCAAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                    <  2:857/150‑1 (MQ=255)                                              ggcggcGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaa                                                                                                      <  2:157598/148‑1 (MQ=255)                                              ggcggcCAAACGGCCCCTCTTTAACACAGG‑TCATTATTTTTGCGCCAAACCATGCGCCAGGGCTGTAAAAAGGAAAATTGCGGATCACCCTG‑TTCATCCCCCCCCACGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTa                                                                                                          <  2:154722/144‑1 (MQ=255)                                               gcggcgCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                    >  1:667775/1‑149 (MQ=255)                                                   cgcAACGGCCCCTTTTTAACACAGT‑TCATTATTTTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCCCACTG‑TTCATTCCCCCCCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                               <  2:966235/150‑1 (MQ=255)                                                   cgcAACGGCCCCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCCCAAAGGCTGTAAAAAGGAAAATTGCGGATCACACTT‑TTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCg                                                                                               <  2:15406/150‑1 (MQ=255)                                                         ggCACCTCTTTAACACAGG‑TCATTATTTTTGCGCAAAACCATGCGCCAAGGCTGTAAAAAGGCAAATTGCGGATCACCCTG‑TTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTg                                                                                        <  2:676037/151‑1 (MQ=255)                                                             ccTCTTTAACACAGG‑TCATTATTGTTGCGCAAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAg                                                                                    <  2:128783/151‑1 (MQ=255)                                                              ctctTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                                   >  1:550551/1‑151 (MQ=255)                                                                 tttAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                >  1:229129/1‑151 (MQ=255)                                                                 tttAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                >  2:10081/1‑151 (MQ=255)                                                                  ttAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                >  2:1002258/1‑150 (MQ=255)                                                                   tAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                                >  2:1080303/1‑149 (MQ=255)                                                                    aaacgcggttgattgttGTGCCGCAAAACCACCCGCCAAGGATGAAATAAGGAAAATGGCGGACCACACTC‑TTCGTCCCCCCCCACGGCGGGGGAGCAAAATAAACCGCCACGCTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAg                                                                              <  1:128029/140‑1 (MQ=255)                                                                    aaCACAGG‑TCTTTTTTGTTGCGCCAAACCAGGCGCCAGGGCTGTAAAAAGGAAATTTGCGGATCCCACTT‑TTCTTCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAg                                                                              <  2:3392/150‑1 (MQ=255)                                                                      ctcaGG‑TAATTTTTGTTGGCCCAACCCATGCGCCAAGGCTGAAAAAGGGAAAATTGCGGATAACACTT‑TTCATTCCCCCCCACGGCGTGGAAAAAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGAc                                                                            <  2:83102/148‑1 (MQ=255)                                                                      cacaGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                           >  2:1041922/1‑151 (MQ=255)                                                                          gg‑TCATTATTGTTGCCCCAAACCATGCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTG‑TTCATCCCCCACCACGGGGTGGGAAGGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                        <  1:798773/150‑1 (MQ=255)                                                                             tCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                      >  1:305089/1‑150 (MQ=255)                                                                             tCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                      >  2:632801/1‑150 (MQ=255)                                                                             tCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaa                                                                      >  1:773689/1‑150 (MQ=255)                                                                                ttattGTTCCGCAAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCCCCCGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaa                                                                     <  2:715121/148‑1 (MQ=255)                                                                                 ttttgttgCGCAAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCTTTCCCCCCCACGGCGGGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                                  <  2:1071084/148‑1 (MQ=255)                                                                                   tttttgCGCCAAACCATGCGCCAAGGCTGTAAAAAGGCAAATTGCGGATCCCCCTG‑TTCATTCCCCCCCACGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                <  2:26584/147‑1 (MQ=255)                                                                                   ttgttgCGCCAGACCATGCCCAAAGGCTATAATAAGGAAAATTGCGAATCACACTC‑TCCATCCCCCCCCACGGCGTGGAAACGAATTTAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAttt                                                                <  2:750266/150‑1 (MQ=255)                                                                                    tggagcgcCAGCCCAGGCGCAGAGGCTATAATAAGGAAAAGTGCGGACCACACTG‑TTCCTCCCCCCCGGCGGCGGGGGAACGAAATTAGCCGCCGCACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                               <  1:583724/146‑1 (MQ=255)                                                                                     gttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtttt                                                               >  1:797136/1‑149 (MQ=255)                                                                                      ttgCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCTCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAtt                                                                 >  1:716656/1‑146 (MQ=255)                                                                                      ttgCGCCAAACCAGGCGCCAAGGCTGTAATAAGGAAAATTGCGAATCACACTT‑TTCATTCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAc                                                             <  2:511873/150‑1 (MQ=255)                                                                                      tggcgcCAACCCAGGCCCCAGGGCTGAAAAAAGGCAAATGGCGGACCACACTG‑TTCATCCCCCCCCACGGCGGGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAAt                                                                  <  1:179623/143‑1 (MQ=255)                                                                                              acccATGCGCCAGGGTTGTAAAAAGGAAAATTGCGGATCACCTTG‑TTCATCCCCCCCCACGGGGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGa                                                    <  2:911980/149‑1 (MQ=255)                                                                                              aaCCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                     >  2:725805/1‑150 (MQ=255)                                                                                                 caggcgcCAAGGCTTAAAAAAGGAAAGTGGCGGACCACCCTG‑TTCATCCCCCCCCGCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGaa                                                   <  2:220062/146‑1 (MQ=255)                                                                                                  aTGCCCAAAGGCTGAAATAAGACAAAGTGCGAATCACTCTG‑TTCTCCCCCCACCGCGGCGGGGGAAAGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                                 <  1:885366/150‑1 (MQ=255)                                                                                                    gcgcCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAAGAAAt                                                 >  1:321546/1‑148 (MQ=255)                                                                                                        cAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAt                                          >  1:327411/1‑151 (MQ=255)                                                                                                           ggCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                       <  2:894748/151‑1 (MQ=255)                                                                                                           ggCTGAAATAAGAAAAATGGCGTATCCCACTG‑TTCATTCCCCACCCGGGCGTGGGACAGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTa                                        <  1:435276/150‑1 (MQ=255)                                                                                                            gttGTAAAAAGGAAAATTGCGGATCACCCTG‑TTCATCCCCCCCCACGGGGGGGAAACAAAATTACCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAAc                                      <  1:151100/149‑1 (MQ=255)                                                                                                            gCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAAc                                      >  1:243420/1‑151 (MQ=255)                                                                                                            gCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAAc                                      >  1:1068735/1‑151 (MQ=255)                                                                                                             cTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACg                                     >  2:1004755/1‑151 (MQ=255)                                                                                                                 aaaaaGGCAAATTGCGAATCACACTT‑TTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                                  <  2:880991/147‑1 (MQ=255)                                                                                                                  aaacggaaaaTTGCGGATCAC‑CTGTTTCATCCCCCCCCCCGGCGTGGAAACAAAATCAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                <  1:119102/144‑1 (MQ=255)                                                                                                                   taaGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                                >  2:538170/1‑150 (MQ=255)                                                                                                                         aaaTGGCGGGCCACCCTG‑TTCATCCCCCACCCGGGGGTGGGAGGGAAATAAGCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                         <  1:781878/151‑1 (MQ=255)                                                                                                                            ttGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac                       >  2:8672/1‑150 (MQ=255)                                                                                                                            ttGCGGATCACACTG‑TTCATCCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac                       <  2:415911/150‑1 (MQ=255)                                                                                                                              gggaaTCACACGG‑TTCATCCCCCCCCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                     <  2:988775/146‑1 (MQ=255)                                                                                                                                   tCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGt               >  2:853284/1‑151 (MQ=255)                                                                                                                                   tCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGt               >  2:301195/1‑151 (MQ=255)                                                                                                                                       acTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg           >  1:226200/1‑151 (MQ=255)                                                                                                                                        cTG‑TTCTTCCCCCCCCACGGGGTGGAAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTgg                <  1:222525/145‑1 (MQ=255)                                                                                                                                         tG‑TTCATCCCCCGCCACGGCGTGGGAACAAAATTATCCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTg           <  2:257349/149‑1 (MQ=255)                                                                                                                                            ttCATCCCCCCCCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCt            <  1:567610/146‑1 (MQ=255)                                                                                                                                            ttCATCCCCCCCCACGCAGTGGGAAAAAATTTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTgg          <  2:37926/148‑1 (MQ=255)                                                                                                                                               aTTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt      >  2:103727/1‑149 (MQ=255)                                                                                                                                                ttCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat   >  1:539330/1‑151 (MQ=255)                                                                                                                                                 cccccACCACGGTGTGGAACCGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGat   <  2:288075/149‑1 (MQ=255)                                                                                                                                                  ccccccccCGGGGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGatc  <  2:297828/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  AGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGG‑TCATTATTGTTGCGCCAAACCATGCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTG‑TTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGTGATC  >  NC_000913/1174463‑1174755

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