Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	192,579	A→T	5.9%	*242L (TAA→TTA)	pyrH →	uridylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	192,579	0	A	T	5.9%	238.5 / 2.7	86	*242L (TAA→TTA)	pyrH	uridylate kinase
Reads supporting (aligned to +/- strand):  ref base A (46/33);  new base T (2/4);  total (48/38)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGA  >  NC_000913/192434‑192724                                                                                                                                                   |                                                                                                                                                  gAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAAtt                                                                                                                                                >  2:502378/1‑149 (MQ=255)        agCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGa                                                                                                                                         <  2:540242/149‑1 (MQ=255)         gCTGTAAGTCATGGACCTGTCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGAt                                                                                                                                        <  1:200447/149‑1 (MQ=255)          cTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATg                                                                                                                                       >  2:988423/1‑149 (MQ=255)          cTGAAAGTAATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGATCAAAGCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GTAATAATTCCCGTGATg                                                                                                                                       <  2:854238/149‑1 (MQ=255)             aaaGGCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGCAGGGACTTTAATCACGGGAATTATT‑CCGTGAttaaa                                                                                                                                    <  1:971141/149‑5 (MQ=255)              aaGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGTTTCGGGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCACGGGAATTATT‑CCGTGAttaaag                                                                                                                                   <  2:47933/149‑6 (MQ=255)               aGTCATGGACCTGGCGGCCTTCACGCTGGCTAGTGCCCATAAATTAGGGAGTCGGGTGTTCAATATGTACAAAGCGGGTGCGCGGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACGTTAATCTC‑GGAATAAATCCCGTGCTGGataa                                                                                                                                  >  2:663835/1‑149 (MQ=255)                   aTGGACCTGGCGGCCTTCACGCTGGCGCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                              <  2:1154658/149‑1 (MQ=255)                   aTGGACCTGGCGGCCCTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGataaataa                                                                                                                              <  1:141942/149‑1 (MQ=255)                       aCCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGATCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGt                                                                                                                          >  2:1043613/1‑149 (MQ=255)                              ggcCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGAtta                                                                                                                   >  2:126509/1‑148 (MQ=255)                               gcCATCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTGTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTcc                                                                                                                  <  2:735611/149‑1 (MQ=255)                                 cTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAAGGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgg                                                                                                                >  2:408529/1‑148 (MQ=255)                                  ttCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCcgcg                                                                                                               <  1:1134641/149‑1 (MQ=255)                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                            <  1:339547/149‑1 (MQ=255)                                     aCGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTaa                                                                                                            >  2:76950/1‑149 (MQ=255)                                       gCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGt                                                                                                          >  1:196807/1‑149 (MQ=255)                                         tGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATgg                                                                                                                                      >  1:839678/1‑119 (MQ=255)                                         tGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATgg                                                                                                                                      <  2:839678/119‑1 (MQ=255)                                           gCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAT‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcg                                                                                                      <  2:604687/149‑1 (MQ=255)                                            cTCGTTAACATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTATTGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgc                                                                                                     <  1:724452/149‑1 (MQ=255)                                            cTCGTGAACATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgc                                                                                                     <  2:381730/149‑1 (MQ=255)                                             tagTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATcgcg                                                                                                    <  1:156798/147‑1 (MQ=255)                                              ggTGACCATAAATTACCGATTCGTGTTTTCAATTTTAATACACCGTGTGCGCTGAGCCGTGTGGTTATGGGTGAAAAAGAAGTGACTTTATTCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgg                                                                                                   <  2:1142275/148‑1 (MQ=255)                                              cGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGCGGTAATGGGGGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCgg                                                                                                   <  1:691686/149‑1 (MQ=255)                                                 gACCCTAAATTACCGAATCGTGTTTGCAATATGAAGAAACCGGGTGCGATGCGCCGTGGGGTAATGGGTGAAAAAGAAGTGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAACATTCCGCGTAAGTATCGCggggg                                                                                                <  2:800715/149‑1 (MQ=255)                                                   ccATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCg                                                                                              >  2:212056/1‑149 (MQ=255)                                                   ccATAAATTACCGATTCGTGTTTTCAAAATGAACAAACCGGGTGCGCTCCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GTAATAATTCCCGTGATGGATAAATAAGGGTAATATTCCGAGTAAGTATCGCGGGGGCg                                                                                              >  2:664225/1‑149 (MQ=255)                                                    cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGt                                                                                             >  1:615960/1‑149 (MQ=255)                                                    cATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTACCGTGATGGATTAATAAGGGTAAGATTCCGAGTAAGTATCGCGGGGCCGt                                                                                             >  1:127456/1‑149 (MQ=255)                                                      tAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTaa                                                                                           >  1:28386/1‑149 (MQ=255)                                                         aTTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGtg                                                                                        >  2:52301/1‑148 (MQ=255)                                                         aTTACCGATTCGTGTTTTAAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGCGCGTAATTc                                                                                        >  2:1069471/1‑149 (MQ=255)                                                          ttACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCt                                                                                       >  2:1029321/1‑149 (MQ=255)                                                              cGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCTTAAGTCTGGtt                                                                                   >  2:200265/1‑149 (MQ=255)                                                              cGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGtt                                                                                   >  1:936467/1‑149 (MQ=255)                                                                aTTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTATGTCTGGTtat                                                                                 <  1:114976/149‑1 (MQ=255)                                                                  tCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                               <  2:196807/149‑1 (MQ=255)                                                                  tCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATaa                                                                               <  2:656347/149‑1 (MQ=255)                                                                     gggtttCAATATGAACAAACCGGGTGCTCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTAAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGCTTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGc                                                                            <  2:791641/146‑1 (MQ=255)                                                                      gtTTTCAATATGAACAAACCGGGGGCGCTGCGCCGGGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGTTGGATATAAAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCg                                                                           <  2:312651/149‑1 (MQ=255)                                                                       ttttCAATATGAACAAACCGGGGGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGt                                                                          <  2:584848/149‑1 (MQ=255)                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTa                                                                        <  1:494556/149‑1 (MQ=255)                                                                         ttCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTa                                                                        >  1:43410/1‑149 (MQ=255)                                                                           cAATATGAACAACCCGGGTGCGCGGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGGTAtt                                                                      <  1:514669/149‑1 (MQ=255)                                                                              tatGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTAttgtt                                                                   >  2:740688/1‑149 (MQ=255)                                                                               atGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGT‑AGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                                 >  2:88008/1‑149 (MQ=255)                                                                                tGAACACACCGGGTGCGCTGCGCCGTGTGGTAATGGGTTAAATAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgc                                                                 <  1:711011/149‑1 (MQ=255)                                                                                 gAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgca                                                                >  1:675246/1‑149 (MQ=255)                                                                                 gAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgca                                                                >  1:553816/1‑149 (MQ=255)                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                               >  2:1061826/1‑149 (MQ=255)                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                               >  1:197477/1‑149 (MQ=255)                                                                                  aaCAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcag                                                               <  1:637658/149‑1 (MQ=255)                                                                                   aCAAACCGGGTGCGCTGCGCCGTGTGGGAATGGGTGAAAAAGAAGGGACTTTAATCAC‑CGAATAATTCACGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcagg                                                              >  2:604614/1‑149 (MQ=255)                                                                                      aaCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTgcaggcag                                                           >  1:847023/1‑149 (MQ=255)                                                                                       aCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGt                                                          >  1:961539/1‑149 (MQ=255)                                                                                       aCCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGt                                                          >  1:753100/1‑149 (MQ=255)                                                                                        ccGCGTGCGCGGAGCCGGGTGGTAATGGGTGAAAAATAATCGACTTTAATCAC‑GGAATTTATCCCGTGATGGATAAATAATGGTAATATACAGCGTAAGGAGCGCGGGGGCGGAAGTCTGGTTTTTAGGCGTTATTGTTGCAGGCAGtt                                                         <  2:765138/149‑1 (MQ=255)                                                                                          ggggGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAAATCCCGGGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                       <  2:215045/149‑1 (MQ=255)                                                                                          gggTGCGCTGCGCCGTGTGGTAATGGGTGAAACAGAAGGGACTTTTATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                       <  1:525189/149‑1 (MQ=255)                                                                                          gggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCGGGCAGTTTg                                                       <  1:988423/149‑1 (MQ=255)                                                                                          gggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATAAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTg                                                       >  2:165287/1‑149 (MQ=255)                                                                                           ggTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGCGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTgg                                                      >  2:88416/1‑149 (MQ=255)                                                                                               gggtggtccGTGTGGTTATGGGTGAAAATTTAGGGACTTTAATCAC‑GGTTTTATTGCAGTGATTGATAAATAAGGGTAAGATTGCGCGTAAGTATCGGGGGGGCGTAAGTCTGGTTATACGGCGTTATTGTTGCAGTCAGTTTGGTGAc                                                  <  1:1049608/142‑1 (MQ=255)                                                                                                 ctgcccCGTGTGGGAATGGGCGAAAACTAAGCGCCTTTAATCAC‑GGAATAATTCCCGTGAGGGATAAATAAGGGTAAGATTCCTCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACgg                                                <  2:1145933/149‑1 (MQ=255)                                                                                                              tAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCgcagca                                   >  1:358823/1‑149 (MQ=255)                                                                                                               aaTGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCaa                                  >  1:358986/1‑149 (MQ=255)                                                                                                                aTGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCGGTTTGGTCACGGCCCGCGCGCAGCAAc                                 >  2:710316/1‑149 (MQ=255)                                                                                                                    gTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGa                             >  2:900706/1‑149 (MQ=255)                                                                                                                      gAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGc                           <  1:1029321/149‑1 (MQ=255)                                                                                                                      gAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGc                           >  2:855219/1‑149 (MQ=255)                                                                                                                       aaaaaGAAGGGACTTTAATCAC‑GGAAAAATTACCGTCAGGTAGAAATAAGGGTAAGATTCAGCGTAAGTATCGAGGCGGAGTAAGTCTGGTTATAAGGAGTTATTGTTGCAGGCAGTTTGGTAACGGCCAGCTCGCAGCAACAGGACCg                          >  1:955332/1‑149 (MQ=255)                                                                                                                       aaaaaGAAGGGACTATAATCAC‑GGAATAATTCCCGTGATGGATAAATAACGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCg                          >  2:517297/1‑149 (MQ=255)                                                                                                                           agaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACa                      <  2:936467/149‑1 (MQ=255)                                                                                                                            gaagGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACaa                     >  2:300591/1‑149 (MQ=255)                                                                                                                               gggACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCACTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAg                  >  2:776974/1‑149 (MQ=255)                                                                                                                                   cTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACg              >  1:906144/1‑149 (MQ=255)                                                                                                                                     tttaTCAC‑GGAATCATGCCCGTGATGTATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGACAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTg            <  1:677222/146‑1 (MQ=255)                                                                                                                                      tAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCACCCGGAGCGTACAAAAGTACGTGa           >  2:876971/1‑149 (MQ=255)                                                                                                                                      ccaTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGa           >  2:1158236/3‑149 (MQ=255)                                                                                                                                       aaTCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAg          <  1:502378/149‑1 (MQ=255)                                                                                                                                        aTCACGGGAATTATT‑CCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg         >  1:754097/1‑149 (MQ=255)                                                                                                                                        aTCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAgg         >  1:24157/1‑149 (MQ=255)                                                                                                                                         tCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGa        >  1:587564/1‑149 (MQ=255)                                                                                                                                           aC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATg      >  1:117161/1‑149 (MQ=255)                                                                                                                                                aataatTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGa  >  1:256841/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GAAAAAGAGCTGAAAGTCATGGACCTGGCGGCCTTCACGCTGGCTCGTGACCATAAATTACCGATTCGTGTTTTCAATATGAACAAACCGGGTGCGCTGCGCCGTGTGGTAATGGGTGAAAAAGAAGGGACTTTAATCAC‑GGAATAATTCCCGTGATGGATAAATAAGGGTAAGATTCCGCGTAAGTATCGCGGGGGCGTAAGTCTGGTTATAAGGCGTTATTGTTGCAGGCAGTTTGGTCACGGCCAGCGCGCAGCAACCGGAGCGTACAAAAGTACGTGAGGATGGCGA  >  NC_000913/192434‑192724

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