Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,083,929	C→T	P392S (CCG→TCG)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,083,929	0	C	T	98.2%	179.2 / ‑6.3	56	P392S (CCG→TCG)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (28/27);  minor base A (0/1);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAG  >  NC_000913/2083690‑2084172                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    cGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGat                                                                                                                                                                                                                                          <  1:438660/251‑1 (MQ=255)  gAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATAGTCAGCATTGCCTCGATAACCTGAAAATTCTGAGTGAAAATCCGCAAGTGCGCGAAGAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGGAACCGAGCCGCGTAATTTATCGGCACTGGat                                                                                                                                                                                                                                          >  1:433432/1‑250 (MQ=255)          tacgcCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACtttt                                                                                                                                                                                                                                 >  1:132177/1‑251 (MQ=255)                    aGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAAc                                                                                                                                                                                                                       <  1:215631/251‑1 (MQ=255)                         cataccgacTGGGAATTAATCGTCAGCATCGCCTCGAGTACCGGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGAGAACGTGGATGCACAGCTTTATACCGGCTGTTTCAGTGTCGCAGATCGTTCAGCAATGAAAATTGTGCTGTAAACCGAGCCGCGTAATTTATCGGCCCTGGATATCACTTTTGTTGATAAACGGAtt                                                                                                                                                                                                                  <  2:542406/248‑1 (MQ=255)                                  tGGGAATTAATCGTCAGCATTGCCTCGTTAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATACCGGCTTTTTCAGTGACGCAGCTCGTGCAGCAATGAAAAATGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtg                                                                                                                                                                                                         >  2:348920/1‑251 (MQ=255)                                   gggAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgt                                                                                                                                                                                                        >  1:546674/1‑251 (MQ=255)                                     gAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgttg                                                                                                                                                                                                      <  2:114668/251‑1 (MQ=255)                                     gAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgttg                                                                                                                                                                                                      <  1:81904/251‑1 (MQ=255)                                      aaTTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCtgttgt                                                                                                                                                                                                     >  1:126949/1‑251 (MQ=255)                                         tAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATAGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGCTTGAAAAGCTGTTGTTc                                                                                                                                                                                                   >  1:275649/1‑250 (MQ=255)                                         tAATCGGCAGCATTGCCTCGATAACCTGAAACTTCTGCGTGAAAATCCGCAAGTGCGCGACAAAGTGGTGGCGATATTCGCGGCAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTc                                                                                                                                                                                                   >  2:494831/1‑250 (MQ=255)                                             cGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAAtt                                                                                                                                                                                               <  1:402835/250‑1 (MQ=255)                                                 aGCATTGCCTCGATAACCTGCAAATTCTGCGTGAAAATCCGCACGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCgg                                                                                                                                                                                          >  1:130997/1‑251 (MQ=255)                                                      tGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCAc                                                                                                                                                                                      <  1:570529/250‑1 (MQ=255)                                                      tGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTGATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCAc                                                                                                                                                                                      <  1:296497/250‑1 (MQ=255)                                                       gCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGc                                                                                                                                                                                    >  1:471118/1‑251 (MQ=255)                                                            gATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACtt                                                                                                                                                                               >  1:497951/1‑251 (MQ=255)                                                               atccTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCTGCAATGAAAATTGTGCTGGAAACCGAGCCGCTTAATTTATCGGCACTGGATATCACTTTGGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTcc                                                                                                                                                                             <  1:502285/248‑1 (MQ=255)                                                                     aaaaTTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGAc                                                                                                                                                                      >  1:508180/1‑251 (MQ=255)                                                                     aaaaTTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGAc                                                                                                                                                                      <  1:507780/251‑1 (MQ=255)                                                                                     aTCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCg                                                                                                                                                        >  1:99500/1‑249 (MQ=255)                                                                                       ccGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGagc                                                                                                                                                     >  1:64892/1‑250 (MQ=255)                                                                                             gTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACg                                                                                                                                              <  2:126949/251‑1 (MQ=255)                                                                                                cgcgAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACgctg                                                                                                                                           >  1:20959/1‑251 (MQ=255)                                                                                                 gcgAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACgctgg                                                                                                                                          >  1:137649/1‑251 (MQ=255)                                                                                                    aaaaaGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACgctggctg                                                                                                                                       <  2:433432/251‑1 (MQ=255)                                                                                                            gtggCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCg                                                                                                                               >  1:158674/1‑251 (MQ=255)                                                                                                                cGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGc                                                                                                                           >  2:304322/1‑251 (MQ=255)                                                                                                                           gAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCAc                                                                                                                <  2:5405/251‑1 (MQ=255)                                                                                                                                  aaCCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCag                                                                                                           >  1:325680/1‑249 (MQ=255)                                                                                                                                          aCGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGc                                                                                                  <  1:348920/250‑1 (MQ=255)                                                                                                                                                             aTGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTg                                                                              <  2:137649/251‑1 (MQ=255)                                                                                                                                                               gCACAGCTTTATAACTGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGTAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGTCACGCAACTTCCCTGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGTAGGGTTATGCTGATGAATTGCa                                                                            <  2:380646/251‑1 (MQ=255)                                                                                                                                                                cacaGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATTTCACTTTTGTTGATa                                                                                                                                                                                                                          <  1:394628/108‑1 (MQ=255)                                                                                                                                                                cacaGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATa                                                                                                                                                                                                                          >  2:394628/1‑108 (MQ=255)                                                                                                                                                                  caGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGc                                                                             >  1:551633/1‑247 (MQ=255)                                                                                                                                                                  caGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGAt                                                                         <  1:452657/251‑1 (MQ=255)                                                                                                                                                                   aGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATg                                                                        <  2:132177/251‑1 (MQ=255)                                                                                                                                                                     cattATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATGTATCGGCACTGGATATCTCTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGTTTATTCTGATGAATTGCAGATGCt                                                                      <  2:130997/249‑1 (MQ=255)                                                                                                                                                                                  ttttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCAc                                                                                                                <  1:79165/196‑1 (MQ=255)                                                                                                                                                                                  ttttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCAc                                                                                                                >  2:79165/1‑196 (MQ=255)                                                                                                                                                                                  ttttCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCGTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAata                                                          <  2:551633/250‑1 (MQ=255)                                                                                                                                                                                      cAGTGAAGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGCTGAATTGCAGATGCTGGTACAACAATATGCCg                                                     >  2:90964/1‑251 (MQ=255)                                                                                                                                                                                           aCGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGAGATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTGCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGGCGCGAGGTCTTCACGCCAGAGTTTGTGCAGGGTTATGCTGATGAATTGCAGATGCTGGGACAACAATATGCCGATGAc                                                <  1:246677/251‑1 (MQ=255)                                                                                                                                                                                                 aTCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAagag                                          <  1:304322/251‑1 (MQ=255)                                                                                                                                                                                                    gTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGGTGACAAAGAGaa                                        >  2:298065/1‑250 (MQ=255)                                                                                                                                                                                                          caATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACAGTCGCCAGGTCTTCACGCCCGAGTTTTTGCAGGGTTATGCTGGTGAATTGCGGATGCTGGTCAAACAATATGCCGATGCCACAGAGAAAGGGgcg                                 >  1:191040/1‑251 (MQ=255)                                                                                                                                                                                                            aTGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTCTCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGCGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCCGGGTTATGCTGATGAATTGCAGCTGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCt                               >  2:292193/1‑251 (MQ=255)                                                                                                                                                                                                                     ttgctggaaagcgatgctggTAATTTAACGGGACTGGATAGCTATTTTGTTGCTAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGAACGCAACGTCCCGGGGACGCTTGATTATGCCGAGCAGCAACGCTGGTTGGAGCACCGTCGCCAGGTCTTGACGCCAGAGTTTT‑GCAGTGTTATGCTGATGATTTGCATATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGc                      <  1:83990/231‑1 (MQ=255)                                                                                                                                                                                                                           gAAACCGAGCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCCGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTg                >  1:571468/1‑251 (MQ=255)                                                                                                                                                                                                                                  aGCCGCGTAATTTATCGGCAATGGATATCCCTTTTGTTGATAAACGGCTTGAAAAGCTGTTGTTCAATTATCGGGCACGCAAATTCCCGGGGCCGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGCATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTAc         >  2:400481/1‑251 (MQ=255)                                                                                                                                                                                                                                   gCCGCGTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTAcg        <  2:99500/251‑1 (MQ=255)                                                                                                                                                                                                                                     cgcgTAATTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTAcgcg      >  1:317001/1‑251 (MQ=255)                                                                                                                                                                                                                                       cgTAATTTGTCGGAACGGGATATCGCTTTTGTTGATTGTAGGATTGATAAGCTGTTGTCCAATAATCGGGCACGCGACGTCCCCGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGGGCACCGTCGGCAGGTCTTCACGTCATAGTTTGTGCAGGGTTATGATGATGAATTGCTGATGCTGGTACAACAATATTGCGATAACAAAGAGAAAGTGGCGCTGGTAATAGCACTTTGGCAGttaaccg     <  2:463270/250‑7 (MQ=255)                                                                                                                                                                                                                                          aaTTTATCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAATCTGTTGTTCAATTAGCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGaac  <  2:312095/250‑2 (MQ=255)                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                    CGAATACGCTACGCCCGGAAGATGCCGACCGACTGGGAATTAATCGTCAGCATTGCCTCGATAACCTGAAAATTCTGCGTGAAAATCCGCAAGTGCGCGAAAAAGTGGTGGCGATATTCGCGGAAGCCGAACCGTTTACGCCTTCAGATAACGTGGATGCACAGCTTTATAACGGCTTTTTCAGTGACGCAGATCGTGCAGCAATGAAAATTGTGCTGGAAACCGAGCCGCGTAATTTACCGGCACTGGATATCACTTTTGTTGATAAACGGATTGAAAAGCTGTTGTTCAATTATCGGGCACGCAACTTCCCGGGGACGCTGGATTATGCCGAGCAGCAACGCTGGCTGGAGCACCGTCGCCAGGTCTTCACGCCAGAGTTTTTGCAGGGTTATGCTGATGAATTGCAGATGCTGGTACAACAATATGCCGATGACAAAGAGAAAGTGGCGCTGTTAAAAGCACTTTGGCAGTACGCGGAAG  >  NC_000913/2083690‑2084172

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