Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,662	0	T	A	43.5%	99.8 / 23.7	62	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (5/30);  new base A (27/0);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 7.35e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATC  >  NC_000913/3508519‑3508802                                                                                                                                                 |                                                                                                                                             gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAAttaaa                                                                                                                                          >  1:262253/1‑146 (MQ=255) gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTcacaca                                                                                                                                       <  1:698143/151‑1 (MQ=255)   ttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTacacac                                                                                                                                      >  1:650087/1‑150 (MQ=255)    tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATcacaca                                                                                                                                       >  2:734499/1‑148 (MQ=255)    tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATca                                                                                                                                           >  2:949784/1‑144 (MQ=255)      tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGt                                                                                                                                   <  2:636639/150‑1 (MQ=255)       aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCaaacaccg                                                                                                                                    >  2:1129059/1‑146 (MQ=255)         aGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGTAGcc                                                                                                                               >  2:642632/1‑151 (MQ=255)           ccGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGCGGGAAAAAttaaacaagggaagcc                                                                                                                               >  1:349511/1‑136 (MQ=255)             ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGGAGccccc                                                                                                                            >  2:27560/1‑150 (MQ=255)               agagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACGGTAGccccccc                                                                                                                          >  2:548139/1‑148 (MQ=255)                     cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAATTCACAAACGGTAGCCCCCCCAACAGa                                                                                                                    >  1:465970/1‑150 (MQ=255)                       atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATTACACACGGGAGCCCCCCCAACAGaaa                                                                                                                  >  2:568508/1‑150 (MQ=255)                         atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGGAGCCCCCCCATCAGaaa                                                                                                                  >  2:742170/1‑148 (MQ=255)                         atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCg                                                                                                               <  2:374556/151‑1 (MQ=255)                               aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAAATCACAAACGGGAGCCCCCACCTCCGaaaaaggggggg                                                                                                         >  1:510959/1‑142 (MQ=255)                                tGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGC‑GGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAACgcg                                                                                                             >  1:705589/1‑144 (MQ=255)                                 gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcgc                                                                                                       <  1:452535/151‑1 (MQ=255)                                   tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAAAAACCGGAACCCCCACCTCAGAAAGCggggggggga                                                                                                      >  1:48723/1‑144 (MQ=255)                                    cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                    <  2:845531/151‑1 (MQ=255)                                    cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGca                                                                                                      <  2:258207/149‑1 (MQ=255)                                    cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGATAGCCCCCACATCAGAAAGCGTGGTGCGcaca                                                                                                    <  1:424956/151‑1 (MQ=255)                                       gggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCg                                                                                                 <  1:883891/151‑1 (MQ=255)                                        ggCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGa                                                                                                <  1:480319/151‑1 (MQ=255)                                            aGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCGAGCt                                                                                             >  2:1119484/1‑150 (MQ=255)                                                 gtAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACAGCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAg                                                                                       <  1:13392/151‑1 (MQ=255)                                                    aCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCGCCCACATCAGAAAGCGTGGTGCTCACACCGAGGTTCGTAGCGa                                                                                    <  1:879469/151‑1 (MQ=255)                                                       cGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGatca                                                                                 <  1:284536/151‑1 (MQ=255)                                                        gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCAAAAACGGGAGCCCCCCCATCCGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGatc                                                                                  >  2:734853/1‑149 (MQ=255)                                                          cGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATcc                                                                              <  1:748285/151‑1 (MQ=255)                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGG‑GGAAAAAATTACACACGGGAGCCCCCACATCAGAAAACGGGGGGCGCACACCGAGCCGCGTAGCGATCATccc                                                                             >  1:1021392/1‑146 (MQ=255)                                                                aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCt                                                                        <  1:322321/151‑1 (MQ=255)                                                                 aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                       <  1:979739/151‑1 (MQ=255)                                                                  aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGGC‑GGAAAAAATCACACACCGGAGGCCCCACATCAGAAAACGGGGGGCGCACACCCAGGCGCGTAGCGGTCATCCCCTGc                                                                         >  2:778634/1‑148 (MQ=255)                                                                  aaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAAAGGGGTGGGCACACCCAGCTGCGTAGCGATCAAccc                                                                             >  1:757823/1‑144 (MQ=255)                                                                     tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATCACACACCGTAGCCCCCCCACCAGAAAACGGGGGGTGCACACCCAGCCGCGTAGCCATCATcccctac                                                                         >  2:335253/1‑141 (MQ=255)                                                                     tCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGGAAAATTAAACACGGGAACCCCCCCATCAGAAAACGTGGTGCGGACACCGCGGTGCGGAGCGATCATCCCTTGCTGATa                                                                    >  1:363349/1‑150 (MQ=255)                                                                      cATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATaa                                                                   <  2:120294/150‑1 (MQ=255)                                                                       aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACCGTAGCCCCCCCATCAGAAAACGGGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGGTaaa                                                                  >  2:570403/1‑150 (MQ=255)                                                                       aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAAc                                                                 <  2:590910/151‑1 (MQ=255)                                                                           ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                             <  1:9161/151‑1 (MQ=255)                                                                                       aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCACACACGGGAGCCCCCCCAACCGAAAACGTGGGGGGCACACCCCGCTGCGGAGGGATCATCCCCTTCTGGTAAAACAGCACTGGcaccac                                                  >  2:755789/1‑150 (MQ=255)                                                                                       aaaCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATAACAAACGGGGGGCCCCCCATCAGGAAGAGGGGGGCGCACACCCAGCGGCGGTGCGGTCATCCCCTGCGGGTAAAACAACCCAGGccccacc                                                 >  1:85409/1‑150 (MQ=255)                                                                                         aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                               <  1:527654/151‑1 (MQ=255)                                                                                         aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTc                                               <  1:1108236/151‑1 (MQ=255)                                                                                            ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGC‑GGAAAAAAAAAAAACAGGGAGCCCCCCCACCAAAAAGCGTGGTGCGGCCACCCAGCCTCGTGAGGCACATCCCCTTCCGATAAACACCCCCTGGCACCACTTccc                                             >  1:1078838/1‑148 (MQ=255)                                                                                             gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTAAACACGGGAGCCCCCCCACCACAAAACGGGGGGCGCACACCCCGCTGCGGAGCGATCATCCCCTTCTGATAAAACAGCACTGGCACCACCTcgc                                             >  1:1136855/1‑149 (MQ=255)                                                                                             gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                           <  2:261873/151‑1 (MQ=255)                                                                                               gggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCggg                                          <  1:145077/150‑1 (MQ=255)                                                                                                    gaaagaaaAGAGCGACCATTTTGTGGCGGACGCGCGC‑GGGAAAATTTAAACACAGGAAGCCCCTCAACAAGAAAAGGGGGGCGCCCAACCCGCTGCCGAGCGGAACTACCCTGCTGATAAAACAGCACCCGCCCCACTTccgggggctg                                      >  1:422847/1‑140 (MQ=255)                                                                                                       agaaaagaGCGACCATTTTGTGGCGGTGGGGCGC‑GGGAAAAATAACACACGGTAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCGAGCCGCGTAGCGGTCATCCCCTGCTGTTAAAACAGCACTGGCACCACTTCGCGGGGAtggtgg                                  >  1:240584/1‑150 (MQ=255)                                                                                                       agaaaagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATCAAAAACCGGAGCCCCCACATCAGAAAGAGGGGTGGGCACACCCAGCTGCGGAGAGAACATCCCTTTCTGATAAAACAGCACTGGCACCACCTCCCGGGGAtggtgg                                  >  2:524467/1‑150 (MQ=255)                                                                                                           aagaGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCtgtg                             <  2:1084809/151‑1 (MQ=255)                                                                                                             gagCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGAc                           <  1:909438/151‑1 (MQ=255)                                                                                                                  aCCATTTTGTGGCGGTGGCGCGC‑GGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAGCGGGGTGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACCTCGCGGGGATGGTGGCTGTGAcc                          >  2:337409/1‑147 (MQ=255)                                                                                                                   ccATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                     <  2:950706/151‑1 (MQ=255)                                                                                                                            ggcggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt            <  2:250436/151‑1 (MQ=255)                                                                                                                             gcggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGtt           <  2:274147/151‑1 (MQ=255)                                                                                                                              cggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg          <  2:1068678/151‑1 (MQ=255)                                                                                                                               ggtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGc         >  2:1110699/1‑151 (MQ=255)                                                                                                                                gtggcgCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGcc        <  2:924478/151‑1 (MQ=255)                                                                                                                                    cgcgc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg     >  1:537716/1‑150 (MQ=255)                                                                                                                                      cgc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc  <  2:612339/151‑1 (MQ=255)                                                                                                                                        c‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc  >  2:640104/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                             GTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCGC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATC  >  NC_000913/3508519‑3508802

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