Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	83,344	G→A	A62A (GCC→GCT)	leuA ←	2‑isopropylmalate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	83,344	0	G	A	95.1%	105.3 / NA	43	A62A (GCC→GCT)	leuA	2‑isopropylmalate synthase
Reads supporting (aligned to +/- strand):  ref base G (0/2);  major base A (21/18);  minor base C (2/0);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

CGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAT  >  NC_000913/83109‑83589                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      cGCAAGAAAATTCAACATCAGCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCAACCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCAc                                                                                                                                                                                                                                         <  2:90083/250‑1 (MQ=255)             cAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATca                                                                                                                                                                                                                            <  1:60044/251‑2 (MQ=255)                            aaTTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACg                                                                                                                                                                                                             >  1:10354/1‑251 (MQ=255)                               tACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGATTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgggg                                                                                                                                                                                                          >  1:453250/1‑251 (MQ=255)                                 cGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGATTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGa                                                                                                                                                                                                         >  1:43203/1‑250 (MQ=255)                                      gcgTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACAGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGCGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCg                                                                                                                                                                                                    >  2:50487/1‑250 (MQ=255)                                        gTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGCGCGCTGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACt                                                                                                                                                                                                 >  2:237427/1‑251 (MQ=255)                                        gTTTCACCATATAGATAGCGCGTTCGATCACTTCATCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACt                                                                                                                                                                                                 <  2:201533/251‑1 (MQ=255)                                               cATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCa                                                                                                                                                                                          <  1:38693/251‑1 (MQ=255)                                                                cgatcccttcatCCAGCGTGCGGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACg                                                                                                                                                                          <  2:139856/239‑1 (MQ=255)                                                                   tcacttcatCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTc                                                                                                                                                                       >  1:361407/9‑250 (MQ=255)                                                                    cacttcatCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCa                                                                                                                                                                      >  1:153038/8‑250 (MQ=255)                                                                     acttcatCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCa                                                                                                                                                                      >  1:298776/7‑249 (MQ=255)                                                                          atCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGcc                                                                                                                                                                >  1:376253/2‑250 (MQ=255)                                                                                            cTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGCTATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCCCCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGttttt                                                                                                                                             >  2:215394/1‑251 (MQ=255)                                                                                                gtggCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCGTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTtttcttt                                                                                                                                         <  2:222038/251‑1 (MQ=255)                                                                                                     gATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGCTATCTTTATCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCacac                                                                                                                                    >  2:304303/1‑251 (MQ=255)                                                                                                      aTGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACt                                                                                                                                   <  1:40990/251‑1 (MQ=255)                                                                                                          gCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCaa                                                                                                                                <  1:9677/250‑1 (MQ=255)                                                                                                              tGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCtt                                                                                                                            <  2:210016/250‑1 (MQ=255)                                                                                                                          aaTAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCt                                                                                                               >  1:415962/1‑251 (MQ=255)                                                                                                                                   gtGAATACGGAAGGCTTCGGCTAATTTCAGGGAGTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCg                                                                                                       <  2:275203/249‑1 (MQ=255)                                                                                                                                     gAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGGCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTcg                                                                                                    <  2:376253/251‑1 (MQ=255)                                                                                                                                                  ttCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTAt                                                                                       >  2:368047/1‑251 (MQ=255)                                                                                                                                                    cGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgg                                                                                                                                                                                                            <  1:453110/132‑1 (MQ=255)                                                                                                                                                    cGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCg                                                                                     <  1:37575/251‑1 (MQ=255)                                                                                                                                                    cGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATATTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACgg                                                                                                                                                                                                            >  2:453110/1‑132 (MQ=255)                                                                                                                                                        gACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTtttct                                                                                                                                           <  1:326480/193‑1 (MQ=255)                                                                                                                                                        aaCTTTCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCCGTTtttct                                                                                                                                           >  2:326480/2‑193 (MQ=255)                                                                                                                                                            ttCAGGGATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAaataat                                                                             >  2:396098/1‑251 (MQ=255)                                                                                                                                                              cAGGGATCCGGCCGCCTCGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGa                                                                           <  2:88233/251‑1 (MQ=255)                                                                                                                                                                 ggATTCGGCCGCCACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGTATTCGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGTGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTATCGCCTGTTCATCGTCGCGCAATGTGGTATCGAAAATAATGACtt                                                                        <  1:395422/251‑1 (MQ=255)                                                                                                                                                                           gccACGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCAt                                                              <  1:445738/251‑1 (MQ=255)                                                                                                                                                                               cGTCGATATCTTTCTCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGtt                                                          <  1:368047/251‑1 (MQ=255)                                                                                                                                                                                     tatCTTTCTCCACGAAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGATCATGGTTTGGGTc                                                    >  1:378624/1‑251 (MQ=255)                                                                                                                                                                                         tttctCCACGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCACACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCCAGCTTGCCTGGAACGCCTGTTCACCGTCGCGCCATGTTGTATCGACAATAATGTCATGCTGGCTCATGGTTTGGGTCCTTg                                                >  2:381252/5‑251 (MQ=255)                                                                                                                                                                                              ccACGCAGCGCGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCAGGAGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGATCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCACGCAATGTGGTATAGAAAATAATGTCTTGCTGGCGCATGGTTTGGGTCCTTGtctgt                                           >  1:466671/1‑249 (MQ=255)                                                                                                                                                                                                aCGCAGCGAGCTAACGCGCATACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCttt                                         <  1:237427/251‑1 (MQ=255)                                                                                                                                                                                                                  caTACGCGGCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCgg                       >  1:193128/1‑251 (MQ=255)                                                                                                                                                                                                                         ggCTGTTTTTAACCTGGCGAGCGATGGTTTGCACTGATTCAACATCACCCGGCGAAGAAACGGGGAAACAGACTTCCATCACGTCAACACCCATACGCTCAAGGGCAAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTACCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATaa                 >  1:435853/1‑250 (MQ=255)                                                                                                                                                                                                                                tttAACCTGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAAcc          <  1:372571/250‑1 (MQ=255)                                                                                                                                                                                                                                     ccTGGCGAGCGATGGTTTGCACTGATTCCAAATCACCCGGCGAAGAAACGGGGACACCGCCTTCCATCACGTAAACACCCATACGCTCAAGGACCAGCGCAATTTGCAGTTTTTCTTTCCCACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGTAAATATTGACTTGTTGGCTCATGGTTTGGGTCCTTGACTCTTTTAGAGCGACTAGATTCGGGCATaaaaaaatcagcgc     >  1:69104/1‑243 (MQ=255)                                                                                                                                                                                                                                       tGGCGAGCGATGGTTTGCACTGATTCAAAATCACCCGGCGAAGAAACGGGGAAACCGACTTCAATCACGTCAACACCCATACGCTCAAGGGCCAGCGAAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCACTGTGGTATCAACACTAATAACTTGCTGGATCATGGTTTGGGTCCATGTCTCTTTTAACGCGCCTCGCTTCGGGCATAAAAAACACCGCGCAAt  >  2:347416/1‑251 (MQ=255)                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                      CGCAAGAAAATTCAACATCATCGGTGTAATTACGGGCGCGTTTCACCATATAGATAGCGCGTTCGATCACCTCGTCCAGCGTGCTGCGCAGCTTGGTGGCGATGTGCATTGGCGAAGTGGCAATAAAGGTATGAATACGGAAGGCTTCGGCGACTTTCAGGGATTCGGCCGCCACGTCGATATCTTTTTCCACGCAGCGAGCTAACGCACATACGCGGCTGTTTTTAACCTGGCGGGCGATGGTTTGCACCGATTCAAAATCGCCCGGCGAAGAGACGGGGAAACCGACTTCCATCACGTCAACACCCATACGCTCAAGGGCCAGCGCAATTTGCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAAT  >  NC_000913/83109‑83589

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