Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,895,032	0	T	C	38.9%	120.3 / 20.3	73	I250V (ATT→GTT)	ygcQ	putative flavoprotein
Reads supporting (aligned to +/- strand):  ref base T (38/6);  new base C (1/27);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 5.75e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCC  >  NC_000913/2894883‑2895181                                                                                                                                                      |                                                                                                                                                      aGGCCAGACCCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAAt                                                                                                                                                       >  2:824839/1‑150 (MQ=255)    ccAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGcc                                                                                                                                                   >  2:679222/1‑151 (MQ=255)     cAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGcca                                                                                                                                                  <  1:122256/151‑1 (MQ=255)     cAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGcc                                                                                                                                                   >  1:91409/1‑150 (MQ=255)         aCCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCcaccac                                                                                                                                              >  2:261095/1‑151 (MQ=255)           cGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACaa                                                                                                                                            <  1:940966/151‑1 (MQ=255)                  aGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCt                                                                                                                                     <  1:1182127/151‑1 (MQ=255)                   ggATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTa                                                                                                                                    <  1:1242510/151‑1 (MQ=255)                        atcCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCg                                                                                                                               >  1:324947/1‑151 (MQ=255)                        atcCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCg                                                                                                                               <  1:727039/151‑1 (MQ=255)                          ccGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGt                                                                                                                              >  1:878680/1‑150 (MQ=255)                               aaCACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcgccc                                                                                                                        >  2:887729/1‑151 (MQ=255)                                  acGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTAcgcccgcc                                                                                                                     >  1:531034/1‑151 (MQ=255)                                    gTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCAt                                                                                                                   >  1:811956/1‑151 (MQ=255)                                      tcttGGCAATCTGCGTTAATGATTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAGCATCTGCCTCTGAAAATACTGCAGCGCTGGCGTCATGATTAACCGCCACCACAATCGTGCTATTGCGCACGCCCGCCAt                                                                                                                   >  2:877793/3‑149 (MQ=255)                                                cTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCtgctgc                                                                                                       >  2:683066/1‑151 (MQ=255)                                                cTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCtgctgc                                                                                                       >  2:401091/1‑151 (MQ=255)                                                 tGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTgcg                                                                                                      >  1:815413/1‑151 (MQ=255)                                                   cGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGcc                                                                                                    >  1:34390/1‑151 (MQ=255)                                                      gAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTgcgc                                                                                                     >  2:410851/1‑147 (MQ=255)                                                         tttttGTCCCCGGTCCCCCCACCCCCACTCTCCAGTCATCACCCCCCCCAACTTCTCCGTGTAAAAATCCCGCGGCGCGGGCGTAATATTTAACCCCCCCCAAAACTTTGTTTTTGCGTCCGCCCGCCATTAACGCTGCTgcgc                                                                                                     <  1:511003/142‑1 (MQ=255)                                                           tttGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAgcgc                                                                                            >  1:953881/1‑151 (MQ=255)                                                           tttGCCCCCAGCCCCCCCAGCCCCACTTCCCACACAACACCCACCCAACCATCTCCCTGAAAAAACCCCGCCGCGCGGGCGTTATTTTAACCCCCCCCCACAAACTTGTTTTTGCGTCCCCCCGCCATTAACGCTGCTGCGCCAGAAg                                                                                               <  2:1298774/148‑1 (MQ=255)                                                            ttGTCCCCATCGCCCCCAGCCCCCCTTCCCAGTCCCCACCCCCGCCAACATCTGCCTGGGAAAACACCGCAGCGCTGGCGTAATGATTAACCCCCACCAAAACTTTGGTATGGCTTCCGCCCGCCATTAACGCTGCTGCGCCa                                                                                                   <  2:1244176/143‑1 (MQ=255)                                                              gTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCca                                                                                          >  1:871105/1‑150 (MQ=255)                                                               tCAACAGTGCCTCCAGCACCACTTTTCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCGGCGCTGGCTTCATGATTAATCGCCACCTTAAACTTGCTATTGCGTACGCCCGCCATTAACGCTTCCGCCCCAGa                                                                                                 >  2:587220/1‑142 (MQ=255)                                                                ccccAGTCCCCCCACCCCCACTTCCCACCCATCACCCACCCCAACATCTGCCTGAAAAAATCCCGCAGCGGTGGCGTCATGATTAACCCCCACCAAAAACTTGCTATGGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCcaaca                                                                                       <  1:477897/149‑1 (MQ=255)                                                                cACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCcaac                                                                                        >  1:1217198/1‑150 (MQ=255)                                                                cACCAGCCCCCCCACCCCCTCTTCCCAGCCACCACCCCCCCCACCATCTGCCGGAAAAAACCCCGCGGCGCTGGCGTCATATTTAACCCCCACCAAAACCTTTCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCcaa                                                                                         <  1:828819/149‑1 (MQ=255)                                                                        ccTCCAGCCCCCCTTCCCACTCACCACCCCCCCAACCATCGCCCTGAAAAAACCCCGCGGCGGGGGGGTCATGATAAACCCCCACCAAAACCTTGCTTTTGCGTCCACCCGCCATTAACGCTGCTGCGCCAGAAGCGCcaacaa                                                                                      <  1:241023/144‑1 (MQ=255)                                                                          cccacccccacTTCCCAGTCACCACCCCCCCAAACCTCTCCCTATAAAAACACCGCGGCGCTGGCACAATGATTACCCCCCCCCACAAGCTTGCTATGGCCTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATacaca                                                                             <  2:732725/144‑1 (MQ=255)                                                                           ccacccccccTTCCCAGTAACCACCCCCGCCAACTTCTGCCTGTAAAAACCCCGCGGCGCGGGCGTAATGATTAACCCCCACCAAAACCTTGCTATTGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATacaca                                                                             <  1:1341845/141‑1 (MQ=255)                                                                           ccAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATacacac                                                                            >  1:1022999/1‑151 (MQ=255)                                                                               ccccacTTCCCACTCACCACCCCCCCCAACATCTCCTTGTAAAAACCCCGCAGCGCTGGCGTCATGATTAACCCCCACCACAACCTTGTTTTTGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCa                                                                        <  1:1271635/149‑1 (MQ=255)                                                                               caccacTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACGAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTc                                                                         >  2:1259408/1‑150 (MQ=255)                                                                                 ccacTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAgg                                                                      >  1:949033/1‑151 (MQ=255)                                                                                     ttgccacacatcaACCACCCCACCTTCTGCTTGTAAAAATCCCGCGGCGCGGCCGTAATAATAAGCCCCCACCACAAACTTTCTATTGCGACCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGt                                                                     <  1:1328626/140‑1 (MQ=255)                                                                                      ttCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGcc                                                                  >  2:934060/1‑150 (MQ=255)                                                                                          aGTCATCAACCACGCCAACATCTACCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCcagc                                                               >  1:1356464/1‑149 (MQ=255)                                                                                             catcaACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATg                                                          >  1:937921/1‑151 (MQ=255)                                                                                             caccaACCACCCCAACTTCTCCCGATAAAAGTCCCGCGGCGCTGGCTTCATTATTAACCCCCACCACAAATTTTTTAGTGACTACCCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGa                                                            <  2:716708/146‑1 (MQ=255)                                                                                                   ccccGCCACCCTCCGCCTGTAAAACCCCCTCAGCGCTGGCGTCATAATACGACCCCCCCACAAACTTGCTTTTCCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGAt                                                           <  2:682603/141‑1 (MQ=255)                                                                                                   ccACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAACCTTGCTATTGCGTACGCCCGCCATTCACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCgg                                                     >  2:568227/1‑150 (MQ=255)                                                                                                   ccACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCgg                                                     >  1:1312382/1‑150 (MQ=255)                                                                                                    ccccccAACACCTGCCTGGGAAATCCCCGCAGCGCGGGCTCAATAATCAACCCCCCCCACAAACTTCATATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGaa                                                   <  2:262409/147‑1 (MQ=255)                                                                                                    cACCCCAACTTCTCCCTGAAAAAATCCCGCGGCGTTGGCGTAATGATTAACCCCCCCCAAAAATTTGTTTTGGCGTACCCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGaa                                                   <  1:659195/151‑1 (MQ=255)                                                                                                      ccccAACCTCGGCCTGAAAAAACCCCGCCGCGCTGGCGTCATGATTACCCCCCACCACAAACTTGTTATGGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGaaa                                                  <  1:1227296/148‑1 (MQ=255)                                                                                                      cGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAAt                                                 >  2:1252441/1‑151 (MQ=255)                                                                                                        ccAACCTCTGCCTATAAAAATCCCGCGGCGCGGGCGTCATAATTACCCCCCCCCACAAACTTTTTATGGCGCCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATa                                                <  2:177597/150‑1 (MQ=255)                                                                                                            ccgcTGCCTGTAAAAACCCCGCACCGCTGGCTTCATGATAAACCCCCCCCACAAACTTGCTATTGCCTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATa                                                <  2:177915/143‑1 (MQ=255)                                                                                                                  ccTGTAAAAACCCCGCAGCGCTGGCTTCATGATAAACCCCCACCACAACCTTGCTATTGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCa                                         <  1:1090086/147‑1 (MQ=255)                                                                                                                     gtgAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTc                                  >  1:892171/1‑151 (MQ=255)                                                                                                                      tgAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTc                                  >  2:877643/1‑150 (MQ=255)                                                                                                                      aaaaaaccccccAGCGCGGGCCTCTTAATAATTCGCCCCCACAACCTTTCTCTGGCCTCCCCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCa                                         <  1:520640/132‑1 (MQ=255)                                                                                                                       gAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTc                                  >  1:1419/1‑149 (MQ=255)                                                                                                                       aaaaaccccGCGGCGCTGGCTTCATGATTAACCGCCACCACAAACTTGCTTTTGCGCCCCCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCa                                 <  1:975847/143‑1 (MQ=255)                                                                                                                       aaaaaccccGCGGCGCGGGCGTAATAATTACCCCCCCCCACAACCTTTTTATTGCTTCCCCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTc                                  <  1:1202903/142‑1 (MQ=255)                                                                                                                       aaaaaTCCCGCGGCGCTGGCGTCATAATTACCCCCCACCACAAACTGGCTATTGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCa                                 <  1:709300/149‑1 (MQ=255)                                                                                                                       aaaaaTCCCGCAGGGTTGGCGTAATATTAAACCCCCACAAAAAACTTTCTATTGCGTACCCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTc                                  <  1:1325571/148‑1 (MQ=255)                                                                                                                        aaacccccGCGGCGCGGGCGTAATGATAAACCCCCACCAAAACCTTGCTTTGGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACtttt                                    <  1:1282116/140‑1 (MQ=255)                                                                                                                        aaaaaaCCCCGGCGCGGGCATCATAAATACCCCCCCCAAAAAACTTGTTAGTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCa                                 <  2:1150771/149‑1 (MQ=255)                                                                                                                               cGCAGCGGGGGCGTAATTATAAACCCCCCCCAAAATCTTGTTTTTGCTCCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCaa                          <  1:679402/149‑1 (MQ=255)                                                                                                                                gcggcggTGGCTTCATAATTAACCCCCACACCAACCTTTCTATTGCCACCGCCCGTCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACg                        <  2:324859/147‑1 (MQ=255)                                                                                                                                 cAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGcc                      >  2:1126362/1‑151 (MQ=255)                                                                                                                                  aGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGcc                      >  1:1126919/1‑150 (MQ=255)                                                                                                                                  aGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGcc                      >  1:842805/1‑150 (MQ=255)                                                                                                                                    cgcTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCAcc                   >  2:511225/1‑151 (MQ=255)                                                                                                                                     gggggCGTAATGATAACCCCCCCCCACACCCTTGTTATGGCGCCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCa                     <  1:261158/145‑1 (MQ=255)                                                                                                                                     gcTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCg                  >  1:465001/1‑151 (MQ=255)                                                                                                                                        ggCGTCATAATTAACCCCCACCAAAATCTTGTTTTGGCGTCCGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTtc               <  1:401188/151‑1 (MQ=255)                                                                                                                                             cATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCAc          >  1:449855/1‑151 (MQ=255)                                                                                                                                                  atAATCGCCACAATAAACTTGCTATTGCGTACGACCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTcgc      >  2:1275911/2‑150 (MQ=255)                                                                                                                                                   tAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCgcgc    >  1:843872/1‑151 (MQ=255)                                                                                                                                                     aTCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGccc  >  1:91267/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AGGCCAGAACCGTAGGCAGGATCATCCGGCAACACGTTATTGGCAATCTGCGTGAATGTTTGTCACCAGTGCCTCCAGCACCACTTTCCAGTCATCAACCACGCCAACATCTGCCTGTGAAAATACCGCAGCGCTGGCGTCATGATTAATCGCCACCACAAACTTGCTATTGCGTACGCCCGCCATTAACGCTGCTGCGCCAGAAGCGCCAACAACAATACACACCTCAGGTGCCAGCAGATGTCCGGAAATACCAATCACTTTTTCAGCATCAACGCCACCGTTCATCACTCGCGCCC  >  NC_000913/2894883‑2895181

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