Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,508,657	0	A	G	40.4%	85.6 / 14.5	58	F386S (TTC→TCC)	yhfW	phosphopentomutase‑related metalloenzyme superfamily protein
Reads supporting (aligned to +/- strand):  ref base A (11/23);  new base G (23/0);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 5.59e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508513‑3508804                                                                                                                                                   |                                                                                                                                                    cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAAtt                                                                                                                                               <  2:648349/151‑1 (MQ=255) cACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGGGGGGCGGC‑GGaa                                                                                                                                                    >  2:371862/1‑147 (MQ=255)     ttgttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GCGGGAAAAAtta                                                                                                                                              >  1:78902/1‑148 (MQ=255)       gttTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGgaaaa                                                                                                                                                 >  1:848253/1‑143 (MQ=255)          tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATaacaa                                                                                                                                          >  2:662094/1‑146 (MQ=255)          tttAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAATTAACACACg                                                                                                                                       >  1:792564/1‑150 (MQ=255)                   ggAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCACACACGGGAGccccc                                                                                                                              >  1:865903/1‑150 (MQ=255)                      gagaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCATACACGGGAGccccc                                                                                                                              >  1:548098/1‑147 (MQ=255)                       agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATTACACACGGTAGCCCCcaca                                                                                                                           >  1:28538/1‑149 (MQ=255)                       agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATc                                                                                                                         <  1:193112/151‑1 (MQ=255)                       agaTCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATc                                                                                                                         <  2:315214/151‑1 (MQ=255)                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCACACACGGGAGCCCCCCCATCAg                                                                                                                       >  1:412093/1‑149 (MQ=255)                           cAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGaa                                                                                                                     <  2:784321/151‑1 (MQ=255)                            aaTATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGaaa                                                                                                                    <  1:397981/151‑1 (MQ=255)                                  cAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCACACACGGGAGCCCCCCCATCAGAAAGCgggg                                                                                                              >  1:576416/1‑151 (MQ=255)                                     aTGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTgc                                                                                                           <  1:882173/151‑1 (MQ=255)                                       gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGCAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGG‑GG‑GGGAAAAATCACACACGGTTGCCCCCCCATCAGGAAACGGGGGgcg                                                                                                          >  1:603335/1‑150 (MQ=255)                                       gCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAAATTACACACCGTAACCCCCCCATCAGAAAGCgtggt                                                                                                             >  1:719597/1‑147 (MQ=255)                                               gCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAg                                                                                                 <  2:274209/151‑1 (MQ=255)                                                 gAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATTACACACGGGAGCCCCCCCATCAGAAAGCGTGGGGCGCACAccca                                                                                                  >  2:924803/1‑146 (MQ=255)                                                         aaCCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCg                                                                                       <  1:23780/151‑1 (MQ=255)                                                              gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCACACACGGTAGCCCCCCCATCAGAAAACGTGGTGCGCACACCCCGCTGCGTAGCGatc                                                                                    >  1:869552/1‑149 (MQ=255)                                                              gTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAACGGGGTGGGCACACCCAGCTGCGTAGCGatca                                                                                   >  1:485717/1‑150 (MQ=255)                                                                  gATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCt                                                                              <  1:535259/151‑1 (MQ=255)                                                                   aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATAAAACACGGTAGCCCCCACATCAGAAAGCGGGGGGCGCGCACCGAGCTGCGTAGCGATCATccc                                                                               >  2:199714/1‑149 (MQ=255)                                                                   aTCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCtt                                                                             <  2:379488/151‑1 (MQ=255)                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCACACACAGTAGACCCCACACAAGAAAGCGGGGGGCGCACACCCAGCGGCGGAGCGATAATCCCCTGCt                                                                          >  2:560850/1‑151 (MQ=255)                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCACACACCGTAGCCCCCCCATCAGAAAACGTGGGGCGCACACCGAGCTGCGTAGCGATCATCCCTTGc                                                                           >  1:895693/1‑150 (MQ=255)                                                                      aaaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCCGAAAACGGGGGGGGCCCACCGAGCTGGGGAGCGATTATcccct                                                                             >  1:469167/1‑146 (MQ=255)                                                                       aaaCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTg                                                                         <  2:834685/151‑1 (MQ=255)                                                                         aCTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTGAt                                                                       >  1:233118/1‑151 (MQ=255)                                                                               aGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCa                                                                 <  1:868589/151‑1 (MQ=255)                                                                                 ggTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGc                                                               <  1:167248/151‑1 (MQ=255)                                                                                        ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTgg                                                          <  1:419512/149‑1 (MQ=255)                                                                                        ccTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTg                                                           <  1:894758/148‑1 (MQ=255)                                                                                         cTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGAGGGGG‑GCGGGAAAAATTAAAACCGGGAAGCCCCC‑CATCCGAAAGCGGGGGGGGCGCACCGCGCCGCGGAGGGATCCACCCCTTCTGATAAAACCGCACCGGcac                                                       >  2:81976/1‑151 (MQ=255)                                                                                          tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGAAAAATTCACAAACCGTAGCCCCCACAACAGAAAGGGCGGTGCGGACACCGGGCTGCGGAACGAAACTCCCTTGCTGTTTAAACAGCCCTTGca                                                        >  1:55023/1‑149 (MQ=255)                                                                                          tGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAAATTACACACGGTAGCCCCCACATCAGAAAACGGGGGGCGCACACCGAGCTGCGTAGCGATCATCCCCTGCTGATAAACCAGCACTGGgacc                                                      >  1:260104/1‑151 (MQ=255)                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacc                                                      <  1:328713/150‑1 (MQ=255)                                                                                           gCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGcacc                                                      <  2:329764/150‑1 (MQ=255)                                                                                               aCCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATCACACACGGGAGCCCCCACATCAGAAAGCGGGGTGCGCACACCCAGCTGCGTAGCGATCATCCCTTGCTTATAAAACAGCACTGGccccacct                                                  >  1:495119/1‑148 (MQ=255)                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGCGAAATTCACACACAGGAGCCCCCCCATCAGAAAGCAAGGGGCGGCCACCGCGCTGAGTAGCGGTCCTCCCTTTCTGATAAACCAGCACCGGcaccac                                                    >  2:333099/1‑146 (MQ=255)                                                                                                 cGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAATTTACAACCGGGAGACCCCCCATCACAAAGCGCGGGGGGCGCACCCCGGTGCGCAGCGATTCTCCCTTGCTTGTAAACCACCACCgg                                                          >  2:454672/1‑140 (MQ=255)                                                                                                  ggagggagGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAATTTACACAAGGTAGCCCCCACATCCGAAAGCGTGGGGGGCACACCCCGCTGCGGAGCGATCCTCCCTTTCTCGTAAAACAGCACTGGCACCACTTcgc                                               >  2:398319/1‑150 (MQ=255)                                                                                                       gaggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATTAAACACGGTAGCCCCCACACCAGAAAGCGGGGTGGGCACACCGAGCTTCGTAGCGATCCTCCCTTTCTGAACAAACAACACCGGCCCCACTTCGCggg                                            >  2:49371/1‑148 (MQ=255)                                                                                                         ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGGGC‑GC‑GGGAAAAATTACACACGGGAACCCCCACATCAGAAAGCGTGGGGCGCACACCGAGCTGCGGAGCGATCATCCCTTGCTGATAAAACAGCACCGGCACCACTTCGCGGGGAtgg                                       >  2:716942/1‑151 (MQ=255)                                                                                                         ggaAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTAtgg                                       <  2:797209/151‑1 (MQ=255)                                                                                                            aagaaaagaGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAATTTACACACGGGAGCCCCCACATCAGAAAACGTGGTGGGCACACCGAGCTGCGTAGCGATTATCCCTTGCTGATAAAACAACACCGGCACCACTTCTCGGGGAtggggg                                    >  1:875411/1‑151 (MQ=255)                                                                                                              gaaaagaGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGGAAAATTTACCCAAGGTAATCCCCCCCTCCGAAAGCGAGGAGGGCCCACCGCGCTGCGTAGCGACCATCCCTTGCTGATAAACCACCACCGGCACCCCCTCTCggg                                            >  1:432575/1‑141 (MQ=255)                                                                                                                    agCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAAATCACAAACGGTAGCCCCCACATCAGAAAAAGGGGGGGGCACACCCAGCTGCGTTTCGATCATCCCTTGCTGATAAAACAACACTGGCACCACTTCGCGGGGATGGTGGGTGTGGcc                            >  1:107961/1‑151 (MQ=255)                                                                                                                      cGACCATTTTGAGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCaa                          <  2:414357/151‑1 (MQ=255)                                                                                                                           aTTTTGTGGCGGTGGGGC‑GC‑GGGAAAATTCACAAACGGGAGCCCCCCCATCAAAAAGCGTGGGGGGCCCACCCAGCTGCGGAGCGATCATCCCCTTCTGATAAACCAGCACTGGCACCACTTCGCGGGGATGGTGGGTGTGGCCaa                          >  1:484467/1‑146 (MQ=255)                                                                                                                                 tggcggtggcgC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATc                <  1:36875/150‑1 (MQ=255)                                                                                                                                  ggcggtggcgC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCGCCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt              <  2:170377/151‑1 (MQ=255)                                                                                                                                  ggcggtggcgC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGt              >  2:306913/1‑151 (MQ=255)                                                                                                                                        ggcgC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCg         <  1:452586/150‑1 (MQ=255)                                                                                                                                         gcgc‑gc‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTg       <  2:174677/151‑1 (MQ=255)                                                                                                                                            c‑gc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc    >  2:357849/1‑151 (MQ=255)                                                                                                                                            c‑gc‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATc    >  1:776350/1‑151 (MQ=255)                                                                                                                                               c‑gGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCg  >  2:791576/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CACGTTGTTTTTAAAGCCGGAGAGATCAATATGCAGATGCTCATGGGCGAGGGTGTAACCCGTCGGATCAAAACTCATAGGGTGTCTCCTGCAAACCGGAGGGAGGAAAGAAAAGAGCGACCATTTTGTGGCGGTGGCGC‑GC‑GGAAAAATTCACACACGGTAGCCCCCACATCAGAAAGCGTGGTGCGCACACCGAGCTGCGTAGCGATCATCCCTTGCTGATAAACCAGCACTGGCACCACTTCGCGGGTATGGTGGCTGTGACCAATGGTCGGATCGTTGCCGTGATCCG  >  NC_000913/3508513‑3508804

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