Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,672,484	A→T	6.2%	L175F (TTA→TTT)	ydgD →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,672,484	0	A	T	6.2%	219.3 / 4.9	64	L175F (TTA→TTT)	ydgD	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (20/40);  new base T (2/2);  total (22/42)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

TTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTT  >  NC_000913/1672341‑1672632                                                                                                                                                |                                                                                                                                                     ttAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGc                                                                                                                                                 <  2:567916/149‑1 (MQ=255)   aaaaGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCgg                                                                                                                                               <  2:186801/149‑1 (MQ=255)       gCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAgg                                                                                                                                           <  2:2567/149‑1 (MQ=255)         aGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTc                                                                                                                                         <  3:301782/149‑1 (MQ=255)                gATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGt                                                                                                                                  <  3:434713/149‑1 (MQ=255)                gATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGt                                                                                                                                  <  1:256181/149‑1 (MQ=255)                  tGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTCCGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGa                                                                                                                                <  1:102867/149‑1 (MQ=255)                     ttGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTc                                                                                                                             >  3:233187/1‑149 (MQ=255)                     ttGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTc                                                                                                                             >  3:322374/1‑149 (MQ=255)                      tGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTca                                                                                                                            <  4:41929/149‑1 (MQ=255)                        gATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTcagg                                                                                                                          <  4:241767/149‑1 (MQ=255)                             tACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCaggcaggc                                                                                                                     <  4:41331/149‑1 (MQ=255)                             tACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCaggcaggc                                                                                                                     >  4:444900/1‑149 (MQ=255)                                   ccGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCt                                                                                                               >  2:118781/1‑149 (MQ=255)                                    cGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTg                                                                                                              >  3:556841/1‑149 (MQ=255)                                     gCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGa                                                                                                             <  3:586553/149‑1 (MQ=255)                                      cAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                                                                                            <  3:444905/149‑1 (MQ=255)                                      cAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                                                                                            <  3:235503/149‑1 (MQ=255)                                      cAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGaa                                                                                                            <  1:525073/149‑1 (MQ=255)                                      cAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCCACCCTGaa                                                                                                            <  3:539551/149‑1 (MQ=255)                                         ccGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAGGTGACTCAGGCAGGCTACCCTGAAGat                                                                                                         <  3:546492/149‑1 (MQ=255)                                         ccGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGat                                                                                                         <  2:326801/149‑1 (MQ=255)                                                 tGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGa                                                                                                 <  1:118781/149‑1 (MQ=255)                                                      cTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGt                                                                                            <  2:537888/149‑1 (MQ=255)                                                      cTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGt                                                                                            <  3:24599/149‑1 (MQ=255)                                                       ttCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGtt                                                                                           <  4:302329/149‑1 (MQ=255)                                                               gattgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAACGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtc                                                                                   <  1:542040/149‑1 (MQ=255)                                                                attgTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGtca                                                                                  >  1:185142/1‑149 (MQ=255)                                                                       tACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAAc                                                                           >  1:226282/1‑149 (MQ=255)                                                                         cGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtg                                                                         >  4:203481/1‑149 (MQ=255)                                                                            aaTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgtg                                                                       >  4:176767/1‑148 (MQ=255)                                                                            aaTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACtgtg                                                                       <  3:176767/148‑1 (MQ=255)                                                                             aTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGaa                                                                     <  1:411693/149‑1 (MQ=255)                                                                             aTCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGaa                                                                     <  3:203481/149‑1 (MQ=255)                                                                              tCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTaaa                                                                                                                                    <  4:163958/85‑1 (MQ=255)                                                                              tCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTaaa                                                                                                                                    >  3:163958/1‑85 (MQ=255)                                                                                           aTTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAg                                                                                                                                                  >  3:331861/1‑58 (MQ=255)                                                                                           aTTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAg                                                                                                                                                  <  4:331861/58‑1 (MQ=255)                                                                                            ttACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGgcg                                                      >  1:150439/1‑149 (MQ=255)                                                                                                gccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCaaa                                                  <  4:302347/149‑1 (MQ=255)                                                                                                gccgttgccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCaaa                                                  <  1:297835/149‑1 (MQ=255)                                                                                                      gccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCgg                                            <  4:322374/149‑1 (MQ=255)                                                                                                       ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGt                                           >  1:297505/1‑149 (MQ=255)                                                                                                       ccgttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCGGAAGAGCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTTGGCGCAAACGGCGGt                                           <  2:413696/149‑1 (MQ=255)                                                                                                          ttATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGAt                                        >  3:147347/1‑149 (MQ=255)                                                                                                            aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGt                                      <  3:444900/149‑1 (MQ=255)                                                                                                            aTTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGt                                      <  2:150439/149‑1 (MQ=255)                                                                                                             tttGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTc                                     >  3:479407/1‑149 (MQ=255)                                                                                                                   ggAGATAAATCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGAGTCAGGCAGGCTACCCTTAAGATCATCACGATACTCTGTACAGTCATCAAAACTGTTAAGTGACTGCCTGGGCGCAAACGTCGGTGATATCACATCa                               >  2:131104/1‑149 (MQ=255)                                                                                                                      gaTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTg                            >  1:256541/1‑149 (MQ=255)                                                                                                                       aTAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGc                           <  2:256541/149‑1 (MQ=255)                                                                                                                          aaGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGAt                        <  4:513329/149‑1 (MQ=255)                                                                                                                           aGCCGCGCTTACCGCCGCATTATAAGCGGCAGGTCGTAGAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATa                       <  1:364149/149‑1 (MQ=255)                                                                                                                             ccGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATAcc                     >  3:227217/1‑149 (MQ=255)                                                                                                                                   gacctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGAc                                                                >  3:512232/6‑100 (MQ=255)                                                                                                                                   gacctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGAc                                                                <  4:512232/95‑1 (MQ=255)                                                                                                                                      ctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCaaa                                                  <  4:23810/109‑1 (MQ=255)                                                                                                                                      ctgccgcTTTTAATGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCaaa                                                  >  3:23810/3‑111 (MQ=255)                                                                                                                                      ccgccgCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGAGATGTCACATCAGTGCGATACCTTGCCGGGt            <  4:233187/149‑1 (MQ=255)                                                                                                                                        gccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGa          >  1:32509/1‑149 (MQ=255)                                                                                                                                        gccgcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGa          <  2:383262/149‑1 (MQ=255)                                                                                                                                           gcATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAg       <  4:296646/149‑1 (MQ=255)                                                                                                                                              ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCgg    >  4:279320/1‑149 (MQ=255)                                                                                                                                              ttAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCgg    <  2:297505/149‑1 (MQ=255)                                                                                                                                                aaaaGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGCGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGtt  <  1:521275/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     TTAAAAGCAGATGGGGATGGTTGGATTGTACCTCCCGCAGCCGCGCCGTGGGACTTCGGATTGATTGTGCTACGTAATCCCCCTTCTGGCATTACGCCGTTGCCGTTATTTGAGGGAGATAAAGCCGCGCTTACCGCCGCATTAAAAGCGGCAGGTCGTAAAGTGACTCAGGCAGGCTACCCTGAAGATCATCTCGATACGTTGTACAGTCATCAAAACTGTGAAGTGACTGGCTGGGCGCAAACGTCGGTGATGTCACATCAGTGCGATACCTTGCCGGGTGACAGCGGTT  >  NC_000913/1672341‑1672632

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