Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,435,727	0	G	A	41.9%	90.1 / 46.5	74	intergenic (‑108/+33)	uspF/ompN	stress‑induced protein, ATP‑binding protein/outer membrane pore protein N, non‑specific
Reads supporting (aligned to +/- strand):  ref base G (36/7);  new base A (0/31);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.95e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGAC  >  NC_000913/1435579‑1435874                                                                                                                                                         |                                                                                                                                                    gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAg                                                                                                                                                     <  2:648002/149‑1 (MQ=255) gAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGa                                                                                                                                                    >  1:1162895/1‑150 (MQ=255)   aTCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGaaaa                                                                                                                                                 <  2:645388/151‑1 (MQ=255)     cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGaaaaa                                                                                                                                                <  2:780103/150‑1 (MQ=255)     cGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGaaaa                                                                                                                                                 >  1:1098437/1‑149 (MQ=255)       gAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGcc                                                                                                                                             <  1:366700/151‑1 (MQ=255)             tCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCAt                                                                                                                                       >  1:629297/1‑151 (MQ=255)               gATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATa                                                                                                                                      >  2:153294/1‑150 (MQ=255)                   gggACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCTCGCCATAACAGc                                                                                                                                 <  2:676398/151‑1 (MQ=255)                    ggACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGc                                                                                                                                 >  2:299571/1‑150 (MQ=255)                       cAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGc                                                                                                                             >  2:240104/1‑151 (MQ=255)                        aaaaaTCGTTCTGTTCAAAAACCCTCT‑TTTTTGAGTATCCGAAAAAATTATAACGCGAGAGGACACTTTTTGGATAATGCCGGATTCCATTTTTTAAATTGAGAT‑CAAAGAAGG‑AAAAAT‑AACAAAAAAAAGCCCGCCATAACAGCGGGCa                                                                                                                            <  2:836551/148‑1 (MQ=255)                         aGAATGGTTCTGTTAAAAAAACCCCT‑TTTTTTATTACCCAAAAAAAGTATAACGCAAAATGACTTTTTTGGTTTAAAGAGGGAGTCCCCTTTTTTAATTTAAAT‑C‑AAGGAAG‑AAAAATGAAAAAAAAAAAGCCCGCCATAACAGCGGGCa                                                                                                                            <  1:121402/150‑1 (MQ=255)                          gAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCagg                                                                                                                          >  2:728823/1‑151 (MQ=255)                           attCGTTCTGTTCAAAAACCCTCT‑TGTTTTAGTACCCGAAAAAAGTGTACCGCAAGATGACATTTTTTGTGTAATGAGGGAGTTCACTTTTTTAATTTAGAT‑AAAAGGAGA‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCagg                                                                                                                          <  1:767521/148‑1 (MQ=255)                           aatggttgtttcaaaaaacCCTCT‑TGTTTTATTACCCGAATAAAGTGAAACGCAAGATGACATTTTTTGTGAAAGGAGTAATTTCACTTTTTTATTTGAGAC‑AAAAGGAGA‑AAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCagga                                                                                                                         <  2:243350/136‑1 (MQ=255)                            aTCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCagg                                                                                                                          >  1:1327357/1‑149 (MQ=255)                               gTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGAt                                                                                                                     >  1:1219877/1‑151 (MQ=255)                                tttttttCATAAACCCTCC‑TTTTTTAGTACCCGCATAAAGTGTACCACCGGATAACCCTTTTTGTATAATGGCGG‑GTTCCAGTTTTTAAATT‑GATGAAAAGGGGGAAAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCaggagg                                                                                                                       <  1:1130230/143‑1 (MQ=255)                                   tGTTCAAAAACCCCCC‑TTTTTTATTATCCGAATAAAGTGTAACCCCAGATGACTTTTTTTGTGTAATGACGGATTTCACATTTTTAATTTAGAT‑AAAAGGAGA‑AAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTa                                                                                                                  <  1:1020005/150‑1 (MQ=255)                                      tCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAAc                                                                                                              >  1:220513/1‑151 (MQ=255)                                        aTAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTg                                                                                                            >  1:1324670/1‑151 (MQ=255)                                          aaaaCCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGAt                                                                                                          >  2:928350/1‑151 (MQ=255)                                           aaaCCTCCTTGTTTTAGTATCCAAAAAAAGTGAAACGCAAGATGACATTTTTTGTGTAATGGGGGAGTGCAAATTTTTAATTTAGAA‑AAAAGGAGG‑AAAAAA‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGa                                                                                                           <  2:224705/150‑1 (MQ=255)                                                       tttAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCt                                                                                               >  2:222711/1‑150 (MQ=255)                                                        ttAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTa                                                                                              >  2:851844/1‑150 (MQ=255)                                                          aGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc                                                                                                      >  2:631983/1‑140 (MQ=255)                                                          aGTACCCGAAAAAAGTGTACCGCAAGATAACTTTTTTTGTGTAATGAGGGAGTTCATTTTTTTAATTAAGAT‑AAAAGGAGA‑AAAAAT‑AAGAAAAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTaa                                                                                             <  1:155185/149‑1 (MQ=255)                                                            tATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                         >  1:432648/1‑151 (MQ=255)                                                             aTCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCg                                                                                         >  1:1279606/1‑150 (MQ=255)                                                              tCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                                                                       >  2:901939/1‑151 (MQ=255)                                                               ccGCAAAAAGTGTACCGCCAGATCACACTTTTTGTAAAAGGAGGGATTACACTTTTTAAATTAAAAA‑AAAGGGAGA‑AAGAAA‑AACAAAAAAACGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGAc                                                                                       <  2:268757/150‑1 (MQ=255)                                                               ccGAATAAAGTGAAACGCAAGATGACTTTTTTTGTGTAATGACGGAGTTCACTTTTTTAATTTAGAT‑AAAAGGAGA‑AAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAAc                                                                                                      <  1:631993/135‑1 (MQ=255)                                                                cgaaTAAAGTGTAACGCCAGATGACTTTTTTTGTTTAATGACGGATTTCACTTTTTTAATTTAGAA‑AAAAGGAGG‑AAGAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACaa                                                                                     <  1:206480/148‑1 (MQ=255)                                                                  cATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                    >  2:195086/1‑150 (MQ=255)                                                                  cATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                    >  2:998519/1‑150 (MQ=255)                                                                  aaaaaaGTGTAACGCAAGATGACTTTTTTTGTGAAAGGAGGGATTTCAATTTTTTAATTAAGAA‑AAAAGGAGA‑AAAAAA‑AAAAAAAAAAGCCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGa                                                                                        <  1:360941/143‑1 (MQ=255)                                                                   aaaaaGTGTACCGCAAAATGCCGTTTTTTGTGAAATGACGGATTTCATTTTTTAAATTAAGAT‑AAAAGGAGA‑AAAAAA‑AAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAat                                                                                    <  1:1218080/147‑1 (MQ=255)                                                                      aaGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatca                                                                                >  2:882941/1‑150 (MQ=255)                                                                          gtAACCCGAGATGACTTTTTTTGTGAAAGGAGGGTTTTCAATTTTTAAATTTAGAA‑AAAAGGAGA‑AAAAAT‑AAAAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATatcatc                                                                              <  2:591285/148‑1 (MQ=255)                                                                           tAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                           >  1:1139814/1‑150 (MQ=255)                                                                            aaCCCCAGATGACTTTTTTTGTGTAATGGCGGAGTTCACTTTTTTAATTAAGAT‑AAAGGGAGA‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTa                                                                          <  1:134428/150‑1 (MQ=255)                                                                             aCGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGt                                                                           >  2:854990/1‑148 (MQ=255)                                                                                  agagGACATTTTTGGTGAAAGGGGGAATTTCACTTTTTTAATTAAAAG‑AAAAGGAGA‑AAGAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCg                                                                             <  1:1277596/141‑1 (MQ=255)                                                                                    aTGACCTTTTTTGGGAAAGGAGGGATTCCACTTTTTTAATTAAGAT‑AAAAGGAGA‑AAAAAT‑AAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGcc                                                                  <  2:689712/150‑1 (MQ=255)                                                                                    aTGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGcc                                                                  >  2:226782/1‑150 (MQ=255)                                                                                      gacttttttgGTGTAAGGAGGGAGTTCACTTTTTTAATTAAAAT‑AAAAGGAGA‑AAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGc                                                                   <  1:627197/142‑1 (MQ=255)                                                                                          cTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTg                                                              >  1:1268087/1‑148 (MQ=255)                                                                                           ttttgGTGAAAGGACGAAGTACATATTTTTAATTTAAAG‑AAAGGGAGA‑AAAAAT‑AGGAAGAAAAGGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCa                                                           <  2:141362/150‑1 (MQ=255)                                                                                            ttttGTGAAATGACGGAGTTCACATTTTTAATTAAAAT‑AAAAGGAGA‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAg                                                          <  1:742801/150‑1 (MQ=255)                                                                                               tgtgTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGt                                                       >  1:1298656/1‑150 (MQ=255)                                                                                                    aaTGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGc                                                 >  2:16645/1‑151 (MQ=255)                                                                                                      tGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTg                                               >  1:130606/1‑151 (MQ=255)                                                                                                          ggAGTTCACATTTTTAATTAAGAA‑AAAAGGAGG‑AAAAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTg                                               <  1:1277914/147‑1 (MQ=255)                                                                                                                    tttttAATTGAAAT‑AAAAGGAGG‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCa                                  <  1:1313234/150‑1 (MQ=255)                                                                                                                    tttttAATTAAGAA‑AAAGGGAAA‑AAAAAT‑AAAAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGt                                              <  2:514119/138‑1 (MQ=255)                                                                                                                       taaaTTTAAAA‑AAAGGGAGA‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACa                               <  2:1353759/148‑1 (MQ=255)                                                                                                                        ttatttaagacaaagggaga‑aAAAAA‑AAGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACa                               <  1:790239/135‑1 (MQ=255)                                                                                                                        tAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                             >  1:354553/1‑151 (MQ=255)                                                                                                                         aaTTTAGAA‑AAAAGGAGG‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                            <  1:487816/151‑1 (MQ=255)                                                                                                                         aaTTAAGAA‑AAAAGGAGA‑AAAAAT‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                            <  1:908951/151‑1 (MQ=255)                                                                                                                           tttAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTg                          >  1:482802/1‑151 (MQ=255)                                                                                                                           tttAGAA‑CAAAGGAGG‑AAAAAA‑AAGAAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGtt                           <  1:1025506/150‑1 (MQ=255)                                                                                                                            ttAGAT‑CAAAGGAGA‑AAAAAT‑AAGCAAAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGt                            <  2:248108/148‑1 (MQ=255)                                                                                                                              aGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtt                       >  2:77261/1‑151 (MQ=255)                                                                                                                                 t‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGAtttt                     >  2:1004284/1‑150 (MQ=255)                                                                                                                                  aaaagggaga‑aAAAAA‑AAGCAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTta                    <  2:137185/146‑1 (MQ=255)                                                                                                                                    ataggaag‑AAAAAT‑ACGAAAAAAAACCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAGCCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAgg                             <  2:379183/138‑1 (MQ=255)                                                                                                                                    aaaGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtat                   >  1:738846/1‑150 (MQ=255)                                                                                                                                    aaaGGAAA‑AAAAAA‑AAGAAAAAAAAACCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTtata                  <  2:350351/151‑1 (MQ=255)                                                                                                                                      aggagg‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt                >  1:447548/1‑151 (MQ=255)                                                                                                                                             aagaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGt                <  2:21805/145‑1 (MQ=255)                                                                                                                                              agaaT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACg           >  2:1319496/1‑149 (MQ=255)                                                                                                                                                      gCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGAc  >  2:972961/1‑151 (MQ=255)                                                                                                                                                      aaaGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGa   <  2:54956/148‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                                    GAATCGGAAATATCGATAGGGACAAGAATCGTTCTGTTCATAAAACCTCC‑TGTTTTAGTATCCGCATAAAGTGTAACGCCAGATGACACTTTTTGTGTAATGACGGAGTTCACATTTTTAATTTAGAT‑CAAAGGAGG‑AAGAAT‑AAGCAGAAAAAGCCCGCCATAACAGCGGGCAGGAGGATTTAGAACTGATAAACCAGACCTAAAGCGACAATATCATCGGTAGAGATGCCATTGGCAGCGTAGAAGCTGTCATCTTCATCCAACAGGTTGATTTTATAGTCAACGTAGGTGGAC  >  NC_000913/1435579‑1435874

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