Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	715,756	G→T	5.3%	Q115H (CAG→CAT)	ybfP →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,756	0	G	T	5.3%	234.7 / 2.2	84	Q115H (CAG→CAT)	ybfP	lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (35/44);  new base T (2/3);  total (37/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTA  >  NC_000913/715607‑715905                                                                                                                                                      |                                                                                                                                                      gATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAg                                                                                                                                                       <  2:490053/150‑1 (MQ=255)  aTAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAgg                                                                                                                                                      <  2:665947/150‑1 (MQ=255)   tAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAggg                                                                                                                                                     >  2:1084883/1‑150 (MQ=255)   ataCTTACCTGATATTACCCGCCTTTGGGTGTTGTCCTGCGGGAACGAAAAAGAACTTACTATTCTTCATAATATGTCCGCACAAGAAAAGGTCGGTTTAAATCCCGGTTTGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAggg                                                                                                                                                     <  1:140340/148‑1 (MQ=255)          ccTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGccc                                                                                                                                              <  1:810150/150‑1 (MQ=255)          ccTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGccc                                                                                                                                              >  2:29001/1‑150 (MQ=255)           ccgATAATACCCGACGCTGGCAGATGTCAGTCGGGAACGAAAAAGATCTTACTTATCTTCATAATATGTCAGCACAAGAAAAGCTCGGTTTAACTACCGGATCGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCATACCCg                                                                                                                                             <  1:152984/148‑1 (MQ=255)           cTGATAATACCCGCCGCTGGCAGTTGTCATGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCg                                                                                                                                             <  2:908360/150‑1 (MQ=255)            tGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCgg                                                                                                                                            >  2:145403/1‑150 (MQ=255)             gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGc                                                                                                                                           <  2:1308764/150‑1 (MQ=255)               taataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCAc                                                                                                                                         <  2:1312383/150‑1 (MQ=255)               taataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCAc                                                                                                                                         >  1:444721/1‑150 (MQ=255)                  taCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                      <  1:812906/150‑1 (MQ=255)                    cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGAc                                                                                                                                    >  2:716710/1‑150 (MQ=255)                          gcTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATTCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAATCCCGGCACTCGACCGATAg                                                                                                                              >  2:1028623/1‑150 (MQ=255)                            tGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCATCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTc                                                                                                                            <  2:1017721/150‑1 (MQ=255)                               cAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAAt                                                                                                                         >  2:567530/1‑150 (MQ=255)                                aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATg                                                                                                                        >  2:496397/1‑150 (MQ=255)                                 gTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCCCAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGa                                                                                                                       <  1:306121/150‑1 (MQ=255)                                  ttGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAg                                                                                                                      <  2:1232842/150‑1 (MQ=255)                                  gtGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAg                                                                                                                      <  1:756666/149‑1 (MQ=255)                                      cAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGAtt                                                                                                                  <  1:144357/150‑1 (MQ=255)                                       aggaggGAAAGAAAAAGAACTAACTTATCTTCATAATATGTCAGCAAAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGAAAGTCAATGAGGATTa                                                                                                                 <  1:1288664/150‑1 (MQ=255)                                       aGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGAAAGCCCGGCACTCGACCGATAGTCAATGAGGATTa                                                                                                                 <  2:1042576/150‑1 (MQ=255)                                        ggCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTAt                                                                                                                >  2:718678/1‑150 (MQ=255)                                         gCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATg                                                                                                               <  1:989753/150‑1 (MQ=255)                                            ggAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCa                                                                                                            <  1:984646/150‑1 (MQ=255)                                            ggAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCa                                                                                                            <  1:911710/150‑1 (MQ=255)                                               aCGAAACAGAACTTACATATCTTCATAATATGTCAGCCCAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCg                                                                                                         <  2:782968/150‑1 (MQ=255)                                               aCGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCACGGCACTCGACCGATAGTCAATGAGGATTATGTCAGcc                                                                                                         >  2:179715/1‑149 (MQ=255)                                                  aaaaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTATGTCAGCGaaa                                                                                                      >  1:84183/1‑150 (MQ=255)                                                   aaaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCTATGCTACACCTCAATTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                                     >  1:209104/1‑150 (MQ=255)                                                        aaCTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGAt                                                                                                <  1:1253780/150‑1 (MQ=255)                                                             acttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAACTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCTACCTATAGTCAATGAGGATTATGTCAGCGAAAGCGGAtttttt                                                                                           >  2:803662/1‑150 (MQ=255)                                                               ttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCTATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTgg                                                                                         <  1:121816/150‑1 (MQ=255)                                                               ttaTCTTCATAATATGTCACCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTgg                                                                                         <  2:41036/150‑1 (MQ=255)                                                                taTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTCTACATCCCGGTCGGCAATGCTACACCTCATTTTTTATGCGGGATTGCCAGGGCAAGCCCGGCACTCGACCGATACTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGc                                                                                        <  1:43632/150‑1 (MQ=255)                                                                  tCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAACTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTc                                                                                      >  2:1132115/1‑150 (MQ=255)                                                                   cTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCg                                                                                     >  2:237880/1‑150 (MQ=255)                                                                   cTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCg                                                                                     >  2:1155265/1‑150 (MQ=255)                                                                   attCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCg                                                                                     <  1:961482/149‑1 (MQ=255)                                                                     tCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgat                                                                                   >  1:916329/1‑150 (MQ=255)                                                                     acATAATATGTCACCACAAGAAAAGTGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgat                                                                                   <  1:1144349/149‑1 (MQ=255)                                                                      cATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatg                                                                                  >  1:23425/1‑150 (MQ=255)                                                                      cATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTGTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATTGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatg                                                                                  >  2:524557/1‑150 (MQ=255)                                                                       ataataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatga                                                                                 >  1:187216/1‑150 (MQ=255)                                                                        taataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatgat                                                                                >  2:423620/1‑150 (MQ=255)                                                                        taataTGTCAGCACAAGAAAAGGTCGGATTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatgat                                                                                <  2:63346/150‑1 (MQ=255)                                                                           tatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGTGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCa                                                                             <  2:254010/150‑1 (MQ=255)                                                                           tatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCa                                                                             <  1:403098/150‑1 (MQ=255)                                                                             tGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCTGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATg                                                                           <  2:553411/150‑1 (MQ=255)                                                                                 aGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCgg                                                                       <  2:392812/150‑1 (MQ=255)                                                                                     caAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTAGGCGGGTTTGCCAGGGCCAGCCCGGCACTCGTCCGATAGTCAATGAGGATGAGGGCTGCGAAAGAGGATTTTTTGGCTCGATGATGCATGTCGGGatt                                                                   >  2:162279/1‑150 (MQ=255)                                                                                         aaaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                               <  1:1202402/150‑1 (MQ=255)                                                                                         aaaaGGGCTGATTAACTCCCGGTTGGCAATGCTACACCTCATTTTTTATCCGGGTTTGCCATGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                               <  2:979934/150‑1 (MQ=255)                                                                                           aaGGTCGGTTTAAATGCCGGGTGGCAATGCTACAGCTGATTTTTTATGCGGGTTTGCCAGGGCTAGCCCGTCAGTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGGATGTCGGGATTATTGAg                                                             <  1:111093/150‑1 (MQ=255)                                                                                            aGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGt                                                            >  2:665988/1‑150 (MQ=255)                                                                                              gTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTa                                                          >  2:954586/1‑150 (MQ=255)                                                                                                 ggTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacga                                                       >  1:842973/1‑150 (MQ=255)                                                                                                  gTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacgac                                                      >  1:608649/1‑150 (MQ=255)                                                                                                     tAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGcg                                                   >  1:1203601/1‑149 (MQ=255)                                                                                                        aTCCCGTTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGAGTATGTCAGCGAAAGCGGATTTTTTGGCTCGAGGATGGATGTCGGGATTATTGAGTTACGACGCTGcc                                                <  1:736410/150‑1 (MQ=255)                                                                                                        aTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCTCTCGCCCGATAGTCAATGAGGATTATGTCAGCGAATGCGGATTTTTTGGCTCGTTGGTGGATGGCGGGTGTGTTGAGTTACGACGCGGcc                                                >  2:228173/1‑150 (MQ=255)                                                                                                        aTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGcc                                                <  1:496397/150‑1 (MQ=255)                                                                                                        aTCCCGGTTGGCAATGCTACACCTCATTTATTATGCGGGTTTGCCAGGGCCAGCCCGGCACTCGACCGATCGTCAATGAGGATTATGTCAGCGAAAGCGGATGTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGcc                                                <  1:562777/150‑1 (MQ=255)                                                                                                           ccGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGTTTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGt                                             >  1:972286/1‑150 (MQ=255)                                                                                                           ccGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGt                                             <  1:579209/150‑1 (MQ=255)                                                                                                            cGGTTGCCAATGCTACACCTCATTTATTATGCGGGTCTGCCAGGGCAAGCCAGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAACTCGGATTTTTTGGCTCGACGATGCCTGTCGGGATTATTGAGTTACGACGCTGCCAGTc                                            <  1:653955/150‑1 (MQ=255)                                                                                                             ggTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATTATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCa                                           <  1:75562/150‑1 (MQ=255)                                                                                                                  gCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaaa                                      >  2:283710/1‑150 (MQ=255)                                                                                                                       gCTACACCTCATTTTTTATGCGTGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAAAAAATTGcc                                 <  2:564104/150‑1 (MQ=255)                                                                                                                       gCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGcc                                 <  1:29001/150‑1 (MQ=255)                                                                                                                         tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTGGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCTTGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAg                               <  2:241054/150‑1 (MQ=255)                                                                                                                            acCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCaa                            <  1:1027999/150‑1 (MQ=255)                                                                                                                                    ttttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa                    <  1:30443/150‑1 (MQ=255)                                                                                                                                      ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAg                  <  2:497186/150‑1 (MQ=255)                                                                                                                                      agtgccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGAGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAg                  >  1:407578/6‑150 (MQ=255)                                                                                                                                       aaTGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                 <  1:642197/149‑1 (MQ=255)                                                                                                                                         tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGATAATTGCCAGCAAGAATTGAAAGCGa               >  1:278014/1‑150 (MQ=255)                                                                                                                                         tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa               <  1:716710/150‑1 (MQ=255)                                                                                                                                          gCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt              <  1:924127/150‑1 (MQ=255)                                                                                                                                          ccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATTCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGAGAt              >  2:1252713/2‑150 (MQ=255)                                                                                                                                                 tGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca       <  1:1084883/150‑1 (MQ=255)                                                                                                                                                  gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGCCGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac      >  2:1080523/1‑150 (MQ=255)                                                                                                                                                   ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca     >  2:413984/1‑150 (MQ=255)                                                                                                                                                   ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca     >  1:326341/1‑150 (MQ=255)                                                                                                                                                   ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca     >  1:1086020/1‑150 (MQ=255)                                                                                                                                                    cAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAt    >  1:962489/1‑150 (MQ=255)                                                                                                                                                     tgggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACAtt   <  1:257267/149‑1 (MQ=255)                                                                                                                                                      gggCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTa  >  1:839017/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTA  >  NC_000913/715607‑715905

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