Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,718	0	G	C	37.2%	125.3 / 30.0	79	P292R (CCG→CGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (15/34);  new base C (29/0);  total (45/34)
Fisher's exact test for biased strand distribution p-value = 1.12e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT  >  NC_000913/2984569‑2984867                                                                                                                                                      |                                                                                                                                                       caccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg                                                                                                                                                       <  1:586224/151‑1 (MQ=255) caccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCgg                                                                                                                                                       <  1:1420989/151‑1 (MQ=255)   ccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgcg                                                                                                                                                     <  1:13241/151‑1 (MQ=255)        gacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCaact                                                                                                                                                >  2:830036/1‑149 (MQ=255)         acattgacCCGAGGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTc                                                                                                                                               <  1:1443146/151‑1 (MQ=255)           attgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGTGCATACAACCgcgcccatccc                                                                                                                                             >  2:67184/1‑144 (MQ=255)            ttgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCACAAACAGCCGCAATTctcc                                                                                                                                            >  2:618437/1‑148 (MQ=255)              gacCCGACTCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTCCACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCAt                                                                                                                                           <  1:878011/150‑1 (MQ=255)                    aCGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCATTCCCCTCAATAcc                                                                                                                                    >  1:723673/1‑151 (MQ=255)                    aCGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTTCCATCAATGcc                                                                                                                                    >  1:448616/1‑151 (MQ=255)                      gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAACCCGCGGAAATCCCAttcaagccg                                                                                                                                   >  1:1420828/1‑142 (MQ=255)                          aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCt                                                                                                                              <  1:969648/151‑1 (MQ=255)                            cccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑a                                                                                                                            <  2:623765/150‑1 (MQ=255)                            cccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCACACACCCGCGCAAATCCCCCCCATGCCGGACCtcc                                                                                                                            >  2:36157/1‑149 (MQ=255)                             ccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCAGCCCCATTACCACCAAATCCGGATCt                                                                                                                              >  1:1329387/1‑148 (MQ=255)                              cACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAAACCGCGCCAATCCCCCCAATGCCGGCGCT‑acac                                                                                                                         >  2:1172382/1‑151 (MQ=255)                                  cAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACAACCGCGCAATTCCCCTCCATGCCGGAGCC‑ACACCAAc                                                                                                                     >  2:26870/1‑151 (MQ=255)                                    aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGc                                                                                                                   <  1:1344408/151‑1 (MQ=255)                                      cATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGGATACACCCGCGCCATTCCCAACCATGCCGGATCT‑Acccccaacc                                                                                                                   >  1:1408998/1‑144 (MQ=255)                                          ttGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCAATGCCGGGTCT‑ACACCAACGCTgg                                                                                                                >  1:100701/1‑148 (MQ=255)                                           tGCCAACCGAAAGCGCCTGTGCAGCAAGCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGc                                                                                                               <  1:1337314/148‑1 (MQ=255)                                           tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCGGCGCAATTCCCATCAAAGCCCGATCT‑ACACCAACCCTgg                                                                                                                >  1:610772/1‑147 (MQ=255)                                           tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCa                                                                                                              <  1:495491/149‑1 (MQ=255)                                           tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAAATCCCCTCAAAGCCGGGACT‑ACAcccaacctgg                                                                                                                >  2:246487/1‑139 (MQ=255)                                           tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCGCAATTTCCATCAAAGCCGGATCT‑ACACCCACGCTGGCaa                                                                                                             >  1:461652/1‑150 (MQ=255)                                             ccAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAAt                                                                                                           <  1:1202373/150‑1 (MQ=255)                                             ccAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGCCGCATACACCCGCGCAATTACCATCAATGCCGGGTCT‑TCACCAACGCTGGCaaa                                                                                                            >  2:1202009/1‑149 (MQ=255)                                              cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGGTCT‑ACACCAACGCTGGCAAAAgc                                                                                                         >  1:456056/1‑151 (MQ=255)                                                  cGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGTGCATAAAACCCCGCCATTCCCCACAATGCCCGATCC‑ACCCCCACACCGGCAAAtac                                                                                                         >  2:1301535/1‑145 (MQ=255)                                                   gAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTTGCCTACACCCCCGCCCTTCCCCCCAATGCCGGGTCT‑CCACCAACGCCGGCTAAAGCGCGa                                                                                                     >  2:284775/1‑150 (MQ=255)                                                      aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCCAAACCCGCGCAATTCCCATCAATGCCGGGTCT‑ACACCCACCCTGGCAAAAGcgcgga                                                                                                    >  2:77881/1‑146 (MQ=255)                                                       gcgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCg                                                                                                <  1:365981/151‑1 (MQ=255)                                                        cgcCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCGGCGCAAATCACCTCAATGACCGGGCT‑ACACCAACGCCGGCAAAAGCGCCAAAGCgg                                                                                               >  2:1385349/1‑151 (MQ=255)                                                           cTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCCCGCCCATCCCCCCCATGCCCGAACG‑TGCCCCACACCGGCAAAAGCGCGAAAGCggggc                                                                                            >  2:903968/1‑149 (MQ=255)                                                            tgtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGGTCT‑ACACCAACGCTGGCAAAAGCGCGAATGCggg                                                                                              >  1:611336/1‑148 (MQ=255)                                                            tgtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAAATCCCCTCACTGCCGGATCT‑ACACCAACGCTGGCAAATTCGCGAAAGCggggc                                                                                            >  2:606765/1‑148 (MQ=255)                                                             gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCATCAATGCCGGATCT‑ACAACCACGCTGGCAAAAGCGCGAATGGGGGcc                                                                                            >  1:845958/1‑149 (MQ=255)                                                             gtgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACGTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCa                                                                                           <  2:752855/150‑1 (MQ=255)                                                                cagcaAACGCTTCATTAGCCGCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATATTCCGGCGCAATTCCCCGCCATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCAc                                                                                       <  2:1215487/151‑1 (MQ=255)                                                                   caAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACggg                                                                                    <  1:98008/151‑1 (MQ=255)                                                                     aaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCAATTCCCCCCCATGCCGGGTCT‑TCACCAACGCTGGCCAAAACGCCAATGCCGGCCCGCACggg                                                                                    >  2:654304/1‑149 (MQ=255)                                                                      aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCAATTCCCATCCATGCCCGGTCT‑ACACCCAAGCCGGCAAAAGGGCGAAAGCGGGCCAGCACgggg                                                                                   >  2:767968/1‑148 (MQ=255)                                                                       cGCTTCATTAGCCGCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCGCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTaaa                                                                                <  2:415995/151‑1 (MQ=255)                                                                        gCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAAt                                                                               <  1:467486/151‑1 (MQ=255)                                                                             aTTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTAAACCGCCCCCATTTCCCTCAATTCCGGGACT‑ACCCCAAAGCCGGCAAAATTGCGCATGCCGGCCCCCGCCGGGAAATTcca                                                                           >  1:180985/1‑148 (MQ=255)                                                                                aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCGGCGCAATTCCCATCCATGCCGGATCT‑CCACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATTCAACGcc                                                                       >  2:910478/1‑150 (MQ=255)                                                                                aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTAGCAAAAACCCGCGCCATTCCCATCACACCCGGAACT‑ACACCAACACTGGCAAAAGCGCGACAGCGGGCCAGCAAGGGGaaaattaaa                                                                          >  2:1346989/1‑142 (MQ=255)                                                                                   cTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGGTCAAACCCGCGCCATTCCCATCAATGCCGGATCTAACACCAA‑GCTGGCAAAAGCGCGAAAGCGGGCCCGCCCCGGGAAATTCAACGc                                                                        >  2:1309356/1‑147 (MQ=255)                                                                                     cGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTg                                                                   <  1:472266/150‑1 (MQ=255)                                                                                      gATGAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCGCCATTCCCCTCCATGCCGGGTCT‑ACACCAACGCTGGCCAAAACGCGAATGCGGGCCCGCACCGGGAAATTTAACACCCGTg                                                                   >  2:86319/1‑149 (MQ=255)                                                                                      gATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCACACCCCGCCGCAATTCCCCCCCATGCCCGATCC‑ACACCAACGCCTGCAAAAGCGCGAAAGCGGGCCCGCCCCGGTAAATACAACGCTCg                                                                     >  2:78114/1‑147 (MQ=255)                                                                                       aTAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCtt                                                                <  1:284784/151‑1 (MQ=255)                                                                                                  ttCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCTCATACACCGGCGCAATTCCCCTCAAAGCCGGAGCT‑CCACCAACGCTGGCAAAAGCGCGGATGCGGGGCAGCACGGGGAAATTCAACGCTCGTGCGTTTGCTTCGc                                                       >  2:665911/1‑149 (MQ=255)                                                                                                   tCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtc                                                     <  1:312368/150‑1 (MQ=255)                                                                                                       ccAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCCAAACCCGCGCCATTCCCACCAATTCCGGGTAC‑ACCCCAAACCTGGGAAATGCGGGAATGCCGGGCCGCCCGGGGGAAATTAACCCCCCTGGTTTTGCTTccccccaca                                                 >  2:1120311/1‑142 (MQ=255)                                                                                                        cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatc                                               <  1:980063/151‑1 (MQ=255)                                                                                                        cAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTAACCCCGCGCCAATCCCCCCAATGCCGGGTCT‑ACACCCACGCCGGCGAAAGCGCGGAAGCGGGCCCGCACCGGGAAAATCCACCCCCGCGCTTTTGCTTCGCccctcatc                                               >  2:20967/1‑151 (MQ=255)                                                                                                         aaCTGCCAGCCAACCCGCTCCAGGCAACGGCGGGTCGCATACTCCGGAACGATTCCCCGCACTCCGGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatcatca                                              <  2:28838/151‑1 (MQ=255)                                                                                                          aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCGTACAACCGCGCAAATCCCCTCCATGCCGGGTCT‑TCACCAAACCTCGGAAAAACGCGAAAGCGGGGCCGCACGGGGAAATTCAAAGCCCGGGCTTTTGCTTTGCtcactctca                                              >  2:193126/1‑144 (MQ=255)                                                                                                          aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATAAACCCGCGCCATTCCCATTCAAGTCGGGTCT‑CCACCAACGCTGGCCAATGCGCGCAACCGGGCCAGGACGGGGAAATTTAACACTCGGGCTTTGGCCTCCCtcatcatc                                               >  2:400987/1‑149 (MQ=255)                                                                                                             gCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCAtt                                            <  1:1388705/149‑1 (MQ=255)                                                                                                               cAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCg                                        <  1:1253970/151‑1 (MQ=255)                                                                                                                aGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCGCCATTCCCATCAATGCCGGAACT‑ACACCAACGCTGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAATTCAACGCTCCTGCCTTGGCTTCGCTCCACATCATTa                                           >  1:1105720/1‑147 (MQ=255)                                                                                                                 gCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCa                                      <  1:88038/151‑1 (MQ=255)                                                                                                                    tACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGGTCT‑ACACCCACGCTGGCAAATGCGCGAATGGGGGCCAGCACGGGGAAAATCAACGCCCCGGCTTTTGCTTTGCCCATCATCATTACCGc                                       >  1:140329/1‑147 (MQ=255)                                                                                                                    tACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCGCCATTCCCCTCAAAGCCGGATCT‑ACACCAAAGCTGGCAAATGCGCGGAAGCGGGCCCGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCTCATCATCATTTCCGc                                       >  1:53060/1‑147 (MQ=255)                                                                                                                     acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTg                                  <  1:163257/151‑1 (MQ=255)                                                                                                                        cGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcgc                               <  1:980116/151‑1 (MQ=255)                                                                                                                        cGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgcgc                               <  1:1303330/151‑1 (MQ=255)                                                                                                                              aGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGt                         <  2:142598/151‑1 (MQ=255)                                                                                                                                gCAACGGCGGGTCGCATTCTCCGGAGCGATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGttt                       <  2:603403/151‑1 (MQ=255)                                                                                                                                gCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGttt                       <  1:1143190/151‑1 (MQ=255)                                                                                                                                  aaCGGCGGGTCGCATACACCGGCGCGATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTa                      <  2:933643/150‑1 (MQ=255)                                                                                                                                  aaCGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTAt                     <  1:749948/151‑1 (MQ=255)                                                                                                                                    cggccgGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATgg                   >  1:12189/1‑151 (MQ=255)                                                                                                                                      gcggGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat                 <  1:580369/151‑1 (MQ=255)                                                                                                                                          gTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGatgatg             <  1:1347394/151‑1 (MQ=255)                                                                                                                                            cGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCa           <  1:819556/151‑1 (MQ=255)                                                                                                                                              cATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt         >  1:538013/1‑151 (MQ=255)                                                                                                                                              cATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAt          >  2:119358/1‑150 (MQ=255)                                                                                                                                                     cGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt  >  2:571022/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCGCAATTCCCATCAATGCCGGATCT‑ACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT  >  NC_000913/2984569‑2984867

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