Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,604,350	0	C	T	38.6%	81.7 / 36.0	71	S64F (TCT→TTT)	bglJ	bgl operon transcriptional activator
Reads supporting (aligned to +/- strand):  ref base C (8/35);  new base T (27/0);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 2.59e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAT  >  NC_000913/4604203‑4604499                                                                                                                                                     |                                                                                                                                                        aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGc                                                                                                                                                      <  2:289909/151‑1 (MQ=255)   aaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGACGCCGATCTCGTTTGCGGCGG‑TCATCTTTTGATTTTTTGcc                                                                                                                                                     >  1:567605/1‑150 (MQ=255)       gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAGCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATga                                                                                                                                                 <  1:1141611/150‑1 (MQ=255)       gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATg                                                                                                                                                  >  2:236016/1‑149 (MQ=255)       gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATgag                                                                                                                                                <  1:719846/151‑1 (MQ=255)       gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATga                                                                                                                                                 <  1:911033/150‑1 (MQ=255)       gCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCAGTTTTTCTTTTTTTGACAt                                                                                                                                                   >  1:207209/1‑148 (MQ=255)         gTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGaa                                                                                                                                              <  2:453511/151‑1 (MQ=255)           cATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGt                                                                                                                                            <  2:667487/151‑1 (MQ=255)                  agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAgcgc                                                                                                                                      <  1:822962/150‑1 (MQ=255)                    taTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑GCATTTTTTTTTTTTTTGCCATTGGAAGGGGgcgcc                                                                                                                                     >  1:19970/1‑148 (MQ=255)                       cGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTGAGAAGTGAgcgccg                                                                                                                                    >  1:539893/1‑145 (MQ=255)                        ggTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTTTGCCCTGAGAAGTGGGCGCAGAGAgg                                                                                                                               >  2:529991/1‑151 (MQ=255)                         gTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTCTGTCGTGCGAAGAGAGCGACGAGAggg                                                                                                                              >  1:820316/1‑151 (MQ=255)                              gagaGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCGTTTTTTCTTTTTTTTCCATTAGAAGGGAACGCagag                                                                                                                                  >  1:966922/1‑142 (MQ=255)                               agagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATc                                                                                                                        <  2:1329223/151‑1 (MQ=255)                                gagTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TGAATTTTTTTTTTTCTCACCTGGAAAGTGGACGCAGAGggggggtatt                                                                                                                        >  1:720694/1‑144 (MQ=255)                                 agTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCtt                                                                                                                      <  2:153020/151‑1 (MQ=255)                                  gTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑GCATTTTTTTTTTTTCTTCCATGCGAAGGGAGGGCAGAGAGGGGTTATCtt                                                                                                                      >  2:989606/1‑150 (MQ=255)                                    ttATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGc                                                                                                                    <  1:492650/150‑1 (MQ=255)                                      aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTCTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCct                                                                                                                  >  2:112578/1‑150 (MQ=255)                                      aTTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctg                                                                                                                 <  1:108455/151‑1 (MQ=255)                                         cAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCctgact                                                                                                              <  2:1100532/151‑1 (MQ=255)                                          aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCCTTTTTTTTTTTTTTTCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCctgactg                                                                                                             >  2:160911/1‑151 (MQ=255)                                          aGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGGGAGCGCAGAGAGGGGTTATTTTTCctgact                                                                                                              >  2:990489/1‑150 (MQ=255)                                          aGAAAGTTTGCGGGTAACCCGTAAAAGCTCCATACCTATACCAGTCAGGAGTCAATTCAGGATAACATGTCGAGGAAAACGTTTGCGGCGA‑TCATTTTTTCTTTTTCTGACCTGCGAAGTGAGCGCAAAGCGGGATTATCTTACctggcgg                                                                                                             <  1:374251/151‑5 (MQ=255)                                              aGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑CCGTTTTTTCTTTTTGTTCCATTGGAAAATTGGGCAAAAAGGGGTTATATTGCctatctgta                                                                                                           >  2:667841/1‑142 (MQ=255)                                                    cGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TTATTTTTTTTTTTTTTGCCTTGAGGAGTGGGCGCAGAGAGGGGTTCACTTGTCTGGCCGAACg‑ggg                                                                                                     >  2:1117217/1‑147 (MQ=255)                                                      ggTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTa                                                                                                <  2:1120831/151‑1 (MQ=255)                                                         aaCCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCCTTTTTTTTTTTTTTGCCCTGAGGAGTGGGCGCAGAGAGGGGTTATCTTGCCTGACTGAACT‑GGCGATTaa                                                                                               >  2:140592/1‑149 (MQ=255)                                                           cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTTAGAAGGGAGGGCAGAGAGGGATTTTCTTTCCCTGCTGAACT‑GGCGGTTaaagtt                                                                                           >  2:542281/1‑147 (MQ=255)                                                           cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGttt                                                                                           <  2:108845/151‑1 (MQ=255)                                                           cccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGttt                                                                                           <  2:139640/151‑1 (MQ=255)                                                            ccTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTc                                                                                          <  2:1245212/151‑1 (MQ=255)                                                               tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTTAGAAGGGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGGTTAAGttttc                                                                                         >  2:1079822/1‑147 (MQ=255)                                                               tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTTCCTGACTGAACT‑GGCGATTAAGttttc                                                                                         >  2:1193769/1‑147 (MQ=255)                                                               tataAGCTCCATACCTATACCACTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTcc                                                                                         <  1:536384/149‑1 (MQ=255)                                                                ataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTCTTCCATTAGGAGGTGGCGCCGAGGGGGGGTATTTTTCCTTAATTAACC‑GGCGATTTAGttt                                                                                           >  1:713906/1‑146 (MQ=255)                                                                  aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TTATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTTAACT‑GGCGATTTAGGTTTcgc                                                                                       >  1:201149/1‑148 (MQ=255)                                                                  aaGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTTTTCCCTGGGAAGTGAGAGCAGAGAGGGGTTATTTTGCCTGGCTTAACT‑GGCGATTAAATTTTcg                                                                                        >  1:898130/1‑147 (MQ=255)                                                                    gCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGCGCAGGGAGGGGTTATTTTTCCTGTCTGAACT‑GGCGATTAAGTTTTCGCGG‑Ac                                                                                  >  1:1185821/1‑150 (MQ=255)                                                                     cTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTCATTTTTT‑TTTTTCTGCCCTTTGAAGTGAGCGGTGAGGGGGGATTTTTTACCTGACTTGACT‑GGCGATTTAGattttccctaacc                                                                                 >  2:852905/1‑141 (MQ=255)                                                                       ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTTTGCCATGAGAAGTGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACT‑GGCGGTTAAGTTTTCGCGT‑Accccg                                                                              >  2:1383728/1‑149 (MQ=255)                                                                       ccATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCgg                                                                              <  2:1387194/151‑1 (MQ=255)                                                                         attccTATGCCAGTCAGGAGTTATTTCAAGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTGTCTTTATCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTCCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGc                                                                             <  1:1071879/147‑1 (MQ=255)                                                                           aCCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTTAGCAGTGAGCGGAGAGGGGGATTCTTCTGCCCTACTGAACG‑GGCGGCTTAGTTTCCGCGC‑Acccc                                                                               >  1:604857/1‑145 (MQ=255)                                                                              tataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCATGAGGAGTGAGGGCAGGGAGGGGTTATTTTGCCTGACTGAAAT‑GGCGGTTTAGTTTTTGCGG‑ACCCCGCGgtt                                                                         >  1:1034613/1‑149 (MQ=255)                                                                               ataCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTa                                                                        <  1:212415/149‑1 (MQ=255)                                                                                 aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTTCCCTGAGAAGGGAGCGCGGGGGGGGGTTTTTTTTCCCGACTGGAAT‑TGCGGTTAAATTTTCGCGG‑ACCCCGGGttttg                                                                       >  1:447347/1‑146 (MQ=255)                                                                                 aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGtt                                                                     <  1:221932/150‑1 (MQ=255)                                                                                 aCCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTa                                                                    <  2:1015764/151‑1 (MQ=255)                                                                                  ccAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTGAGt                                                                      <  1:885233/148‑1 (MQ=255)                                                                                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCgg                                                              <  1:1053987/151‑1 (MQ=255)                                                                                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCgg                                                              <  1:1140948/151‑1 (MQ=255)                                                                                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCgg                                                              <  1:415842/151‑1 (MQ=255)                                                                                       cAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCg                                                               <  1:304596/150‑1 (MQ=255)                                                                                          gAGTCATTTCAGGAGGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑TCCCGGCGTTTAGTTATTGCGgatg                                                           <  2:995152/151‑1 (MQ=255)                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTGCTTTTTTTCCATTTGAAATGAGCCGCGGGCGGGATTATTTTTCCTGACATAACC‑GGCGAATAAGATTTCGCGT‑ACCCGGGGTTTAGTTTTTTCGgaggat                                                         >  2:275840/1‑148 (MQ=255)                                                                                               aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑GCCTTTTTTCTTTTTTTGCCATGAGGAATGGGGGGCGAGAGGGAATATATTTGCTGTCTGAACTGGGGGTTTAAGTTT‑CGCGT‑TCCCCGCGgttttgtgttgggg                                                              >  1:708741/1‑132 (MQ=255)                                                                                                  tCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCCTGAGAAGTGAGGGCAGAGAGGGGTTATCTTGCCTGACTGAAAT‑GGCGGTTTAGTTTTCGCGT‑ACCCGGCGTTTTGTTATTTCGGGTGATGGTAtt                                                   >  2:643277/1‑151 (MQ=255)                                                                                                   cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCCTTTTTTTTTTTTTTGCCCTGAGAAGAGAGCGCAGAGAGGGGATATTTTTCCTGACTGAAAT‑GGGGATTAAGTTTCCGCGTAACCCGGGGTTT‑GTTATTTCGGGTGATGATAtt                                                   >  2:1365512/1‑150 (MQ=255)                                                                                                   cAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTg                                                  <  2:730852/151‑1 (MQ=255)                                                                                                     ggATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGaa                                                <  2:80933/151‑1 (MQ=255)                                                                                                         gCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTc                                            <  1:115151/151‑1 (MQ=255)                                                                                                          ccATGTCGCGGATCTCGTTTGCGGCGG‑TCCTTTTTTCTTTTTTTTCCATGAGAAGGGAGCGCAGAGAGGGGTTATCTTGCCTGACTGAACT‑GGCGATTTAGTTTTCGCGT‑ACCCGGCGGTTAGTTATTGCGGGTGATGATATTTAAGCTc                                            >  2:889077/1‑150 (MQ=255)                                                                                                              gTCGCGGATCTCGTTTGCAGCGG‑TCATTTTTTCTTTTTCTGCCAAGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGATCT‑GGCTATTAAGTTTCCGCGT‑ACCCGGCGGTTAGTTATTGCGGATGATGATATTGAAGCGCGGca                                        <  1:192099/150‑2 (MQ=255)                                                                                                                   ggATCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTTTGCCATGAGAAGGGAGGGCAGAGAGGGATTATCTTGCCTGACTGAAAT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTTGTTATTGCGGATGATGATATTGAAGCTCGGCTGttt                                    >  2:1024122/1‑149 (MQ=255)                                                                                                                     aTCTCGTTTGCGGCGG‑TCATTTTTTTTTTTTCTGCCCTGAGAAGTGAGCGCAGAGAGGGATTATCTTTCCTTACTGAACT‑GGCGGTTAAGTTTCCGCGT‑ACCCCGCGTTTAGTTATTTCGGGTGATGATATTGAAGCTCGGCTGATTTGtt                                >  2:352669/1‑151 (MQ=255)                                                                                                                         cGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGtt                                <  2:528712/147‑1 (MQ=255)                                                                                                                                  cgg‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATc                   >  1:97909/1‑151 (MQ=255)                                                                                                                                       cATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCAcc                <  2:316175/150‑1 (MQ=255)                                                                                                                                          tttttctttttcTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg            >  1:792737/1‑151 (MQ=255)                                                                                                                                           ttttctttttctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTg            >  2:1159059/1‑150 (MQ=255)                                                                                                                                                    tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa   <  1:824590/150‑1 (MQ=255)                                                                                                                                                    tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTa   <  2:1246500/150‑1 (MQ=255)                                                                                                                                                    tctGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAt  >  2:1200382/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                        AAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGG‑TCATTTTTTCTTTTTCTGCCATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACT‑GGCGATTAAGTTTCCGCGT‑ACCCGGCGTTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATTGGTTCATTGTCGCCATCACCGCTGGACGGTGTAT  >  NC_000913/4604203‑4604499

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