Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,718	0	G	C	45.7%	92.5 / 38.3	71	P292R (CCG→CGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (6/32);  new base C (32/0);  total (38/33)
Fisher's exact test for biased strand distribution p-value = 3.57e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCC  >  NC_000913/2984568‑2984864                                                                                                                                                       |                                                                                                                                                   ccaccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCg                                                                                                                                                    <  1:984799/151‑1 (MQ=255)  caccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg                                                                                                                                                   <  1:425969/151‑1 (MQ=255)      attgacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCa                                                                                                                                               <  1:214038/151‑1 (MQ=255)      attgacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAACCCgcgc                                                                                                                                                >  1:323940/1‑150 (MQ=255)         gacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAAtt                                                                                                                                            <  1:303623/151‑1 (MQ=255)          acattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTc                                                                                                                                           <  1:1076335/151‑1 (MQ=255)             ttgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTccc                                                                                                                                         <  1:829915/150‑1 (MQ=255)               gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCAt                                                                                                                                       <  2:106149/150‑1 (MQ=255)                 cccGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCaa                                                                                                                                    <  2:1081410/151‑1 (MQ=255)                   cGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATg                                                                                                                                  <  2:816857/151‑1 (MQ=255)                      cGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCg                                                                                                                               <  1:789067/151‑1 (MQ=255)                       gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAATTTCCCTCAATCCCgg                                                                                                                              >  1:736152/1‑151 (MQ=255)                       gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCCCAACTcccccccatacc                                                                                                                                >  2:839349/1‑140 (MQ=255)                          cATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGa                                                                                                                             <  1:1027667/149‑1 (MQ=255)                              ccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTaca                                                                                                                       <  1:866783/151‑1 (MQ=255)                               cACTCAAGCATATTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCGCAAACCCGCGCAATTCCCCTCAATGCCGGACCTac                                                                                                                        >  2:530634/1‑149 (MQ=255)                                aCTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTAAACCGCGCCATTTCCCTCAATGCCCGGTCTACAcc                                                                                                                     >  2:205004/1‑151 (MQ=255)                                aCTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTACAcc                                                                                                                     >  1:705498/1‑151 (MQ=255)                                   cgaGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAACCCGCGCAATTCCCCTCAATACCGGGTCTTcaacaa                                                                                                                   >  1:679867/3‑150 (MQ=255)                                   cAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCATTCCCATCCAATCCGGATCTACACCaa                                                                                                                   >  1:679876/1‑150 (MQ=255)                                     aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTCCCATCAATGCCCGATCTACACCaa                                                                                                                   >  1:727417/1‑148 (MQ=255)                                      gCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCCGGAATAATTCGGAGCCATGCCAGATCTACACCAACGCt                                                                                                               <  2:416185/151‑1 (MQ=255)                                          cTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCAATGCCGGATCTACACCAACGCTgg                                                                                                             >  1:144983/1‑149 (MQ=255)                                             gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGTAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAAt                                                                                                        <  2:611079/151‑1 (MQ=255)                                             gCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAAt                                                                                                        <  1:1046593/151‑1 (MQ=255)                                               cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCTCCTGAACGATTCCCATCAATGCGAGATCTACACCAACGCTGGCAAATgc                                                                                                      <  2:458938/151‑1 (MQ=255)                                                 aCCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTTTCTGGCATGATTCGCCGCCCTCCGAGATCTACACCAACGCTGGCAAATgcgc                                                                                                    <  2:1074282/151‑1 (MQ=255)                                                  ccGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGcgcg                                                                                                   <  1:238594/151‑1 (MQ=255)                                                      aaGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCCTCCATGCCGGAACTACACCAAAGCTGGCAAAAGCGCGAATg                                                                                               >  1:244793/1‑151 (MQ=255)                                                       aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGc                                                                                              <  1:1094303/151‑1 (MQ=255)                                                         cgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCgg                                                                                            >  1:193970/1‑151 (MQ=255)                                                              gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAg                                                                                       <  1:303306/151‑1 (MQ=255)                                                                gcagcaAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCCATGCCGGATCTACACCAACGCCGGCCAATGCGCGAATGCGGGcc                                                                                         >  2:983780/1‑147 (MQ=255)                                                                    caAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACAACCGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACggg                                                                                 >  2:943493/1‑151 (MQ=255)                                                                     aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAAAGCCGGGTCTACAACAAAGCTGGCAAAAGCGCGAATGCGGGGCAGCACcgg                                                                                 >  1:1035368/1‑150 (MQ=255)                                                                      aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCAATGCCGGGTCTACACCAACGCTGGCAAAAGCGCGAATGCGGGGCAGCACggg                                                                                 >  1:842879/1‑149 (MQ=255)                                                                        cGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACTCCGGCACAATACGGCGCACAGCGAGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTaaa                                                                             <  2:343478/151‑1 (MQ=255)                                                                        cGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCATCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGaa                                                                              >  1:961419/1‑150 (MQ=255)                                                                        cGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCAAATCCCATCAATGCCCGATCTACACCAACCCTGGCAAAAGCGCGAAAGCGGGCCCGCCCGGGGaa                                                                              >  2:164312/1‑150 (MQ=255)                                                                           ttCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCAATCCCCCCAATGCCCGCTCTACACCAAACCCGGCAAAAGCGCGAAAGCGGGCCAGCCCGGGGaaa                                                                             >  2:633314/1‑148 (MQ=255)                                                                                tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCCCCTAAAACCCCGCAATTCCCCCCCATTCCGGAGCCACCACAACCCCGGCAAATGCCCCCAATCGGGCGCGCCCGGGGCAATTCAAccc                                                                     >  2:508015/1‑149 (MQ=255)                                                                                 aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCt                                                                    <  1:771324/151‑1 (MQ=255)                                                                                   ccTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCg                                                                  <  1:47816/151‑1 (MQ=255)                                                                                     tCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCCATTTCCCTCAATGCCGGGTCTACACCAACCCCGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAAATTAACCCTc                                                                   >  2:739425/1‑148 (MQ=255)                                                                                      cGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTg                                                                <  1:210645/150‑1 (MQ=255)                                                                                        aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCAGCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt                                                             <  2:505812/151‑1 (MQ=255)                                                                                             aTCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTAAAACCGCGCAATTCCCCTCCATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAAAGCGGGGCAGCACGGGGAAAATCAAAGCCCGTGCttt                                                            >  1:4236/1‑147 (MQ=255)                                                                                              tCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCCATTTCCATCAATGCCGGATCTACACCAACGCTGGCCAAAGCGCGAATGCGGGCCCGCACGGGGAAATTCAACGCTCGTGCTTTTGCtt                                                       >  1:773977/1‑151 (MQ=255)                                                                                                  cTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTTGCCTACCCCCGCCCCAATCCCCTCAATGCCCGGTCTACACCCACCCTGGCCAATGCGCGGATGCCGGCCCGCACGGGGAAATTCAACCCTCGTGCTTTTGCtt                                                       >  2:965175/1‑147 (MQ=255)                                                                                                   ttCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTTGCATACACCCGCGCAATTCCCCTCAAAGCCCGATCTACCCCAACGCCGGCAAAAGCGCGAACGCGGGCCCGCCCGGGGAAATTCAAAGCCCGGGCttt                                                            >  2:130806/1‑141 (MQ=255)                                                                                                       gCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTTTCGTAGACCCGCGCAATTCCCCCCACTAGCGCATCTACACCAACCCTGGCCAAAGCGCGAATACCGGGCCGCGCGGGGCAAATTAACGCTCGCGATTTGTCTTCGCtcgt                                                >  1:8814/1‑147 (MQ=255)                                                                                                        ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcat                                             <  1:704430/151‑1 (MQ=255)                                                                                                         cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCCATGCCCGATCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCCGCACGGGGAAATTCAACCCCCGTGCTTTGGCtttccccatc                                               >  1:813528/1‑140 (MQ=255)                                                                                                                    cTACACGCTCCAGGCAACGGCGGGGCGCATACCACCGCGCAATTCCCCTCCATTCCCGGTCTACACCCACGCCGGCCAAAGCGCGAATGCGGGCCCGCACCGGGAAATTTAACCCTCGTGCTTTTGCTTTCCTCATCATCATTACCGCAGc                                 >  2:969486/1‑151 (MQ=255)                                                                                                                    cTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCAATTCCCCTCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTTAAAGCCCGTGCTTTTGCTTCGCTCATCATCATTACCGc                                    >  1:41547/1‑148 (MQ=255)                                                                                                                     tACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGATCTACACCAAACCTGGCAAAAGCGCGAATGCGGGCCAGCACCGGGAAATTCAAAGCCCGTGCTTTGGCTTTCCTCCTCATCATTAcccc                                    >  1:296489/1‑145 (MQ=255)                                                                                                                      acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGc                                 <  1:403480/149‑1 (MQ=255)                                                                                                                      acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCt                                <  1:533514/150‑1 (MQ=255)                                                                                                                      acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg                               <  1:158595/151‑1 (MQ=255)                                                                                                                       cacGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATTCCGGATCTACACCAACCCTGGCAAATTCGCGAATTCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCTCATCATCATTACCGCCGCTg                               >  1:176835/1‑150 (MQ=255)                                                                                                                        acGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCAATCCCCACCATGCCCGATCTACACCAACGCTGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAATTCAACGCTCCTGCCTTTGCCTCCCCCCTCATCCTTAACCCCGCTgcg                             >  2:725396/1‑151 (MQ=255)                                                                                                                             ccAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGTTTCGCTCATCATCATTACCGCAGCTGCGCCATc                        <  1:983025/151‑1 (MQ=255)                                                                                                                               aGGCAACGGCGGGTCGCATACACCCGCGCCATTTCCCTCAATGCCGGATCTACACCCACGCCGGCAAATGCGCGAATGCGGGCCAGCACCGGGAAAATTAACGCCCGTGCTTTTGCTTCGCTCATCATCATTACCGCCGCTGCGCCCTc                        >  2:860658/1‑149 (MQ=255)                                                                                                                                     cggcggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATgg                <  1:300872/151‑1 (MQ=255)                                                                                                                                       gcggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat              >  1:433824/1‑151 (MQ=255)                                                                                                                                         gggTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatg             >  2:751772/1‑150 (MQ=255)                                                                                                                                         gggTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatg             >  2:981232/1‑150 (MQ=255)                                                                                                                                           gTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGatgatg          <  2:650826/151‑1 (MQ=255)                                                                                                                                                aTACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt      >  1:489226/1‑150 (MQ=255)                                                                                                                                                  acacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc   <  2:736125/151‑1 (MQ=255)                                                                                                                                                  acacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc   <  1:780075/151‑1 (MQ=255)                                                                                                                                                   cacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAcc   <  2:676215/150‑1 (MQ=255)                                                                                                                                                   cacCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGCGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTAccc  >  2:319679/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                   CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCC  >  NC_000913/2984568‑2984864

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