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.9%	83.6 / 34.3	66	E384D (GAA→GAT)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (9/28);  new base A (29/0);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 3.72e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAT  >  NC_000913/3508518‑3508808                                                                                                                                                   |                                                                                                                                                   tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGG‑CGGGAAAAAttaaa                                                                                                                                                >  1:1132320/1‑147 (MQ=255) tgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAAATcaaa                                                                                                                                               >  2:716774/1‑148 (MQ=255)   tttttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCAACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATcacaca                                                                                                                                             >  2:852480/1‑150 (MQ=255)       tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACgg                                                                                                                                          <  1:258349/149‑1 (MQ=255)       tAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGt                                                                                                                                         <  1:719665/150‑1 (MQ=255)        aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTa                                                                                                                                        <  2:609867/150‑1 (MQ=255)        aaaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAg                                                                                                                                       <  2:386976/151‑1 (MQ=255)         aaGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCACAAACGGt                                                                                                                                         >  1:729800/1‑148 (MQ=255)              ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATTACACACGGGAGccccc                                                                                                                                  >  2:952005/1‑150 (MQ=255)               gagagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCcac                                                                                                                                <  2:208146/151‑1 (MQ=255)                   gaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCa                                                                                                                            <  1:1197623/151‑1 (MQ=255)                    aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAACACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGG‑C‑GGGCAAAATTACCCACAGTCGGCCCCCCAACa                                                                                                                            >  2:703164/1‑150 (MQ=255)                    aTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATTACAAACGGGAGccccccc                                                                                                                                >  2:573231/1‑144 (MQ=255)                       aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAATTTACACACGGTAGCCCCCCCATCaaaa                                                                                                                         >  2:478469/1‑150 (MQ=255)                        atatGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCACAAACGGGAGCCCCCCCCTCAGcaa                                                                                                                        >  1:485691/1‑150 (MQ=255)                          atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                      <  1:101970/150‑1 (MQ=255)                          atGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGc                                                                                                                      <  1:57481/150‑1 (MQ=255)                                  gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgcg                                                                                                              <  1:286206/150‑1 (MQ=255)                                    tCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAATTTACACACGGGAGCCCCCCCATCAGAAAACgtgggggg                                                                                                              >  2:1161258/1‑144 (MQ=255)                                     cATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGAGGGGGGGGcac                                                                                                           >  2:808937/1‑150 (MQ=255)                                      aTGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGG‑C‑GGGAAAAATTACACACGGTAGCCCCCCCATCAGAAAACGCGGGgcgcgc                                                                                                           >  1:982292/1‑147 (MQ=255)                                       tGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCAAACACGGGAGCCCCCACATCAGAAAACGGGGGGGGGACAcc                                                                                                        >  2:96865/1‑151 (MQ=255)                                           cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                    <  2:429779/151‑1 (MQ=255)                                           cGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGc                                                                                                    <  1:154327/151‑1 (MQ=255)                                               ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGGGCGCACACCGAGCTg                                                                                                  >  2:113218/1‑149 (MQ=255)                                               ggTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCACACACCGTAGCCCCCACATCAGAAAGCGTGGGGCGCACACCCAGCTGc                                                                                                 >  2:1178490/1‑150 (MQ=255)                                                      cccGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAAATTACAAAACGTAGCCCCCACCACAGAAAACGTGGTGCGCACCCCGGGCTg                                                                                                  >  2:721341/1‑142 (MQ=255)                                                         gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGGG‑C‑GGGAAAAATTAAACACAGTAGCCCCCACACCAGAAAACGTGGGGCGCCCACCCAgcggcggggcg                                                                                           >  1:618310/1‑141 (MQ=255)                                                          tCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGCGGCGCG‑C‑GGGAAAAATTAAACACGGGAGCCCCCCCACCAAAAAACGGGGGGGGCCCACCCAGGCGCGGAGCGatc                                                                                        >  1:1196059/1‑148 (MQ=255)                                                              aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGGGGGGGGCCCACCCAGCCGCGGAGCGGTCATcccct                                                                                 >  1:393395/1‑149 (MQ=255)                                                              aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                                 <  1:524566/151‑1 (MQ=255)                                                              aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGCGCG‑C‑GGGAAAAATTAAACACGGGAGCCCCCCCCACAGAAAACGGGGTGGGCACACCCAGCCGCGTAGCGATCAAccc                                                                                   >  1:1035118/1‑149 (MQ=255)                                                                cAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                               <  2:677719/151‑1 (MQ=255)                                                                 aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAATTCACACACGGGAGCCCCCCCATCAGAAAGCGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTTCt                                                                              >  2:760983/1‑151 (MQ=255)                                                                     cTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATTACACAACGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCCGCGGAGCGATCATCCCTTTCTGAt                                                                           >  1:1043559/1‑150 (MQ=255)                                                                        aTAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑G‑GGGAAAAATTAAACACGGGAGCCCCCCCAACAGAAAACGGGGGGGGGGCAACCAGCTGCGGAGCGAGCATCCCCTTCTGATaaa                                                                        >  1:38048/1‑150 (MQ=255)                                                                          aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCACACACGGGAGCCCCCCCCTCAAAAAACGTGGTGGGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATaaaac                                                                      >  1:485623/1‑148 (MQ=255)                                                                          aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                     <  1:737835/151‑1 (MQ=255)                                                                              tgtCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGGG‑C‑GGGAAAAATTACACAACGGAGCCCCCCCACCAGAAAAGGGGGGGCGCACACCCAGGTGCGGAGCGATCAACCCTTGCTGATaaa                                                                        >  1:1171905/1‑144 (MQ=255)                                                                                tctcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAATTTAAACACGGGAGCCCCCCAATAAGAAAACGTGGTGGGCCCACCGAGCTGCGTAGCGATCACCCCTTGCTGTTAAAACa                                                                     >  2:108362/1‑145 (MQ=255)                                                                                  tcCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAAATCACACACGGGAGCCCCCCCATCAGAAAACGGGGGGCGCACACCCAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTgg                                                              >  2:1203341/1‑150 (MQ=255)                                                                                   ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACCGGc                                                             >  2:340319/1‑150 (MQ=255)                                                                                   ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                              <  1:33938/149‑1 (MQ=255)                                                                                   ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                               <  2:233266/148‑1 (MQ=255)                                                                                              gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                 <  2:309081/151‑1 (MQ=255)                                                                                              gagggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCgg                                                 <  1:197749/151‑1 (MQ=255)                                                                                               agggaggAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAACGTGGGGGGCACACCCAGCTGCGTAGCGATCATCCCTTTCTGATAAACCAACACTGGCACCACTTCGCggg                                                >  1:1112258/1‑151 (MQ=255)                                                                                                 ggaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCGCC‑GGAAAAATT‑ACACACCGGAGCCCCCCCATACGAAAGAGTGGGGGGCACACCCCGCCGCGGAGCGATCATCCCCTTCTGGTAAAACAACACTGGgacccccttccgggg                                               >  1:134578/1‑139 (MQ=255)                                                                                                   aggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAAATTACAAACCGGAGCCCCCCCATCAGAAAACGGGGGGGGCACACCCGACTGCGTGGCGGTCATCCCCTGCTGATAAAACAGCACCGGCACCACCTCGCggg                                                >  1:592205/1‑147 (MQ=255)                                                                                                     gaaagaaaAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtggt                                          <  1:1090975/151‑1 (MQ=255)                                                                                                      aaagaaaagaGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAAATAAAACACGGTAGGCCCCCCATTAGAAAACATGGGGCGCACACCAAGATGTGGAGCGAATATGTCTTGCTGGTTAACGACAACCGGCACCCCCTCGCGGGTAtgggg                                         >  2:624908/1‑149 (MQ=255)                                                                                                       aagaaaagaGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGGAAAAATTACACACGGGAACCCCCCCATCAGAAAACGGGGTGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAACACTGGCACCCCCTCGCGGGGAtggggg                                        >  2:609851/1‑151 (MQ=255)                                                                                                        agaaaagaGCGACCATTTTGTGGCGGTGGGGCG‑C‑GGGAAAAATTACACACGGTAGCCCCCACATCAGAAAACGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAAACAGCACTGGCACCACTTCGCGGGGAtggggg                                        >  2:39926/1‑150 (MQ=255)                                                                                                                    ccATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATgg                           <  2:276374/151‑1 (MQ=255)                                                                                                                        tttGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCg                        <  2:1115736/150‑1 (MQ=255)                                                                                                                         ttGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCgg                       <  2:894108/150‑1 (MQ=255)                                                                                                                          tgtgGCGGTGGGGCG‑C‑GGAAAAAATCACACACCGGAGCCCCCCCATCAGAAAACGGGGTGGGCACAACCAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACTGGCACCACCTCGCGGGGATGGTGGGTGTGACCAATGGTCgg                       >  1:663218/1‑149 (MQ=255)                                                                                                                           gtggcggtggcgCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGTTGTGACCAATGGTCGGATc                    <  2:986414/151‑1 (MQ=255)                                                                                                                               cggtggcgCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTg                <  1:156049/151‑1 (MQ=255)                                                                                                                                   ggcgCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAATCCAGCACTGGCAACACTTCGCGGGTATGGTGGCTGTGACCAATGGGCGGATCGTTGCCGt            <  1:777475/151‑1 (MQ=255)                                                                                                                                   ggcgCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg             <  1:654602/150‑1 (MQ=255)                                                                                                                                     cgcg‑c‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg           >  2:530265/1‑150 (MQ=255)                                                                                                                                          c‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg      >  2:1168696/1‑151 (MQ=255)                                                                                                                                            ggAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGc     <  2:29960/151‑1 (MQ=255)                                                                                                                                             gAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc    >  2:171749/1‑151 (MQ=255)                                                                                                                                              aaaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGcc    >  2:555071/1‑150 (MQ=255)                                                                                                                                               aaaaTTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAt  <  1:459023/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   TGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGCG‑C‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCGCCAT  >  NC_000913/3508518‑3508808

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