Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	715,755	A→C	5.6%	Q115P (CAG→CCG)	ybfP →	lipoprotein YbfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,755	0	A	C	5.6%	250.1 / 3.2	83	Q115P (CAG→CCG)	ybfP	lipoprotein YbfP
Reads supporting (aligned to +/- strand):  ref base A (37/40);  new base C (2/3);  total (39/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACAC  >  NC_000913/715606‑715901                                                                                                                                                      |                                                                                                                                                   tgATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCa                                                                                                                                                    >  1:1043067/1‑150 (MQ=255)   aTAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAgg                                                                                                                                                  <  1:987869/150‑1 (MQ=255)    tAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAggg                                                                                                                                                 <  2:64704/150‑1 (MQ=255)          aCCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAACAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGGTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGcc                                                                                                                                           <  2:1892281/150‑1 (MQ=255)            cTGCTAATAACCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAATGTCGGTTTAAATCCCGGGTGTCAATGCTACACCTCATTTTTTTTGCGGGCTTGCCCTGGCaaagccg                                                                                                                                         <  2:1418063/150‑6 (MQ=255)              gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGc                                                                                                                                       >  1:31172/1‑150 (MQ=255)              gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGc                                                                                                                                       >  1:820308/1‑150 (MQ=255)              gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGc                                                                                                                                       <  1:1951997/150‑1 (MQ=255)              gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgca                                                                                                                                       <  2:1064164/150‑3 (MQ=255)              gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgca                                                                                                                                       <  2:143695/150‑3 (MQ=255)                  ataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTc                                                                                                                                   >  2:1697795/1‑150 (MQ=255)                   taCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                  >  1:699906/1‑150 (MQ=255)                   aaCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACCTACTTATCTTCATAATATGTCAGCACAAGAAAAGGGCGGGTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                  <  1:317839/149‑1 (MQ=255)                     cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGAc                                                                                                                                >  2:420778/1‑150 (MQ=255)                                cAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCATTCGACCGATAGTCAAt                                                                                                                     <  2:1762823/150‑1 (MQ=255)                                cAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAAt                                                                                                                     >  2:1037385/1‑150 (MQ=255)                                 aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATg                                                                                                                    <  1:1518620/150‑1 (MQ=255)                                 aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATg                                                                                                                    >  2:1949789/1‑150 (MQ=255)                                    tGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAgg                                                                                                                 >  2:208315/1‑150 (MQ=255)                                              gAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCAAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAg                                                                                                       >  2:58368/1‑150 (MQ=255)                                                aCGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCg                                                                                                     <  2:529870/150‑1 (MQ=255)                                                 cGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGa                                                                                                    >  1:807683/1‑150 (MQ=255)                                                 cGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGa                                                                                                    >  1:814432/1‑150 (MQ=255)                                                    aaaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                                 <  2:1872478/150‑1 (MQ=255)                                                    aaaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                                 >  1:483896/1‑150 (MQ=255)                                                     aaaGAACTTACTTATCTTCATAATATGTCCGCACAAGAAAAGGTCGGTTTAAATCCCGGTTCGCAATGCTACACATCATTTTTTATGCGGGTTTTCCCGGGCAAGCCCGGCACACGACCGATAGTCAATGAGGATTATGTAAGCCAAAGc                                                                                                >  1:1243652/1‑150 (MQ=255)                                                     aaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGc                                                                                                >  2:1642367/1‑150 (MQ=255)                                                     aaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTATTTATGCGGGTTTGCCAGCGCAACCCCGCCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGc                                                                                                >  1:1243657/1‑150 (MQ=255)                                                       aGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCgg                                                                                              >  2:833740/1‑150 (MQ=255)                                                        gAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTTGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGATGATTATGTCAGCGAAAGCGGa                                                                                             <  1:1348580/150‑1 (MQ=255)                                                              acttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGTGTTTGCAAGTGCAAGACCGTCACTCGACCTATAGTAAATGATTATTATGTCAGAGAAATAGGAtttttt                                                                                       <  2:1468446/150‑1 (MQ=255)                                                              acttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATg                                                                                                                    >  2:1966423/1‑121 (MQ=255)                                                              acttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATg                                                                                                                    <  1:1966423/121‑1 (MQ=255)                                                              acttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGAtttttt                                                                                       <  2:1821713/150‑1 (MQ=255)                                                               cttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTg                                                                                      <  1:1037385/150‑1 (MQ=255)                                                               attaTCTTCATAATATGTCAGCACAAGAAAAGGGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTg                                                                                      <  1:1929774/149‑1 (MQ=255)                                                                 taTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGc                                                                                                                                                >  1:964112/1‑90 (MQ=255)                                                                 taTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGc                                                                                                                                                <  2:964112/90‑1 (MQ=255)                                                                   tCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTc                                                                                  <  1:957726/150‑1 (MQ=255)                                                                     ttCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                                >  2:296141/1‑150 (MQ=255)                                                                     ttCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                                >  1:161427/1‑150 (MQ=255)                                                                         taataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatgat                                                                            >  1:1140563/1‑150 (MQ=255)                                                                         taataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatgat                                                                            <  1:79079/150‑1 (MQ=255)                                                                          aataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGatgatg                                                                           >  2:446423/1‑150 (MQ=255)                                                                          aataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGatgatg                                                                           <  1:52985/150‑1 (MQ=255)                                                                           atatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGc                                                                          <  1:376743/150‑1 (MQ=255)                                                                           atatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGc                                                                          <  2:1440834/150‑1 (MQ=255)                                                                              tGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATg                                                                       >  1:845917/1‑150 (MQ=255)                                                                                 cAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACATCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCCGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCg                                                                    >  2:362490/1‑150 (MQ=255)                                                                                   gCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCggg                                                                  >  2:1738898/1‑150 (MQ=255)                                                                                    cacaAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGa                                                                 >  1:596238/1‑150 (MQ=255)                                                                                      caAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGatt                                                               >  1:1360641/1‑150 (MQ=255)                                                                                          aaaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                           <  2:807683/150‑1 (MQ=255)                                                                                          aaaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                           <  2:483896/150‑1 (MQ=255)                                                                                          aaaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                           <  2:814432/150‑1 (MQ=255)                                                                                          aaaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                           <  2:1043067/150‑1 (MQ=255)                                                                                            aaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGTGATTATTGAg                                                         <  2:1024341/150‑1 (MQ=255)                                                                                             aGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGt                                                        <  1:459160/150‑1 (MQ=255)                                                                                                   gTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCt                                                                                   <  2:857317/117‑1 (MQ=255)                                                                                                   gTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCt                                                                                   >  1:857317/1‑117 (MQ=255)                                                                                                      tAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCt                                               >  1:1422417/1‑150 (MQ=255)                                                                                                       aaaTCCCTGTTGGCAATGCTACAACTCATTTTTTATGCGGGTTTGCCTTTGGAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCTGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGATg                                              <  2:1243657/150‑1 (MQ=255)                                                                                                       aaaCCCCGTTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTGGACCGTTAGTCAATGAGGTTTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGGCGGGATTATTGAGTTACGACGCTg                                              <  2:1243652/150‑1 (MQ=255)                                                                                                        aaTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGc                                             <  2:672747/150‑1 (MQ=255)                                                                                                          tCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCa                                           >  1:381829/1‑150 (MQ=255)                                                                                                              ggTTGGCAATGCTACACCTCATTTTTTATTCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCa                                       <  1:590572/150‑1 (MQ=255)                                                                                                                       tGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGc                              <  2:1140563/150‑1 (MQ=255)                                                                                                                         cTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                            <  1:296141/150‑1 (MQ=255)                                                                                                                          tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAg                           >  1:1065796/1‑150 (MQ=255)                                                                                                                          tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAg                           <  1:119801/150‑1 (MQ=255)                                                                                                                              ccTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAg                       <  2:845917/150‑1 (MQ=255)                                                                                                                               cTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa                      >  1:818472/1‑150 (MQ=255)                                                                                                                                  aTTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAtt                   >  2:876073/1‑150 (MQ=255)                                                                                                                                     ttttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                >  1:627148/1‑150 (MQ=255)                                                                                                                                     ttttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                <  1:397332/150‑1 (MQ=255)                                                                                                                                       ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGt                                                        >  1:1496612/1‑108 (MQ=255)                                                                                                                                       ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGt                                                        <  2:1496612/108‑1 (MQ=255)                                                                                                                                        tATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc             <  1:448562/150‑1 (MQ=255)                                                                                                                                         aTGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg            >  1:1527087/1‑150 (MQ=255)                                                                                                                                          tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa           <  1:362490/150‑1 (MQ=255)                                                                                                                                          gccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa           >  2:794118/3‑150 (MQ=255)                                                                                                                                              ggTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATaaa       >  1:592036/1‑150 (MQ=255)                                                                                                                                               gTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAAt      >  2:152156/1‑150 (MQ=255)                                                                                                                                                   gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac  <  2:31172/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   TGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACAC  >  NC_000913/715606‑715901

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