Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,482,581	A→G	5.2%	D419D (GAT→GAC)	valS ←	valine‑‑tRNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,482,581	0	A	G	5.2%	238.3 / 2.5	78	D419D (GAT→GAC)	valS	valine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (33/40);  new base G (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGC  >  NC_000913/4482436‑4482726                                                                                                                                                  |                                                                                                                                                  gACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATcgc                                                                                                                                               <  2:1193169/150‑1 (MQ=255) gACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATcgc                                                                                                                                               <  2:1230634/150‑1 (MQ=255)   cGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCgcgc                                                                                                                                             <  1:1582715/150‑1 (MQ=255)      aGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATc                                                                                                                                          >  1:1479862/1‑150 (MQ=255)            acaTCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGgag                                                                                                                                    >  1:995067/1‑150 (MQ=255)                cAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGt                                                                                                                                <  2:1332196/150‑1 (MQ=255)                 aGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTa                                                                                                                               <  2:273426/150‑1 (MQ=255)                    aCCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACAt                                                                                                                            <  1:590185/150‑1 (MQ=255)                     ccGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATg                                                                                                                           <  2:591037/150‑1 (MQ=255)                           ttATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCg                                                                                                                     <  2:50328/150‑1 (MQ=255)                             aTTTTCTTTACGCACTTCGTCTTCGTTTCGGCCAACATAAACGTTACCCGCTTCGTCATAACATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCATTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTa                                                                                                                   <  2:1102429/150‑1 (MQ=255)                             aTTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTa                                                                                                                   >  1:325543/1‑150 (MQ=255)                                   tttACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCAGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCtt                                                                                                             >  2:661016/1‑150 (MQ=255)                                    ttACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTc                                                                                                            >  2:775587/1‑150 (MQ=255)                                      aCGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCgg                                                                                                          >  2:285712/1‑150 (MQ=255)                                         cACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATGTCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTAc                                                                                                       >  1:650296/1‑150 (MQ=255)                                         cACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATGTCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTAc                                                                                                       >  1:332214/1‑150 (MQ=255)                                           cttcgtcttcgtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGa                                                                                                     >  2:195362/1‑150 (MQ=255)                                           cttcgtcttcgtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGa                                                                                                     >  2:647668/1‑150 (MQ=255)                                             tcgtcttcgtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAAc                                                                                                   >  1:740606/1‑150 (MQ=255)                                              cgtcttcgtTGCGGCCAACGTAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACt                                                                                                  >  2:414997/1‑150 (MQ=255)                                              cgtcttcgtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACt                                                                                                  >  2:1200508/1‑150 (MQ=255)                                                tcttcgtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGa                                                                                                >  1:1301040/1‑150 (MQ=255)                                                    cgtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGt                                                                                            >  2:902690/1‑150 (MQ=255)                                                     gtTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTc                                                                                           >  1:272971/1‑150 (MQ=255)                                                      ttGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCg                                                                                          >  1:236901/1‑150 (MQ=255)                                                       tGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGc                                                                                         <  2:928058/150‑1 (MQ=255)                                                         cGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCg                                                                                       <  2:1487223/150‑1 (MQ=255)                                                            ccAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTtc                                                                                    <  1:26570/150‑1 (MQ=255)                                                            ccAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGACGtta                                                                                    >  1:1079352/1‑149 (MQ=255)                                                              aaCATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTtctc                                                                                  <  2:995067/150‑1 (MQ=255)                                                                cATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTAGGTACGAACTGAATGTCGCCGTTCTCaa                                                                                >  2:716136/1‑150 (MQ=255)                                                                  tAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGCCCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAAcc                                                                              >  2:1595047/1‑150 (MQ=255)                                                                  tAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAAcc                                                                              >  1:759200/1‑150 (MQ=255)                                                                  tAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAAcc                                                                              <  1:79038/150‑1 (MQ=255)                                                                  tAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAAcc                                                                              <  1:528106/150‑1 (MQ=255)                                                                  tAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAAcc                                                                              <  2:1548574/150‑1 (MQ=255)                                                                    aaCGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGc                                                                            <  1:213490/150‑1 (MQ=255)                                                                    aaCGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGc                                                                            >  2:274986/1‑150 (MQ=255)                                                                     aCGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACTACTGACGAGAGATACACGAGTACTGAATTTCGCGCATCCAGGAGAAGTACATGTTTTCGTTCTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCt                                                                           <  2:537592/150‑1 (MQ=255)                                                                          aCCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCaag                                                                      <  1:1532394/150‑2 (MQ=255)                                                                          aCCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAAc                                                                      <  2:820965/150‑1 (MQ=255)                                                                           cccGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAccc                                                                     >  1:681130/1‑150 (MQ=255)                                                                           cccGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAAcc                                                                     <  1:647668/150‑1 (MQ=255)                                                                               cTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAAccgccg                                                                 <  1:285712/150‑1 (MQ=255)                                                                               cTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAAccgccg                                                                 <  1:107330/150‑1 (MQ=255)                                                                                 tCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGt                                                               <  1:538428/150‑1 (MQ=255)                                                                                  cGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATGTCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGtt                                                              <  2:650296/150‑1 (MQ=255)                                                                                  cGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATGTCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGtt                                                              <  2:332214/150‑1 (MQ=255)                                                                                   gTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGttt                                                             >  1:356216/1‑150 (MQ=255)                                                                                       tACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGcc                                                         <  1:444321/150‑1 (MQ=255)                                                                                         ccATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAg                                                       >  1:1554297/1‑150 (MQ=255)                                                                                                gggATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCg                                                >  1:1093863/1‑150 (MQ=255)                                                                                                 ttaTACGGTGACCCGACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTAGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCgg                                               <  2:1315307/148‑1 (MQ=255)                                                                                                 ggATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGTAGAAGTAAATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACAGCTTCAACTGTCTGTTTCTCGAGGACATAgg                                               <  2:1286799/150‑1 (MQ=255)                                                                                                 ggATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCgg                                               <  2:598130/150‑1 (MQ=255)                                                                                                    tACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcac                                            <  2:1479862/150‑1 (MQ=255)                                                                                                    tACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcac                                            <  2:221250/150‑1 (MQ=255)                                                                                                      cGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcacgc                                          <  1:232399/150‑1 (MQ=255)                                                                                                      cGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcacgc                                          >  1:534518/1‑150 (MQ=255)                                                                                                      cGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcacgc                                          >  1:1598484/1‑150 (MQ=255)                                                                                                            ccccACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTAc                                    >  1:397444/1‑150 (MQ=255)                                                                                                            ccccACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTAc                                    >  2:934797/1‑150 (MQ=255)                                                                                                                accacAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACt                                <  2:236901/150‑1 (MQ=255)                                                                                                                 ccacAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcacgca                                         >  1:509941/1‑140 (MQ=255)                                                                                                                 ccacAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGgcacgca                                         <  2:509941/140‑1 (MQ=255)                                                                                                                      aCTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTggtcgg                          <  2:1301040/150‑1 (MQ=255)                                                                                                                        tGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTggtcggtc                        <  1:195362/150‑1 (MQ=255)                                                                                                                               agaTACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTg                 <  2:1161086/150‑1 (MQ=255)                                                                                                                                 aTACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGt               >  2:305320/1‑150 (MQ=255)                                                                                                                                   acacCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTc             <  2:468051/150‑1 (MQ=255)                                                                                                                                    cacCAGTCCTGAATATCCCGCATACAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAAACGCCGGTTTCGCCAGGACATCGGCACGCACCTACCACTGGTCGCTCAGCATTGGTTCg            >  2:689683/1‑150 (MQ=255)                                                                                                                                       cAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTAGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATa         <  1:716136/150‑1 (MQ=255)                                                                                                                                            cTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTAc    <  2:740606/150‑1 (MQ=255)                                                                                                                                             tGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTAcg   <  1:661016/150‑1 (MQ=255)                                                                                                                                              gAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTc             <  1:1368885/139‑1 (MQ=255)                                                                                                                                              gAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTc             >  2:1368885/1‑139 (MQ=255)                                                                                                                                              gAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTAcgc  >  1:176401/1‑150 (MQ=255)                                                                                                                                                  |                                                                                                                                                  GACGCAGGACAACATCAGCACCGAGGTTATTTTCTTTACGCACTTCGTCTTCGTTGCGGCCAACATAAACGTTACCCGCTTCGTCATACCATGCCGGGATACGGTGACCCCACCACAACTGACGAGAGATACACCAGTCCTGAATATCGCGCATCCAGGAGAAGTACATGTTTTCGTACTGCTTCGGTACGAACTGAATGTCGCCGTTCTCAACCGCTTCAACCGCCGGTTTCGCCAGGACATCGGCACGCACGTACCACTGGTCGGTCAGCATTGGTTCGATAACTACGC  >  NC_000913/4482436‑4482726

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