Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	718,659	0	T	A	42.4%	102.3 / 34.6	66	D601V (GAT→GTT)	speF	ornithine decarboxylase isozyme, inducible
Reads supporting (aligned to +/- strand):  ref base T (6/32);  new base A (28/0);  total (34/32)
Fisher's exact test for biased strand distribution p-value = 4.05e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGA  >  NC_000913/718511‑718805                                                                                                                                                      |                                                                                                                                                   gCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCa                                                                                                                                                  <  2:975966/151‑1 (MQ=255)  ccTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACa                                                                                                                                                  >  1:480556/1‑150 (MQ=255)        gCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCa                                                                                                                                              <  1:981920/148‑1 (MQ=255)        gCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCAttt                                                                                                                                           <  1:525847/151‑1 (MQ=255)            ccGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAC‑AAAACAAACATTTTCt                                                                                                                                        >  2:337168/1‑150 (MQ=255)                 aGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAAAATGCAttt                                                                                                                                           >  1:285324/1‑142 (MQ=255)                  gACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACaa                                                                                                                                 <  1:1091362/151‑1 (MQ=255)                        ccAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATTCATTTTCTGACACAAATGa                                                                                                                             >  1:644250/1‑149 (MQ=255)                          aGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCa                                                                                                                         <  1:559956/151‑1 (MQ=255)                             tCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAgg                                                                                                                       <  2:781388/150‑1 (MQ=255)                                aCTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGtat                                                                                                                   <  1:938336/151‑1 (MQ=255)                                 cTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCGTAC‑AAAAAATTCATTTTCTGTCACAACTTACGCCggg                                                                                                                      >  1:420270/1‑147 (MQ=255)                                   tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGtata                                                                                                                  <  2:193093/149‑1 (MQ=255)                                   tCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATaa                                                                                                                 <  1:665519/150‑1 (MQ=255)                                      ccGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAAAATACATTTTCTCACCCAACCGACGCAgggggtaa                                                                                                                 >  1:811710/1‑142 (MQ=255)                                        gcgTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCg                                                                                                           <  1:567052/151‑1 (MQ=255)                                            aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAAAAAGCATTTTCTTACACCAATGACGCAGGGGATaa                                                                                                                 >  2:960963/1‑141 (MQ=255)                                            aaaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCAAAA‑AAAACAAGCATTTTCTTACCACACTGAAGCCGGGGAAAAACGGCggcg                                                                                                        >  2:694392/1‑147 (MQ=255)                                             aaTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAgcg                                                                                                      <  1:961621/151‑1 (MQ=255)                                              aTAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATGCATTTTCTGACACAAATGACGCAGGGGATAAACCGCGGAgc                                                                                                       >  1:370519/1‑149 (MQ=255)                                                  gCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCtt                                                                                                 <  2:1049088/151‑1 (MQ=255)                                                    gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCa                                                                                               <  1:805103/151‑1 (MQ=255)                                                    gTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAAAATGCATTTTCTTACAAAACTGACGCCGGGGGTAAACCGCGGAGCGCTCtt                                                                                                 >  1:113833/1‑149 (MQ=255)                                                     tAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGGTAC‑AAAACATGCCTTTCTTTTCACCACTGACGCAGGGGTTAAACGCCGGAGCGCGCtt                                                                                                 >  2:44514/1‑148 (MQ=255)                                                       gTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGc                                                                                            <  1:1035471/151‑1 (MQ=255)                                                        ttGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATGCATTTTCTGTCACAACTGACGCAGGGGATAAACGGCGGAGCGCTCtt                                                                                                 >  1:677704/1‑145 (MQ=255)                                                           gCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGtt                                                                                        <  2:757245/151‑1 (MQ=255)                                                                 tGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCCTAC‑AAAACATGCATTTTCTGACACAAATGACGCAGGGGATAAACGGCGTAGGGCCCTTCATGCTGttt                                                                                       >  1:70301/1‑146 (MQ=255)                                                                   cGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAA‑AAAACATGCATTTTCTGGCACAACTGACGCAGGGGATAAACCGCGGAGCGCTCCTCCTGCTGttttttg                                                                                   >  1:395609/1‑144 (MQ=255)                                                                    ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATgg                                                                               <  1:190957/151‑1 (MQ=255)                                                                    ggATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAA‑AAAATATTCCTTTTCTGACACAAAATACGCCGGGGATAAACGGCGTAACGCTCtttctgtttttttttg                                                                                   >  1:947715/1‑133 (MQ=255)                                                                           tGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAA‑AAAAAAAGAATTTTCTGACCCAACCGAAGCCGGGGATAAACCGCGCAACGCCTTTCATGCTGTTTTTAGAAGGAAGGt                                                                          >  1:254682/1‑149 (MQ=255)                                                                            gCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAg                                                                        <  1:356891/150‑1 (MQ=255)                                                                             cTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATACAAAAACTTCCATTT‑CTTTTACAACTTTCGCAGGGGAAAACCGGCGCAGTGCTCTTCTTTCTCTTTTTTGATGGGAg                                                                            >  1:623209/1‑145 (MQ=255)                                                                                   cgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATGCATTTTCTTACACAAATGACGCCGGGGATAAACCGCGGAGCGCTCTTTATTCTGTTTGTTGATGGAAGGGAGGACtt                                                                   >  1:136204/1‑148 (MQ=255)                                                                                   cgTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATGCATTTCCTGACACAACTGACGCCGGGGGTAAACCGCGTGGCGCTCTTCCTGCTGTTTTTTGATGGGAGGGAGCACtt                                                                   >  1:671502/1‑148 (MQ=255)                                                                                    gTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGAAGCACTTCTGc                                                               <  2:686316/151‑1 (MQ=255)                                                                                         aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                          <  1:190588/151‑1 (MQ=255)                                                                                         aaGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATg                                                          <  2:64407/151‑1 (MQ=255)                                                                                             gCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATACTTTTTCTGCCACAACTATCGCAGGGGGTAACCGGCGGAGCGCCCTTCTTGTTGTTTGTAGATGGGAGGGAGCAGTTCTTCTTAt                                                           >  2:738813/1‑146 (MQ=255)                                                                                              cTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAAGTGACGCAGGGTATAACCGGCGTAGCGCTCTGCATGCTGTTTGTAGATGGAAGGTAGCGCTTCTGCTAATGgcgt                                                      <  2:287582/150‑2 (MQ=255)                                                                                                ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGt                                                    <  2:129396/150‑1 (MQ=255)                                                                                                ccTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTc                                                   <  1:365553/151‑1 (MQ=255)                                                                                                        aaCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATTTATTTTCTTACACAACAGACGCAGGGGGTAACCGGCGGAGCGCTCTTCATGCTGTTTTTTGATGGGAGGGAGCACTTCTTCTAATGGCgcggcg                                                  >  2:748297/1‑144 (MQ=255)                                                                                                         aCATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAg                                          <  1:1045974/151‑1 (MQ=255)                                                                                                           aTCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAAAATTCCTTTTCTTACACAAATGACGCAGGGGATAAACCGCGGAGCGCCCTTTCTGCTGTTTGTTGATGGGAGGGAGCACTTCTTCTAATTGCGCGTCGGGCTCaa                                           >  1:315982/1‑148 (MQ=255)                                                                                                            tctcTTTTTGCAGTTGTTTCACGTTGTGGCGGGGATAC‑AAAAAATGCATTTTCTTACAACACTGAAGCAGGGGATAAACGGCGGAGGGCTCTTCATGCTGTTTGTTGATGGAAGGGAGGACTTCTTCTAATTGGGCGGCGGGCTCaa                                           >  1:570924/1‑147 (MQ=255)                                                                                                             ctctTTTTGCAGTTGTTTCACGTTGTGGCGGGCCAAA‑AAAACATTTATTTTCTTACACAACTGACGCAGGGGGTAAACGGCGGAGCGCTCTTTCTTGTTTTTTTTGATGGAAGGGAGAAATTTTGCTTATGGGGCGTCGGGCTCAAACa                                        >  2:1069502/1‑149 (MQ=255)                                                                                                              tctTTTTGCAGTTGTTTCACGTTGTGGCGGGCTTACAAAATATTGTTTTTCTTTTCAAAAC‑GACGCCGGGGGTAACAGGCGGCGGGCTCCTTATTCTGTTTGTTGATGGGAGGTAGGACTTCTTCTTATGTGGGGGCGGGCCCAAGCAggt                                     >  2:22142/1‑149 (MQ=255)                                                                                                               ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAg                                       <  1:174249/148‑1 (MQ=255)                                                                                                               ctTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAg                                       <  1:1027595/148‑1 (MQ=255)                                                                                                                  tttGCAGTTGTTTCACGTTGTGGCGGGCCTAC‑AAAAAATGTATTTTCTTGTACCACTGACGCCGGGGGGAAACGGCGGAGCGCTCTTCTTTCTGTTTTTTGATGGGAGGGAGCACCTCTGCTCATTGCGCGTCGGGCCCAAACAGttttt                                  >  2:487569/1‑150 (MQ=255)                                                                                                                   ttGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTc                                 <  2:404515/150‑1 (MQ=255)                                                                                                                   ttGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTc                                 <  1:190854/150‑1 (MQ=255)                                                                                                                     gCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAAACATGCATTTTCTTACACAACAGACGCAGGGGATAAACCGCGTAGCGCTCTTCATGCTGTTTTTTGATGGAAGGGAGCACTTCTGCTAATGGCGCGTCGGGCTCCAGCAGttttttta                               >  1:832521/1‑147 (MQ=255)                                                                                                                            tttCACGTTGTGGCGGGCATAC‑AAAACATGCATTTCCTGACACAAATGACGCAGGGGATAAACCGCGTAGCGCTCTTTATGCTGTTTGTTGATGGAAGGGAGCACTTCTTCTAATGGCGCGGCGGGCTCAAGCAGTTTTTTGAAGCGGAcc                       >  1:242098/1‑151 (MQ=255)                                                                                                                                aCGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcagca                   <  2:37470/151‑1 (MQ=255)                                                                                                                                aCGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcagc                    <  2:971923/150‑1 (MQ=255)                                                                                                                                 cGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACcagcag                  <  2:707649/151‑1 (MQ=255)                                                                                                                                   ttGTGGCGGGCATAC‑AAAACATGCATTTCCTGACACAACTGACGCAGGGGATAACCCGCGTAGCGCTCTTCATGCTGTTTGTAGATGGGAGGGAGCACTTCTGCTAATGGCGCGTCGGGCTCAAGCAGTTTTTTGAAACGTACCAGCAggg                >  1:924126/1‑151 (MQ=255)                                                                                                                                       ggCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGcaac            <  1:627295/151‑1 (MQ=255)                                                                                                                                       ggCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTCCGAAGCGTACCAGCAGGGcaa             >  2:740736/1‑150 (MQ=255)                                                                                                                                              aTAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGc     >  1:263754/1‑151 (MQ=255)                                                                                                                                                 c‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTg   >  2:91628/1‑150 (MQ=255)                                                                                                                                                   aaaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa  >  1:407049/1‑150 (MQ=255)                                                                                                                                                    aaTCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGa  >  2:787279/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCCTTCTGCATCCGGCAGACGAACCAGTTCCACTTCACCGCGTAAATAGGCGTAGTTGGCTTCTTGCGGATTCATGCTGACGCGTGGGAAGTGCTCCTTACGGAACATCTCTTTTTGCAGTTGTTTCACGTTGTGGCGGGCATAC‑AAATCATGCATTTCCTGACACAACTGACGCAGGGTATAACCGGCGTAGCGCTCTTCATGCTGTTTGTAGATGGAAGGTAGCACTTCTGCTAATGGCGCGTCGGACTCAAGCAGTTTTTCGAAGCGTACCAGCAGGGCAACAAGTTGCTGA  >  NC_000913/718511‑718805

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