Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	817,739	T→C	5.2%	G232G (GGT→GGC)	moaA →	GTP 3',8'‑cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	817,739	0	T	C	5.2%	233.4 / 2.8	76	G232G (GGT→GGC)	moaA	GTP 3',8'‑cyclase
Reads supporting (aligned to +/- strand):  ref base T (30/42);  new base C (2/2);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GCCCTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTT  >  NC_000913/817594‑817886                                                                                                                                                  |                                                                                                                                                    gCCCTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCcg                                                                                                                                                 >  1:1016874/1‑150 (MQ=255)  cccTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCcgc                                                                                                                                                <  2:1237728/150‑1 (MQ=255)     tATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCa                                                                                                                                             >  1:286711/1‑150 (MQ=255)       tCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAg                                                                                                                                           <  2:69236/150‑1 (MQ=255)           gCTGCGATTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTATTCCGTAAGCATCACATCTCTGGTCAGGTTATGCGTGACCAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCtt                                                                                                                                       >  2:1325666/1‑150 (MQ=255)             tGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCtttt                                                                                                                                     <  2:1375646/150‑1 (MQ=255)              gCGTTTCATCGAACTGTTGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTg                                                                                                                                    <  1:126347/150‑1 (MQ=255)               cGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGc                                                                                                                                   <  1:792574/150‑1 (MQ=255)               cGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGc                                                                                                                                   <  2:1047828/150‑1 (MQ=255)                gTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGcc                                                                                                                                  >  1:616597/1‑150 (MQ=255)                    cATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATcc                                                                                                                              >  1:669027/1‑150 (MQ=255)                      tCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAg                                                                                                                            <  1:109376/150‑1 (MQ=255)                       cGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGa                                                                                                                           >  2:657478/1‑150 (MQ=255)                       cGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGa                                                                                                                           >  2:1237048/1‑150 (MQ=255)                         aaCTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGGCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGAtt                                                                                                                         <  2:924282/150‑1 (MQ=255)                           cTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGCCCCGCGCAAGTCTTTTGCCATCCAGATTAc                                                                                                                       >  1:1106728/1‑150 (MQ=255)                           cTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGCCCCGCGCAAGTCTTTTGCCATCCAGATTAc                                                                                                                       >  1:1101550/1‑150 (MQ=255)                               tGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCg                                                                                                                   >  1:1166796/1‑150 (MQ=255)                                ggAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCgg                                                                                                                  <  1:279728/150‑1 (MQ=255)                                        gCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCGCGCGCAAGTCTTTTGCCAT‑CAGATTACGCCGGAGAGTTTgg                                                                                                         <  2:53907/150‑1 (MQ=255)                                        gCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTg                                                                                                          <  2:389800/150‑1 (MQ=255)                                        gCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTg                                                                                                          >  1:141378/1‑150 (MQ=255)                                               cAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTAt                                                                                                   >  1:925957/1‑150 (MQ=255)                                                       tctTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGt                                                                                           >  2:1213316/1‑150 (MQ=255)                                                          tCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATg                                                                                        <  2:72948/150‑1 (MQ=255)                                                          tCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATg                                                                                        <  2:656055/150‑1 (MQ=255)                                                           ccGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGa                                                                                       <  2:967307/150‑1 (MQ=255)                                                             gTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGaaa                                                                                     >  1:894307/1‑150 (MQ=255)                                                                aGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAg                                                                                  <  1:30940/150‑1 (MQ=255)                                                                aGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAg                                                                                  <  1:517108/150‑1 (MQ=255)                                                                aGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAg                                                                                  <  2:1120378/150‑1 (MQ=255)                                                                aGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAg                                                                                  <  1:934937/150‑1 (MQ=255)                                                                     cacaTCTCTGGTCAGGTTCTGCGTGACGAGCTGCTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTc                                                                             <  2:699105/150‑1 (MQ=255)                                                                     cacaTCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTc                                                                             <  1:1237048/150‑1 (MQ=255)                                                                         tctcTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTgcg                                                                         <  1:123656/150‑1 (MQ=255)                                                                            ctGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCa                                                                      >  1:1014547/1‑150 (MQ=255)                                                                                tCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTg                                                                  >  1:966639/1‑150 (MQ=255)                                                                                 cAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGc                                                                 <  2:432651/150‑1 (MQ=255)                                                                                  aGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCa                                                                <  1:838562/150‑1 (MQ=255)                                                                                   ggTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCTCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCaa                                                               <  2:683256/150‑1 (MQ=255)                                                                                   ggTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCaa                                                               >  1:1393722/1‑150 (MQ=255)                                                                                   ggTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCaa                                                               >  2:281708/1‑150 (MQ=255)                                                                                   ggTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCaa                                                               >  1:450583/1‑150 (MQ=255)                                                                                   ggTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCaa                                                               >  2:6419/1‑150 (MQ=255)                                                                                   ggTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCaa                                                               <  2:1016874/150‑1 (MQ=255)                                                                                      tCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCg                                                            <  1:465058/150‑1 (MQ=255)                                                                                        tGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGcc                                                          <  1:85220/150‑1 (MQ=255)                                                                                                gAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGtt                                                  >  1:430587/1‑150 (MQ=255)                                                                                                  gCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTtc                                                >  1:1346459/1‑150 (MQ=255)                                                                                                    tACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTtcct                                              >  2:84750/1‑150 (MQ=255)                                                                                                      cTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTtcctcg                                            <  2:1393722/150‑2 (MQ=255)                                                                                                      cTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTtcctcc                                            >  2:673518/1‑150 (MQ=255)                                                                                                         cGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCAtt                                         <  1:822147/150‑1 (MQ=255)                                                                                                          gTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTg                                        <  1:657478/150‑1 (MQ=255)                                                                                                           tCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTgg                                       <  2:616597/150‑1 (MQ=255)                                                                                                                gCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAAc                                  >  2:737519/1‑150 (MQ=255)                                                                                                                gCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAAc                                  <  2:894307/150‑1 (MQ=255)                                                                                                                   ggATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGcc                                                          <  2:2629/123‑1 (MQ=255)                                                                                                                   ggATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGcc                                                          >  1:2629/1‑123 (MQ=255)                                                                                                                    gATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCa                              <  1:1348003/150‑1 (MQ=255)                                                                                                                    gATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCa                              >  2:612480/1‑150 (MQ=255)                                                                                                                        caccaATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCAtctc                          <  1:731092/150‑1 (MQ=255)                                                                                                                        caccaATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCAtctc                          >  1:1400802/1‑150 (MQ=255)                                                                                                                          ccaATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTg                        <  2:1346459/150‑1 (MQ=255)                                                                                                                          ccaATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTg                        >  1:805402/1‑150 (MQ=255)                                                                                                                            aaTTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGcc                      >  1:60475/1‑150 (MQ=255)                                                                                                                              ttACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTg                    <  1:894515/150‑1 (MQ=255)                                                                                                                                 cGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCTTATGAAAAAGACTTCTGCGCCACTTGCAACCTCCTGCGCGTTTCCTTCATTGGTAAACTCCATCTATGCCTGttt                 <  2:669027/150‑1 (MQ=255)                                                                                                                                   tCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTgg               <  1:673518/150‑1 (MQ=255)                                                                                                                                   tCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTgg               <  2:1189357/150‑1 (MQ=255)                                                                                                                                    cAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGCCTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGt              >  2:935691/1‑150 (MQ=255)                                                                                                                                      aCGCAGCGACGGCCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGa            <  2:1106728/150‑1 (MQ=255)                                                                                                                                      aCGCAGCGACGGCCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGa            <  2:1101550/150‑1 (MQ=255)                                                                                                                                       cGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTACAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGaa           >  2:881185/1‑150 (MQ=255)                                                                                                                                                ggTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGtt  <  2:141378/150‑1 (MQ=255)                                                                                                                                                ggTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGtt  <  2:1398343/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    GCCCTATCCAGCTGCGTTTCATCGAACTGATGGAAACGGGCGAGGGCAGCGAGCTCTTCCGTAAGCATCACATCTCTGGTCAGGTTCTGCGTGACGAGCTACTGCGTCGCGGCTGGATCCACCAATTACGTCAACGCAGCGACGGTCCCGCGCAAGTCTTTTGCCATCCAGATTACGCCGGAGAGATTGGCCTTATCATGCCGTATGAAAAAGACTTCTGCGCCACTTGCAACCGCCTGCGCGTTTCCTCCATTGGTAAACTCCATCTCTGCCTGTTTGGTGAAGGCGGCGTT  >  NC_000913/817594‑817886

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