Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,914,928	C→T	A203T (GCC→ACC)	proQ ←	RNA chaperone ProQ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,914,928	0	C	T	97.8%	133.5 / ‑3.9	45	A203T (GCC→ACC)	proQ	RNA chaperone ProQ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (21/23);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-01

ACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTC  >  NC_000913/1914686‑1915168                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 aCCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGCTCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATcg                                                                                                                                                                                                                                          >  2:153400/1‑251 (MQ=255)              gcGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCg                                                                                                                                                                                                                              <  1:320929/250‑1 (MQ=255)                          ttAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTc                                                                                                                                                                                                                   <  2:98498/249‑1 (MQ=255)                               aCCGGAATTTGAGCAGCACGCGGTATATCTTCTACAGCGAAGGTCTGTCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTATAAAACATGTTCATGCCTTGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCTAAGTCAGCTGGACGCGGATGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCTCTTTCACCTTCAGGGCtt                                                                                                                                                                                                            <  1:21230/251‑1 (MQ=255)                                                  gcgTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGa                                                                                                                                                                                         <  1:46692/251‑1 (MQ=255)                                                  gcgTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGa                                                                                                                                                                                         <  1:74860/251‑1 (MQ=255)                                                        tatCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTc                                                                                                                                                                                   >  1:181822/1‑251 (MQ=255)                                                             tctACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGaa                                                                                                                                                                               <  2:136176/250‑1 (MQ=255)                                                                       aGGTCTGGCCTGCTAGTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGtgtg                                                                                                                                                                    <  1:350331/251‑1 (MQ=255)                                                                            tGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGtt                                                                                                                                                               >  2:456691/1‑251 (MQ=255)                                                                            tGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGtt                                                                                                                                                               >  1:478671/1‑251 (MQ=255)                                                                              gCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAAGTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGttct                                                                                                                                                             <  2:427157/251‑1 (MQ=255)                                                                               ccTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCa                                                                                                                                                                                                                        >  1:182382/1‑191 (MQ=255)                                                                               ccTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCa                                                                                                                                                                                                                        <  2:182382/191‑1 (MQ=255)                                                                                          aaGCAGGCCAGCGAACGCGGTAAGCCTAAAAAACATGTTCAAGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAAGTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGGGCTGTTCTTCGCGAGGTGCt                                                                                                                                                 <  2:187373/251‑1 (MQ=255)                                                                                                 ccAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACTCCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTg                                                                                                                                          <  1:179164/251‑1 (MQ=255)                                                                                                     cTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGtttt                                                                                                                                      >  1:123898/1‑251 (MQ=255)                                                                                                      tAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCAGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGTACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGATTTGACCCGCTTTCACCTTAAGGGCTTGTCGGACAGTCAGAGCTGCAATGGCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGAGCTTTTACTGttttt                                                                                                                                     >  2:179546/1‑251 (MQ=255)                                                                                                               aaGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTtc                                                                                                                            <  1:64360/251‑1 (MQ=255)                                                                                                                     aaaaaaCATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGACGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGATTTTGGCGCTTTCTCTAcc                                                                                                                      >  1:15894/1‑251 (MQ=255)                                                                                                                              gTTCATGCCTGGCCCGGCCTCCGTTTCAGAACTCCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACGGTTTTTGGCGCTTTCTCTACCGGCTTTTg                                                                                                              <  1:411033/250‑1 (MQ=255)                                                                                                                                        ggCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAgg                                                                                                     <  1:145776/249‑1 (MQ=255)                                                                                                                                          cccGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGttt                                                                                                  >  1:156539/1‑250 (MQ=255)                                                                                                                                             ggCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGt                                                                                                                                         >  1:369918/1‑208 (MQ=255)                                                                                                                                             ggCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGt                                                                                                                                         <  2:369918/208‑1 (MQ=255)                                                                                                                                               ccTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAAACGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTCCGTTc                                                                                            >  1:465367/1‑251 (MQ=255)                                                                                                                                                   cGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAgcg                                                                                        >  2:327843/1‑251 (MQ=255)                                                                                                                                                                    tgtTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCggcgtgg                                                                       >  1:188850/1‑251 (MQ=255)                                                                                                                                                                    tgtTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCggcgtgg                                                                       >  1:218062/1‑251 (MQ=255)                                                                                                                                                                    tgtTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCggcgtgg                                                                       >  2:146396/1‑251 (MQ=255)                                                                                                                                                                                acaATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACg                                                           <  2:149533/251‑1 (MQ=255)                                                                                                                                                                                        aGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTGCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGc                                                   <  2:413147/251‑1 (MQ=255)                                                                                                                                                                                           cATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTCGACACGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGATGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGAGtt                                                >  1:119895/1‑251 (MQ=255)                                                                                                                                                                                              aCCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTCTCCGTCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGCAGTGTCAGAAACCGGGGTGTTCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCCCGAGGTTTCCGTTCAGCGCCTTCTTTGCGGCGTGGCATAGTCGGACGTGGCTTGCGttatg                                             >  2:137290/1‑248 (MQ=255)                                                                                                                                                                                                  gAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcg                                         >  1:219724/1‑251 (MQ=255)                                                                                                                                                                                                  gAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTcgcggcg                                         <  2:218062/251‑1 (MQ=255)                                                                                                                                                                                                    aTTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTGGGACGTGGCTTGCGTTcgcggcgcg                                       >  2:244398/1‑251 (MQ=255)                                                                                                                                                                                                            ggacgcggacgcCGTCTTTTGTGATTTCGAATACGGTGGTATCCATCGCGTTTGGACCCGCTTTCCCCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTTTGCTGTTCTGCGCGAGTTGCTTTTACTGTTTGTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACTAGGCTTGCGTTCGCGGCGCGGTGCGtct                               <  2:478671/251‑1 (MQ=255)                                                                                                                                                                                                             gacgcggacgcCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGATCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtctt                              <  1:99808/251‑1 (MQ=255)                                                                                                                                                                                                               cgcggacgcCGTCTTTGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGtcttct                            >  1:33761/1‑251 (MQ=255)                                                                                                                                                                                                                    acgcCGTCTTTGGTGATTTCTAACACGGTGGTATCCCTCGCGTTTTGACACGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGACATGTCAGAAACCGGGGTGTGCTGATCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGCCCCGGTGGGTCTTTTTtctc                       >  2:278370/1‑251 (MQ=255)                                                                                                                                                                                                                             ttGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCTGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACcagcag               <  1:153400/250‑1 (MQ=255)                                                                                                                                                                                                                             ttGGTGATTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGGCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTGACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCTGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACcagcag               <  2:45653/250‑1 (MQ=255)                                                                                                                                                                                                                             ttGGTGAGTTCTAATACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCGTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACcagc                 <  2:21933/248‑1 (MQ=255)                                                                                                                                                                                                                                         aaTACGGTGGTATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTc  <  2:188850/251‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 ACCGTCGCATGCTGCGTCTCTTCCTTTAATACCGGAATTTGATCAGCACGCGTGATATCTTCTACAGCGAAGGTCTGGCCTGCTATTGCAAGCAGGCCAGCTAACGCGGTAAGCCTAAAAAACATGTTCATGCCTGGCCCGGCCTCCGTTTCAGAACACCAGGTGTTCTGCGCGCACAATCAAAGACATACCCGAATTCAGCTGGACGCGGACGCCGTCTTTGGTGATTTCTAATACGGTGGCATCCATCGCGTTTTGACCCGCTTTCACCTTCAGGGCTTGTCCGACAGTCAGAGCTGAAATGTCAGAAACCGGGGTGTGCTGTTCTTCGCGAGGTGCTTTTACTGTTTTTGGCGCTTTCTCTACCGGCTTTTGCGCACGAGGTTTACGTTCAGCGCCTTCTTTGCGGCGTGGCGTAGTCGGACGTGGCTTGCGTTCGCGGCGCGGTGCGTCTTCTTTCTCACCAGCAGTTGCGGCAGCTTC  >  NC_000913/1914686‑1915168

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