Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	715,762	G→A	8.8%	K117K (AAG→AAA)	ybfP →	lipoprotein YbfP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	715,762	0	G	A	8.8%	266.6 / 12.3	91	K117K (AAG→AAA)	ybfP	lipoprotein YbfP
Reads supporting (aligned to +/- strand):  ref base G (50/33);  new base A (4/4);  total (54/37)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.01e-01

ACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAG  >  NC_000913/715615‑715911                                                                                                                                                    |                                                                                                                                                      aCCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGcc                                                                                                                                                     >  1:1610138/1‑150 (MQ=255)  ccTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGccc                                                                                                                                                    >  2:629904/1‑150 (MQ=255)    tGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCgg                                                                                                                                                  >  1:1154178/1‑150 (MQ=255)     gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGc                                                                                                                                                 <  2:1147508/150‑1 (MQ=255)     gATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgca                                                                                                                                                 <  1:236365/150‑3 (MQ=255)      ataataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCa                                                                                                                                                <  1:1394234/150‑1 (MQ=255)      ataataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgcat                                                                                                                                                <  1:672847/150‑4 (MQ=255)       taataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCAc                                                                                                                                               >  1:486298/1‑150 (MQ=255)       taataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCAc                                                                                                                                               >  1:172034/1‑150 (MQ=255)       taataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCAc                                                                                                                                               >  2:459261/1‑150 (MQ=255)        aataCCCGCCGCTGGCAGTTGTCACGCGGGAACGAAAAAGAACTTACTTATCTTCTTAATATGTCAGCACAAGAAAAGGACGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACt                                                                                                                                              >  2:697954/1‑150 (MQ=255)         ataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgcataaa                                                                                                                                             <  1:101748/150‑7 (MQ=255)          taCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                            <  2:1197706/150‑1 (MQ=255)            cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGAc                                                                                                                                          >  1:741385/1‑150 (MQ=255)            cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGAc                                                                                                                                          >  2:1497246/1‑150 (MQ=255)                ccgcTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                            <  2:529364/144‑1 (MQ=255)                ccgcTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                                            >  1:529364/1‑144 (MQ=255)                 cgcTGGCACTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGACGGTTTAAATCCCGGCTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATa                                                                                                                                     <  2:1550599/150‑1 (MQ=255)                    tGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTc                                                                                                                                  <  1:246566/150‑1 (MQ=255)                        aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATg                                                                                                                              >  1:1312770/1‑150 (MQ=255)                           tGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAgg                                                                                                                           <  2:925100/150‑1 (MQ=255)                           tGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAgg                                                                                                                           <  2:867146/150‑1 (MQ=255)                            gTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGa                                                                                                                          <  1:1079976/150‑1 (MQ=255)                              cAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGAtt                                                                                                                        >  1:411164/1‑150 (MQ=255)                              cAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGAtt                                                                                                                        <  2:897193/150‑1 (MQ=255)                              cAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGAtt                                                                                                                        >  1:166981/1‑150 (MQ=255)                                ggCGGGAACGAAAAAGACCTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTAt                                                                                                                      >  2:1085613/1‑150 (MQ=255)                                ggCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTAt                                                                                                                      >  2:131052/1‑150 (MQ=255)                                ggCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTAt                                                                                                                      >  2:34216/1‑150 (MQ=255)                                ggCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTAt                                                                                                                      >  2:79530/1‑150 (MQ=255)                                       aCGAAAAAGAACTTACTTATCTTCATAATATCTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGATAAACCACATTTTTTATGAGGGTTTGCCAGGGCAAGCCCGGCACTCGACCCATACTAAATGAGGATTATGTCAGCg                                                                                                               >  2:1419514/1‑150 (MQ=255)                                        cGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGa                                                                                                              >  2:583614/1‑150 (MQ=255)                                        cGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGa                                                                                                              >  2:930316/1‑150 (MQ=255)                                             aaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTGTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCg                                                                                                         <  2:819884/150‑1 (MQ=255)                                               gAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGa                                                                                                       >  2:1001935/1‑150 (MQ=255)                                                 acttacttATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGAtt                                                                                                     >  2:758663/1‑150 (MQ=255)                                                  cttacttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                                           <  2:960095/143‑1 (MQ=255)                                                  cttacttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                                           >  1:960095/1‑143 (MQ=255)                                                        taTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGc                                                                                              <  1:491016/150‑1 (MQ=255)                                                          tCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTc                                                                                            <  1:1397639/150‑1 (MQ=255)                                                           cTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCg                                                                                           <  2:1346676/150‑1 (MQ=255)                                                            ttCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                                          >  1:201104/1‑150 (MQ=255)                                                             tCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgat                                                                                         <  2:185951/150‑1 (MQ=255)                                                             tCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgat                                                                                         >  2:663650/1‑150 (MQ=255)                                                              cATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatg                                                                                        >  1:188493/1‑150 (MQ=255)                                                                  atatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGc                                                                                    <  1:663650/150‑1 (MQ=255)                                                                   tatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCa                                                                                   >  1:335675/1‑150 (MQ=255)                                                                   tatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCa                                                                                   <  1:79530/150‑1 (MQ=255)                                                                        cAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCg                                                                                                         <  1:749917/123‑1 (MQ=255)                                                                        cAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCg                                                                                                         >  2:749917/1‑123 (MQ=255)                                                                        cAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCg                                                                              >  1:677048/1‑150 (MQ=255)                                                                         aGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCgg                                                                             <  1:1001935/150‑1 (MQ=255)                                                                           cacaAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGa                                                                           <  1:1497246/150‑1 (MQ=255)                                                                           cacaAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGa                                                                           <  1:258255/150‑1 (MQ=255)                                                                              aaGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGatta                                                                        >  2:633170/1‑150 (MQ=255)                                                                              aaGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGatta                                                                        <  1:830862/150‑1 (MQ=255)                                                                               aGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGattat                                                                       <  1:68892/150‑1 (MQ=255)                                                                                   aaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAg                                                                   <  2:1610138/150‑1 (MQ=255)                                                                                     ggTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGtt                                                                 >  2:978837/1‑150 (MQ=255)                                                                                         ggTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacga                                                             >  2:995403/1‑150 (MQ=255)                                                                                         ggTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacga                                                             >  1:287506/1‑150 (MQ=255)                                                                                          gTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTacgac                                                            >  1:361090/1‑150 (MQ=255)                                                                                                aTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGcc                                                      <  2:619095/150‑1 (MQ=255)                                                                                                      gTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTAAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAg                                                <  2:670549/150‑1 (MQ=255)                                                                                                         ggCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaa                                             >  1:1074020/1‑150 (MQ=255)                                                                                                         ggCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaa                                             >  1:82073/1‑150 (MQ=255)                                                                                                          gCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaaa                                            >  2:462898/1‑150 (MQ=255)                                                                                                          gCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaaa                                            >  2:754856/1‑150 (MQ=255)                                                                                                          gCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGaaaa                                            >  2:474527/1‑150 (MQ=255)                                                                                                              tGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGc                                        >  1:1420961/1‑150 (MQ=255)                                                                                                                cTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                                      >  1:263407/1‑150 (MQ=255)                                                                                                                          ttttttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg                            <  1:1314137/150‑1 (MQ=255)                                                                                                                             tttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa                         <  1:754856/150‑1 (MQ=255)                                                                                                                               tATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                       <  1:758663/150‑1 (MQ=255)                                                                                                                               tATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGc                       <  2:1154178/150‑1 (MQ=255)                                                                                                                                 tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                     >  1:244792/1‑150 (MQ=255)                                                                                                                                 tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                     >  1:433572/1‑150 (MQ=255)                                                                                                                                 gccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                     >  1:1042077/3‑150 (MQ=255)                                                                                                                                 gccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                     >  1:65371/3‑150 (MQ=255)                                                                                                                                 gccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                     >  2:1593895/3‑150 (MQ=255)                                                                                                                                 gccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTAGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa                     >  2:1593901/3‑150 (MQ=255)                                                                                                                                  gCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGAt                    <  2:263407/150‑1 (MQ=255)                                                                                                                                   cGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATa                   <  1:1610578/150‑1 (MQ=255)                                                                                                                                         tgccctggCAAACCGGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca             <  1:196813/144‑1 (MQ=255)                                                                                                                                         tGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca             >  1:677477/1‑150 (MQ=255)                                                                                                                                         tGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATaca             >  1:852385/1‑150 (MQ=255)                                                                                                                                          gCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATacac            >  1:696062/1‑150 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca           >  2:716969/1‑150 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca           >  1:1478324/1‑150 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca           >  1:1472579/1‑150 (MQ=255)                                                                                                                                           ccAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATAcaca           >  1:1047887/1‑150 (MQ=255)                                                                                                                                                    gCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAg  >  1:670032/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                      ACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATAAATACACATTAAAATAG  >  NC_000913/715615‑715911

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