Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,987,860	A→G	12.0%	intergenic (‑102/+13)	yecJ ← / ← azuC	DUF2766 domain‑containing protein YecJ/uncharacterized protein AzuC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,987,860	0	A	G	12.0%	238.5 / 20.0	92	intergenic (‑102/+13)	yecJ/azuC	DUF2766 domain‑containing protein YecJ/uncharacterized protein AzuC
Reads supporting (aligned to +/- strand):  ref base A (39/41);  new base G (6/5);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTAT  >  NC_000913/1987713‑1988009                                                                                                                                                    |                                                                                                                                                       ccACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGc                                                                                                                                                      <  1:115681/150‑1 (MQ=255)  cacaTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCa                                                                                                                                                     <  2:1020220/150‑1 (MQ=255)        aGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCGGCAGG‑ttttt                                                                                                                                             >  1:525018/1‑150 (MQ=255)          aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGTCAGG‑tttttt                                                                                                                                            <  1:982918/150‑1 (MQ=255)          aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑tttttt                                                                                                                                            >  1:918708/1‑150 (MQ=255)           aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGTACACAAAAAAACCTGGCGGAGCAGG‑TTTTTTg                                                                                                                                           <  2:393111/150‑1 (MQ=255)           aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGG‑ttttttt                                                                                                                                           <  2:1396831/150‑2 (MQ=255)           aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTg                                                                                                                                           <  1:1057989/150‑1 (MQ=255)            aCCAATTCGTGATCGGCATTCAGCGGCTGTGACATCGTAAACTCCTGATAGCATTTTGAAAGCAGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCTTCATCTGTACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGt                                                                                                                                          <  2:1416767/150‑1 (MQ=255)             ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCGGCAGG‑tttttttgtg                                                                                                                                        <  1:909056/150‑5 (MQ=255)              cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGttt                                                                                                                                        >  2:1152784/1‑149 (MQ=255)              cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTa                                                                                                                                        >  1:992671/1‑150 (MQ=255)              cAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTa                                                                                                                                        >  1:1258956/1‑150 (MQ=255)                aTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATc                                                                                                                                      <  2:244742/150‑1 (MQ=255)                aTTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATc                                                                                                                                      <  2:581164/150‑1 (MQ=255)                 ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCg                                                                                                                                     <  2:1189784/150‑1 (MQ=255)                   ccTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGG‑tttttttg                                                                                                                                          >  2:1393779/1‑141 (MQ=255)                   ccTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGG‑tttttttg                                                                                                                                          <  1:1393779/143‑3 (MQ=255)                   ccTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGa                                                                                                                                   >  1:667971/1‑150 (MQ=255)                   ccTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGa                                                                                                                                   >  2:583053/1‑150 (MQ=255)                    cTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGaa                                                                                                                                  <  1:1359852/150‑1 (MQ=255)                     tGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAAc                                                                                                                                 >  1:1104113/1‑150 (MQ=255)                      gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACa                                                                                                                                >  1:1270564/1‑150 (MQ=255)                      gATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACa                                                                                                                                <  1:513770/150‑1 (MQ=255)                          ggCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGG‑tttttttg                                                                                                                                          >  1:1567379/1‑134 (MQ=255)                          ggCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCTCCGGCAGG‑tttttttg                                                                                                                                          <  2:1567379/136‑3 (MQ=255)                           gCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTg                                                                                                                           >  1:117721/1‑150 (MQ=255)                             aTTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGcc                                                                                                                         >  1:1376325/1‑150 (MQ=255)                             aTTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGcc                                                                                                                         <  2:287492/150‑1 (MQ=255)                             aTTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGcc                                                                                                                         >  1:450999/1‑150 (MQ=255)                                 aGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGc                                                                                                                     <  2:926942/150‑1 (MQ=255)                                  gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCg                                                                                                                    <  2:1217713/150‑1 (MQ=255)                                  gCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCg                                                                                                                    <  2:918708/150‑1 (MQ=255)                                   cGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCgg                                                                                                                   <  1:583053/150‑1 (MQ=255)                                    ggCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGt                                                                                                                  <  1:1421491/150‑1 (MQ=255)                                        gggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGt                                                                                                              >  2:140227/1‑150 (MQ=255)                                           aCATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCtt                                                                                                           <  1:183855/150‑1 (MQ=255)                                            cATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                                          <  1:654770/150‑1 (MQ=255)                                            cATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                                          >  2:1591141/1‑150 (MQ=255)                                             aTCGTAAACTCCTTATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTg                                                                                                         <  2:365713/150‑1 (MQ=255)                                               cGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGaa                                                                                                       >  2:1331436/1‑150 (MQ=255)                                                     cTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCtt                                                                                                 >  1:1141813/1‑150 (MQ=255)                                                     cTCCTGATACCATTTTGAAAGCGGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGTAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCtt                                                                                                 <  2:342838/150‑1 (MQ=255)                                                       ccTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                               >  1:1094701/1‑150 (MQ=255)                                                       ccTGATAGCATTTTGAAAGCCGTTATAGTAGAGACATCACATCTTCAGCGATAGTCACATCCACCGTCATCATGACACAAAAAAACCTTCCGGCGCATG‑TTTTTTGTTATCGGAACAAATTGCCTGGCGGTACTTCTTTGAACGTCTTGc                                                                                               <  2:1578725/150‑1 (MQ=255)                                                         tGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACa                                                                                                                                <  2:106096/115‑1 (MQ=255)                                                         tGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACa                                                                                                                                >  1:106096/1‑115 (MQ=255)                                                           aTAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATa                                                                                           >  1:1471513/1‑150 (MQ=255)                                                              gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                               >  1:125311/1‑143 (MQ=255)                                                              gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                               <  2:131212/143‑1 (MQ=255)                                                              gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                               <  2:125311/143‑1 (MQ=255)                                                              acATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGc                                                                                               >  1:131212/2‑143 (MQ=255)                                                                aTTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAt                                                                                      >  1:1237930/1‑150 (MQ=255)                                                                   ttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACGGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGa                                                                                   <  2:979535/150‑1 (MQ=255)                                                                   ttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGa                                                                                   >  2:644289/1‑150 (MQ=255)                                                                   ttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGa                                                                                   >  2:300659/1‑150 (MQ=255)                                                                       aaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAAGCTGGCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTGTTGCAATAGTTATTGAACAt                                                                               <  2:220060/150‑1 (MQ=255)                                                                          ccGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACt                                                                            >  2:1325351/1‑150 (MQ=255)                                                                                  gtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGa                                                                    >  2:220804/1‑150 (MQ=255)                                                                                   tagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAt                                                                   >  1:1194628/1‑150 (MQ=255)                                                                                      cGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtttt                                                                <  2:1471513/150‑1 (MQ=255)                                                                                       gACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTg                                                               <  2:1598843/150‑1 (MQ=255)                                                                                             acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGt                                                         <  2:891679/150‑1 (MQ=255)                                                                                                 cTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCa                                                     <  2:667971/150‑1 (MQ=255)                                                                                                     aGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGGCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATcgc                                                 <  1:780941/150‑1 (MQ=255)                                                                                                      gCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATcgcg                                                >  1:388224/1‑150 (MQ=255)                                                                                                           aGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACt                                           >  1:583145/1‑150 (MQ=255)                                                                                                            gTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTc                                          >  1:225489/1‑150 (MQ=255)                                                                                                              cacaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCg                                        >  1:1027382/1‑150 (MQ=255)                                                                                                              cacaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCg                                        >  1:163590/1‑150 (MQ=255)                                                                                                              cacaTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCg                                        >  2:1536739/1‑150 (MQ=255)                                                                                                                caTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                      >  2:262996/1‑150 (MQ=255)                                                                                                                caTCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                                      <  2:1376325/150‑1 (MQ=255)                                                                                                                       cGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGt                               <  1:1425964/150‑1 (MQ=255)                                                                                                                        gTCATCAGGACACAAAACAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTTGCGGGACGTCTTTGAACTTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGTCACTCCTACCATTTGtt                              <  2:10759/150‑1 (MQ=255)                                                                                                                         tcatcaGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTta                             <  1:220804/150‑1 (MQ=255)                                                                                                                          catcaGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTtat                            <  2:451931/150‑1 (MQ=255)                                                                                                                           atcaGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTtata                           >  1:1423070/1‑150 (MQ=255)                                                                                                                             caGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACa                         <  2:1258956/150‑1 (MQ=255)                                                                                                                               ggACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAgg                       <  1:1591141/150‑1 (MQ=255)                                                                                                                                 acacAAAAAAACCTGCCGGAGGAGG‑TTTTTTGTTATCGGAACATATAGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtg                     <  2:1613558/150‑1 (MQ=255)                                                                                                                                 acacAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtg                     >  1:526532/1‑150 (MQ=255)                                                                                                                                  cacaAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAAGATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgtgt                    <  2:1094701/150‑1 (MQ=255)                                                                                                                                    caAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa                  >  2:635400/1‑150 (MQ=255)                                                                                                                                    caAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa                  <  1:300659/150‑1 (MQ=255)                                                                                                                                     caaaaaaCCTGCTCCGGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt                 >  1:262548/2‑150 (MQ=255)                                                                                                                                     aaaaaaaCCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt                 >  1:925436/1‑150 (MQ=255)                                                                                                                                     aaaaaaaCCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt                 >  1:931189/1‑150 (MQ=255)                                                                                                                                      aaaaaaCCTGCTCCGGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt                >  2:65494/1‑150 (MQ=255)                                                                                                                                       aaaaaCCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg               <  2:1270564/150‑1 (MQ=255)                                                                                                                                        aaaaCCTGCTCCGGCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg               >  1:908134/1‑150 (MQ=255)                                                                                                                                            ccTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCttt          >  1:1495711/1‑150 (MQ=255)                                                                                                                                                    aGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTtat  <  2:450999/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       CCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGCCGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTGCGCTTAT  >  NC_000913/1987713‑1988009

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: ‑