Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	69,752	T→A	G99G (GGA→GGT)	araB ←	ribulokinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	69,752	0	T	A	98.0%	142.8 / NA	50	G99G (GGA→GGT)	araB	ribulokinase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base A (30/18);  minor base C (0/1);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCC  >  NC_000913/69506‑69992                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  gCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                                                                                                                                                                               >  1:120786/1‑251 (MQ=255)                 gATGCGGCAGATTGCGCCACGGCGCTGTCCTGTCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGc                                                                                                                                                                                                                               >  2:40920/1‑251 (MQ=255)                     cGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTCCGGGc                                                                                                                                                                                                                           >  2:334486/1‑251 (MQ=255)                        cAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTc                                                                                                                                                                                                                        >  1:169799/1‑251 (MQ=255)                           aTTGCGCCCCGGCGCGGCCCTGGCGAGTCAAATGAAGGTTTTTGGCCCAAAACCATTCGCTGAAAAAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAg                                                                                                                                                                                                                     <  1:11674/251‑1 (MQ=255)                                  cACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGTCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTAc                                                                                                                                                                                                              <  1:482452/251‑1 (MQ=255)                                    cGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTg                                                                                                                                                                                                            >  1:30256/1‑251 (MQ=255)                                       cgcTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTc                                                                                                                                                                                                          >  1:82449/1‑250 (MQ=255)                                                  cGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGAGCGGGCGTCGAGCCGGTACTGTCAACGCCAATccc                                                                                                                                                                                              >  1:55103/1‑251 (MQ=255)                                                                gATTTTTGCCCAGAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTgcgc                                                                                                                                                                               >  2:114230/1‑251 (MQ=255)                                                                       gCCCAGGACCATTCGCTGGAATCAATACCGACACTATAGCGGGAGTAGTCAACATTGCACGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCCTCT‑TTC‑CACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTATCGGCAAACTACGGGCGCAGCGCCAGCAAGTTACCGTCGGTTTCAATCGGTGAGGGCACCGAGCCGGGCCTGTCAACGCCAATCCCG‑ACCCCAGCTGCGCGCTGttt                                                                                                                                                                        >  2:561342/1‑250 (MQ=255)                                                                            gAACCATTCGCTGGAATAAATACCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGc                                                                                                                                                                   >  1:591652/1‑251 (MQ=255)                                                                             acccttccgcggaaaaaaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCt                                                                                                                                                                  <  1:377896/235‑1 (MQ=255)                                                                                    cGCTGGAATAAATACCGCCAATATAGCGGGCGTAGTCAACATTGCCCGGCGCGCGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGCACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCCGCGTCAGCACGTTACCGTCGGCATCAATCGGTGCGGGTGTCGAGCAGGTACTGGCAACGCCAATCCCA‑ACCACGGCTGCGCGCTCTTCCACGCTAAActct                                                                                                                                                           >  2:240288/1‑251 (MQ=255)                                                                                          ataatatcccgccaATGTAGCGGGAGTAGTCAACATTGCCCGTCGTGTGGCACAACCGGGTAATCTCTTC‑GCTTCT‑TGC‑AACCGCAGTGTGGTCTTTCCACAATCCTAACATCGCGTGCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACTCTAAGCTCTGCAAGc                                                                                                                                                     <  2:254461/240‑1 (MQ=255)                                                                                            taaataCCGCCCATATAGCGGGAGTAGTCAACATTGCCCAGCGCGTGGCCCAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTCCTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGCCGCTAAGCTCTGccagacc                                                                                                                                                   >  2:353139/1‑245 (MQ=255)                                                                                             aaataCCGCCAATATAGCGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGt                                                                                                                                                                                                                                               >  2:235697/1‑159 (MQ=255)                                                                                              aaaacccccaATATAGCGGGAGTAGGCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCCCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAACCGGGGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACgg                                                                                                                                                 <  2:185512/244‑1 (MQ=255)                                                                                                              cGGGAGTAGTCAACATTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTc                                                                                                                                 >  1:151853/1‑251 (MQ=255)                                                                                                                         aacattGCCCGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCa                                                                                                                       >  1:552361/5‑250 (MQ=255)                                                                                                                                cccGGCGCGTGGCACAAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTc                                                                                                               >  1:279913/1‑251 (MQ=255)                                                                                                                                       cgTGGCACAAACGGGCAATCTATTC‑GCTGCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTGGGCAAACTCCGGGCGCAGCGACAGCACGTTACCGTCGGCATCAATCGCCGCGGGCGTCGGGCGGGTACTGTCAACGCCAATCCCGCACC‑CAGCTGCGCGCTGTTCGACGCTACGATCTGCAAGCACGGTTTTCAGAGCCGCTTACATTGCCTCATTGTAGTccggcgg                                                                                                         >  2:176268/1‑244 (MQ=255)                                                                                                                                        gTGGCACATACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCACACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGCGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCCCGGTTTTCAGTGCCGCTGCCATTGCCTCAATGTAGTTACGCGgat                                                                                                       >  2:471955/1‑251 (MQ=255)                                                                                                                                            cacaAACGGGTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGatgatg                                                                                                   <  2:526489/251‑1 (MQ=255)                                                                                                                                                   gggTAATCTCT‑GCGCTCCTGT‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACGCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAAc                                                                                            <  1:629945/251‑1 (MQ=255)                                                                                                                                                    ggTAATCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGAGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGTGCGTCGAGCCGGTACTGTCAACGCCAATCTCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGAATTATGACGGAACt                                                                                           >  1:167012/1‑251 (MQ=255)                                                                                                                                                       aaTCTCTTC‑GCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGt                                                                                        <  1:114230/251‑1 (MQ=255)                                                                                                                                                       aaTCTCTTC‑GCTTCT‑TGC‑AACCGCAGTGTGGTCTTTCCACAATACGAACCTCGCGTTCGGGTTTTCGGCAACCTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGt                                                                                        <  1:56896/251‑1 (MQ=255)                                                                                                                                                            tcttcGCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAtt                                                                                    >  1:386380/5‑251 (MQ=255)                                                                                                                                                             cttcccttct‑t‑CCAACCGCAGTTTGGTCTTTCCACAATACGACCACCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCCCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGtt                                                                                       <  2:410978/242‑1 (MQ=255)                                                                                                                                                               tcGCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCg                                                                                  <  1:587025/249‑1 (MQ=255)                                                                                                                                                                cGCTTCT‑TTC‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCggg                                                                                <  1:498320/251‑1 (MQ=255)                                                                                                                                                                     tt‑t‑c‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATc                                                                           >  1:400173/2‑250 (MQ=255)                                                                                                                                                                     tt‑t‑c‑AACCGCAGGGTGGTCTTTCCAGAATACGAACATCGCGTTCGGGTTTTCGACATACTCCGGGCGCAGCGCCAGCACTTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTAGTCGGGGCATc                                                                           <  1:482143/249‑1 (MQ=255)                                                                                                                                                                      cagattcaCCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCAGTTACCGGCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGTCAATCCCG‑CCCACAGCTGCGCGCTTTTCGACGATAAGCTCTGCAAGCCCGGTTTTCAGTGCCGCTTCCATTGACTCTATGTAGTCCCTCGGAGGTTGACGAACCTGGTTATTGGGGGCAt                                                                            >  2:624491/8‑251 (MQ=255)                                                                                                                                                                          c‑AACCGCAGTGTGGTCTTTCCACAATACGCACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCCTCAATCGGTGCGGGCGTCGAGCCGGTCCTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGCTGATGACGGAACTGGTTATTCGGGGCATcacc                                                                        >  2:697414/1‑249 (MQ=255)                                                                                                                                                                                       ggTCTTTCCACAATACGAACATCGCTGTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCtttt                                                           <  2:279913/251‑2 (MQ=255)                                                                                                                                                                                       ggTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCtttt                                                           >  1:249457/1‑250 (MQ=255)                                                                                                                                                                                        gTCTTTCCACAATACGACCATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCCCGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCttttt                                                          <  2:591652/251‑3 (MQ=255)                                                                                                                                                                                          cTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGACAACGCCAATCCAG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCCTTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGc                                                        >  2:599957/1‑251 (MQ=255)                                                                                                                                                                                              ccctcatacgagcgtggtgTTCGGGTTCTCGGCATGCGCCGGCCGCAGCACCAGCCCGTTCCCGTCGGCATCAACCGGTGCGGGCGTCGAGCCGGTTCTGTCAACGCCAATCCCG‑ACCCCAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTTACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATGCGGGGCAGCGCAAAATTGCCCTTGTTGCCAAc                                                    <  2:486132/233‑1 (MQ=255)                                                                                                                                                                                               cacaATACGAACATCGCGTTCGGGTTTTCGGCAAACCCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACg                                                   >  1:589552/1‑251 (MQ=255)                                                                                                                                                                                                  aaTACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACgggg                                                <  1:158367/251‑1 (MQ=255)                                                                                                                                                                                                      cGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGGACTGGGTATTCGGGGAATCACAGAATTGCCCTTTTTGCCAACGGGGATAc                                            >  2:489844/1‑251 (MQ=255)                                                                                                                                                                                                                                   cctgggCGCGGCGCCGGCACGGTACCGTCGGCGTCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCGCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTGCGCTGGTGGCGATCTCTTCAcc               <  2:223932/248‑1 (MQ=255)                                                                                                                                                                                                                                    ccGGGCGCAGCGCCAGCACGGTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCAcc               <  2:386380/250‑1 (MQ=255)                                                                                                                                                                                                                                         cgcAGCGCCAGCACGTTACCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTAGCAAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCgatggc         >  1:516151/1‑246 (MQ=255)                                                                                                                                                                                                                                            aGCGCCAGCACGTTACCGCCGGCATCAATACGTGCGGACATCGCGCCGGTCCTGTCAACCCCCATCCCA‑CCCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGACGCTTCAATTGACTAAATGTAGTCACGCGGATGATGACGGACCTGGTTATTCGGGGCATCACAAAATTGCCCTCTTTGCAAACGGGGATACAACTCTACTCTGGTGGCGATCTCTTCAccactggcgcc      >  2:468087/1‑242 (MQ=255)                                                                                                                                                                                                                                               gcCAGCACGTTACCGTCGGCAGCAATCGGTGCGGGCGTCGAGCCGTGACGGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGCTGGCGCAGTc   <  2:400173/251‑1 (MQ=255)                                                                                                                                                                                                                                                ccAGAACGTTACCGTCGGCATCAATCGGTTCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCAGGGATTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAACATTGCCCTTTTTGCCAACGGGGATACCACTCTACGCTGGTGGCGATATCTTCACCGCTAGCGAAGtac  >  2:354895/1‑249 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                  GCACAGCTCAATCCACGATGCGGCAGATTGCGCCACGGCGCTGTCCTGGCGAGTCACATGCAGGATTTTTGCCCAGAACCATTCGCTGGAATAAATACCACCAATGTAGCGGGAGTAGTCAACGTTGCCCGGCGCGTGGCACAAACGGGTAATCTCTTCCGCTTCT‑T‑C‑AACCGCAGTGTGGTCTTTCCACAATACGAACATCGCGTTCGGGTTTTCGGCAAACTCCGGGCGCAGCGCCAGCACGTTTCCGTCGGCATCAATCGGTGCGGGCGTCGAGCCGGTACTGTCAACGCCAATCCCG‑ACCACAGCTGCGCGCTGTTCGACGCTAAGCTCTGCAAGCACGGTTTTCAGTGCCGCTTCCATTGACTCAATGTAGTCACGCGGATGATGACGGAACTGGTTATTCGGGGCATCACAAAATTGCCCTTTCTGCCAACGGGGATACCACTCTACGCTGGTGGCGATCTCTTCACCGGTAGCGCAGTCC  >  NC_000913/69506‑69992

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