Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,174,605	0	T	C	39.2%	127.0 / 11.4	79	G187G (GGA→GGG)	ycfT	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (35/13);  new base C (0/31);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 6.10e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGT  >  NC_000913/1174455‑1174751                                                                                                                                                        |                                                                                                                                                   aaaCAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                                                                                                                    <  1:1024152/151‑1 (MQ=255) aaaCAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                                                                                                                    <  2:347588/151‑1 (MQ=255) aaaCAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCAtt                                                                                                                                                    >  1:166418/1‑151 (MQ=255)        cAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCccacc                                                                                                                                             <  1:1124825/151‑1 (MQ=255)             accgaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGc                                                                                                                                        <  2:1039308/151‑1 (MQ=255)                gaccgCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTg                                                                                                                                     >  1:1082986/1‑151 (MQ=255)                    gCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGa                                                                                                                                  >  2:1335731/1‑150 (MQ=255)                           cGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAAtt                                                                                                                          <  1:37379/151‑1 (MQ=255)                                aGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCg                                                                                                                     >  1:636947/1‑151 (MQ=255)                                   aaaGAAGCCGTCAGCAGATGGCGGCGCAACGGCCCCTCTTTAACACAGGGCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCca                                                                                                                  <  1:107454/151‑1 (MQ=255)                                   aaaGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCca                                                                                                                  <  2:387728/151‑1 (MQ=255)                                    aagaagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCcac                                                                                                                 >  2:429293/1‑151 (MQ=255)                                       cagCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTACATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACt                                                                                                              >  1:273623/2‑151 (MQ=255)                                                  agATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTaca                                                                                                    >  1:670101/1‑150 (MQ=255)                                                  agATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTaca                                                                                                    >  1:921532/1‑150 (MQ=255)                                                   gATGGCGCCGCAACGCCCCCTCTTTAACACAGGTCATTATTGTGGCCCAAAACCAT‑GCGCCAAGGCTGTAATAAGGAAAATTGGGGATCACCCTGTTCTTCCCCCCCCCCGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTacaa                                                                                                   <  2:898352/150‑1 (MQ=255)                                                      ggcggcGCAACGGCCCCTCTTTAACACGGGTAATTATTGTGGCGCCAAACCAT‑GCGCAAAGGCTGAAATAAGCAAATTTGGGGTCCACTCTGTTCATCCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGt                                                                                                       <  2:753635/143‑1 (MQ=255)                                                      ggcggcGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                >  1:1133688/1‑150 (MQ=255)                                                      ggcggcGAAACGGCTCCTCTTTAACACAGGTAATTTTTGTTGCGCCAAACCTG‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACCCTTTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTacaaaca                                                                                                <  2:994985/150‑1 (MQ=255)                                                        cggcgCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAgg                                                                                             >  2:1223552/1‑151 (MQ=255)                                                         ggcaCAACGCCCCCTCTATACCACAGTTTATTATTTTGGGGCAAAACCAG‑GCCCCAAGGCTGTAATAAGAAAAATTGCGGACCCCCCTTTCCATCCCCCCCCACGGCGGGGGAAAGAAATTAGCCGCCACACTCAGCAGTACAAACAAGGc                                                                                            <  1:1245757/151‑1 (MQ=255)                                                               aCGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAAt                                                                                        >  2:977684/1‑149 (MQ=255)                                                                cGGCCCCTTTTTAACACAGGTCTTTTTTGTTGCGCAAACCCAT‑GCCCAAAGGCTGTAATAAGGAAAATTGCGGATCCCACTGTTCATCCCCCACCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATa                                                                                       <  2:1083767/149‑1 (MQ=255)                                                                cGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGt                                                                                     >  2:1319012/1‑151 (MQ=255)                                                                       tctTTAACACAGGTCCTTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgc                                                                               >  2:1348216/1‑150 (MQ=255)                                                                          ttAACACAGTTCATTATTGTGGCGCCAACCCAT‑CCCCCAAGGCTGAAATAAGGCAAATTGCGGATCACACTGTTCATCCCCCCCCCCGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAgcg                                                                              <  1:61600/148‑1 (MQ=255)                                                                          ttAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                            >  2:967667/1‑150 (MQ=255)                                                                           tAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                            >  2:164223/1‑149 (MQ=255)                                                                            aaCACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGcc                                                                            >  1:387740/1‑148 (MQ=255)                                                                            aaCACAGGTAATTATTGTGGCGCAAAACCAT‑GCGCAAAGGTTGTAATAAGGAAAATTGGGGATCCCACTGTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCa                                                                           <  2:424659/149‑1 (MQ=255)                                                                             acacAGTTCATTATTGTTGCGCCAAACCAG‑GCCCCAAGGCTGTAATAAGGAAAATGGCGGATCACACTTTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGa                                                                         <  2:1228218/150‑1 (MQ=255)                                                                              cacaGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                       >  1:564846/1‑151 (MQ=255)                                                                               acaGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                      >  1:122491/1‑151 (MQ=255)                                                                               acaGGTAATTTTTGTGGCGCAAACCCAT‑GCGCCAAGGCTGAAAAAAGGAAAATGGCGGATCCCACTGTTCATCCCCCCCCACGGGGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                       <  1:1065905/150‑1 (MQ=255)                                                                                caGGTCTTTTTTTTTGCGCC‑AACCATGGCGCCAGGGTTGAAAAAAGGAAAAGGGCGGATCACATTGTTCATCCCCCCCCAGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                       <  1:863547/149‑1 (MQ=255)                                                                                  ggTTATTTTTGTTGCGCCAAACCAC‑GCGCCAAGGATAAAAGAAGGAAAATGGCGAATCACACTGTTAATCCCCCCCCACGGCGTGGAAAAGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACgg                                                                      <  1:472650/148‑1 (MQ=255)                                                                                  ggTCATTTTTGTTGCGCCAAACCAT‑GCGCAAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCCCCCCGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGc                                                                     <  2:309694/149‑1 (MQ=255)                                                                                  ggTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCt                                                                    >  2:716005/1‑150 (MQ=255)                                                                                  ggTCATTATTGTGGCGCCAACCCAT‑GCGCCAAGGCTGTAAAAAGGAAAATTGCGGATCACCCTGTTCATCCCCCACCCCGGCGTGGAAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACg                                                                       <  2:1198793/147‑1 (MQ=255)                                                                                       ttattGTTGCGCAAAACCAG‑GCCCCAGGGCTGTAAAAAGGCAAATTGCGGATCACACTGTTCATCCCCCCCCACGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaaa                                                               <  2:300123/150‑1 (MQ=255)                                                                                       ttattGTGGCGCAAAACCAG‑GCGCCAAGGCTTAAAAAAGGAAAATGGCGAATCACACTGTTCATCCCCCCCCCCGGCGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                <  1:831592/149‑1 (MQ=255)                                                                                        tattTTTGCGCCAAACCAT‑GCGCCAAGGCTGTAAAAAGACAATTTGCGACCCACACTTTTCATTCCCCCCCACGGGGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTaaaa                                                                <  1:153414/148‑1 (MQ=255)                                                                                            gttgCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGAAAATTGCGGATCACACTGTTCATCCCCCCCCACGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTa                                                          <  2:944625/150‑1 (MQ=255)                                                                                             ttgCGCCAACCCTT‑GCGCCAAGGCTGTAAAAAGGCAATTGGCGGACCCCCCTTTTCATTCCCCCCCAGGGCGGGGGAACGAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTa                                                          <  1:69684/149‑1 (MQ=255)                                                                                                 gcCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                    >  1:170878/1‑151 (MQ=255)                                                                                                 gcCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATAcc                                                    >  2:970745/1‑151 (MQ=255)                                                                                                   cAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCAc                                                  >  1:433306/1‑151 (MQ=255)                                                                                                    aaaCCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                 >  2:73426/1‑151 (MQ=255)                                                                                                    aaaCCAG‑GCCCAAGGGCTGTAAAAAGAAAAATTGCGGATCCCCCTGTTCATCCCCCACCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCa                                                   <  1:624108/149‑1 (MQ=255)                                                                                                     aaCCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACg                                                 >  1:690809/1‑150 (MQ=255)                                                                                                        cAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGCCGGCTAAAAATTTTACATACCACGaaa                                              >  1:1179591/1‑150 (MQ=255)                                                                                                          t‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAAt                                             >  2:240291/1‑149 (MQ=255)                                                                                                             cgcCAAGGCTGTAATAAGGAAAATTGCGGATCACCCTGTTCATCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                                          <  2:232997/150‑1 (MQ=255)                                                                                                             ccccAAGGCTGTAATAAGGAAAATTGCGGATCACCCTGTTCATCCCCCCCCCCGGCGTGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGa                                          <  2:1378395/148‑1 (MQ=255)                                                                                                                cAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACaa                                       >  1:1244964/1‑150 (MQ=255)                                                                                                                  agcctttaatagGGCAAATTGCGGACCACCCTTTTCATCCCCCCCCCCGGGGTGGGAACGAAATACACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAAtt                                     <  1:307978/144‑1 (MQ=255)                                                                                                                   ggCTGTAATAAGGAAAATTGGGAATCACACTGTTCATCCCCCCCCACGGCGGGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                   <  2:595741/151‑1 (MQ=255)                                                                                                                   ggCTGAAAAAAGGCAAATTGCGGATCACACTGTTCATCCCCCCCCACGGGGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTaa                                   <  1:1358510/151‑1 (MQ=255)                                                                                                                         aataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                              >  2:333065/1‑150 (MQ=255)                                                                                                                         aataaGGAAAATTGCGGATCCCCCTGTTCATCCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCa                               <  2:950280/149‑1 (MQ=255)                                                                                                                         aataaGGAAAATTGCGGATCACACTGTTCATTCCCCCCCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCAt                              <  2:289964/150‑1 (MQ=255)                                                                                                                         aaaaaGGCAAATTGCGGATCACACTGTTCATCCCCCCCCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATa                             <  2:46222/148‑1 (MQ=255)                                                                                                                          ataaGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                            >  2:767432/1‑151 (MQ=255)                                                                                                                            aaGGCAAATGGCGGATCACACTGTTCATCCCCCCCCCGGGCGTGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACa                           <  1:382696/150‑1 (MQ=255)                                                                                                                              ggCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATAc                            >  2:1230716/1‑147 (MQ=255)                                                                                                                                cAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAc                      >  1:457605/1‑151 (MQ=255)                                                                                                                                 aatttGCGGATCACCCTTTTCTTTCCCCACCCGGGGGTGAGAGCAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGAt                        <  2:452779/145‑1 (MQ=255)                                                                                                                                 aaaTTGCGGATCACACTGTTCATTCCCCCCCCCGGCGGGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                       <  2:215017/149‑1 (MQ=255)                                                                                                                                 aaaTTGCGGATCACACTGTTCATCCCCCACCCCGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATa                       <  2:401787/149‑1 (MQ=255)                                                                                                                                 aaaTTGCGAACCACCCTGTTCTTCCCCCACCACGGGGTGGGAAGAAAATGAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAg                          <  1:1239783/146‑1 (MQ=255)                                                                                                                                  acttGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATAcc                     >  1:1258790/3‑150 (MQ=255)                                                                                                                                  aaTTGCGGATCACACTTTTCATCCCCCACCACGGGGGGGGAACAAAATAAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACca                    <  2:1145228/151‑1 (MQ=255)                                                                                                                                    ttGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACcac                   >  2:44452/1‑150 (MQ=255)                                                                                                                                      gCGGACCACACTGTTCATTCCCCCCCACGGCGTGGAAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAg                <  2:1067461/151‑1 (MQ=255)                                                                                                                                        ggATCCCCCTGTTCTTCCCCCACCCGGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACaa                 <  2:545825/148‑1 (MQ=255)                                                                                                                                        ggAGCACACTTTTCTTTCCCCCCCAGGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGc               <  1:819804/150‑1 (MQ=255)                                                                                                                                         gATCACACTGTTCATCCCCCCCCACGGCGGGGAAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTg             <  2:637133/151‑1 (MQ=255)                                                                                                                                              cacTGTTCTTCCCCCCCCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGc         <  1:1329213/150‑1 (MQ=255)                                                                                                                                                 tGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTAAATACCACGACATAGACAATTAACGCAAACAGATACCAAAAGCTGGTGCTgg      >  1:876121/1‑150 (MQ=255)                                                                                                                                                   ttCATCCCCCACCACGGCGTGGGAACAAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGc     <  2:932816/149‑1 (MQ=255)                                                                                                                                                    tCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCgg   >  2:863950/1‑150 (MQ=255)                                                                                                                                                      aTTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCgg   >  2:660098/1‑148 (MQ=255)                                                                                                                                                       ttCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCgg   >  2:718833/1‑147 (MQ=255)                                                                                                                                                       ttCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGt  >  1:601467/1‑148 (MQ=255)                                                                                                                                                        |                                                                                                                                                   AAACAACCAGGCACCGACCGCCAGAACGGTAAGCAAAGAAGCCATCAGCAGATGGCGGCGCAACGGCACCTCTTTAACACAGGTCATTATTGTTGCGCCAAACCAT‑GCGCCAAGGCTGTAATAAGGCAAATTGCGGATCACACTGTTCATTCCCCACCACGGCGTGGGAACGAAATTAACCGCCACACTCAGCAGTACAAACAAGGCGAATAGTGGCAGCGCCAGACGGCTAAAAATTTTACATACCACGAAATAGACAATTAACGCATACAGATACCACAAGCTGGTGCTGGCGGT  >  NC_000913/1174455‑1174751

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