Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,378,802	0	A	C	41.9%	110.0 / 15.4	62	V48G (GTA→GGA)	ampC	penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (32/4);  new base C (0/26);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.31e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378655‑4378952                                                                                                                                                        |                                                                                                                                                       ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accg                                                                                                                                                     <  2:356359/151‑1 (MQ=255)   tCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgcc                                                                                                                                                   >  1:734694/1‑151 (MQ=255)    cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgcca                                                                                                                                                  <  1:823631/151‑1 (MQ=255)      aGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccacc                                                                                                                                                >  1:761267/1‑151 (MQ=255)         aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGTGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccacc                                                                                                                                                <  1:468519/148‑1 (MQ=255)         aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccaccgcc                                                                                                                                             >  1:452582/1‑151 (MQ=255)          aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccaccgcc                                                                                                                                             >  2:738187/1‑150 (MQ=255)              cGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATAc                                                                                                                                         >  2:510494/1‑150 (MQ=255)                tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCg                                                                                                                                      >  1:1054103/1‑151 (MQ=255)                  gCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCggg                                                                                                                                    >  2:1004265/1‑151 (MQ=255)                     accaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATc                                                                                                                                 >  1:214939/1‑151 (MQ=255)                      ccaAGCACGCCAGTAAATTTTTTGCTGACCGAACCTAACTAAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTCCCC‑TG‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATc                                                                                                                                 <  1:881841/150‑1 (MQ=255)                         aGCCCGCAAGTAATTTTTTTGCGGACCGACCCTAACCAAAACACCTTTTGTTTTGTGAGGGGCTGTTTTTTGGCTATGCCCGCATAGCCCCAGATAAAGTAATAAGTTTTCCCC‑TA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCtt                                                                                                                               <  2:827778/149‑1 (MQ=255)                               caagaaagtGTTTTGCACACCGACCAAAACAAAAACAACGTTTGGTGAGGGAGGGGTTTTTTTGGGGCGATGTCCCAAGCGCCCAGGAAAAAATAATATGTTTTACCC‑TG‑ATAAA‑TA‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                        <  2:1006852/142‑1 (MQ=255)                                   tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGT‑AAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                   >  2:904794/1‑150 (MQ=255)                                    aaaTGTTTTGCTGCCCGAACCTAACAAAAACAACTTTTGTTGTGTGAGGGGTTGTTTTTGGGCGATGCCCGCATACCCCAAGGTAAAGTAATATTTTTTCCCC‑TA‑ATAAA‑TT‑ACCCCCACCGCCATACCCGGGATCTTttgttg                                                                                                                         <  2:607332/144‑1 (MQ=255)                                         ttttGCTGCCCGACCCTAACTAAAACAACGTTTGCTGTGTGCCGGGTTGTTTTTTGGCGATGCCCCCATAGCCCCAGGAAAAGTAATAAGTTTTCCCC‑TA‑AAAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGc                                                                                                               <  2:820533/149‑1 (MQ=255)                                              cTGCCCGAACCAAACACAAACAGCTTTTGCTTGGGGAGGGGCTCTTTTTTGGCATTGCCCGCAACCCCCCAGAAAAAATAAAAAGTTTTCCCT‑TA‑AAAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                           <  2:828623/148‑1 (MQ=255)                                               tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                       >  1:1011319/1‑151 (MQ=255)                                                  ccGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                    >  1:1060625/1‑151 (MQ=255)                                                       ccTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                               >  1:103607/1‑151 (MQ=255)                                                        cTAACCCAAACAACGTTTGCGGTGTGGGGGGCTTTTTTTTGGTGCTGTCCTCATCCCCCCGGTAAAAGTAATAAGGTTTCCCA‑TG‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgt                                                                                                    <  1:695351/145‑1 (MQ=255)                                                          aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                             >  2:655493/1‑150 (MQ=255)                                                            cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                           >  1:894835/1‑150 (MQ=255)                                                            atCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                           >  2:233307/2‑150 (MQ=255)                                                              cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                        >  1:359332/1‑151 (MQ=255)                                                                 aaactgTTTGCTGTGGGTCGGTTTTTTTTTTGGGGTTGCCCGCATCGCCCAAGGTAAAGTAAT‑AGGTTTTCCC‑TT‑ATAAA‑TTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                     <  2:730562/146‑1 (MQ=255)                                                                   aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                                   >  1:458073/1‑151 (MQ=255)                                                                       tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                               >  2:393211/1‑151 (MQ=255)                                                                               gtgAGGGGCTGTTTTTTGGCGATGCCCGCATAGCCCCAGGTAAAGTAAAAAGTTTTCCCC‑TG‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                            <  1:869121/146‑1 (MQ=255)                                                                               gtgACGGGCTGTTTTTTGGGGATGCCCGCAAACCCCAAGGAAAAGAAATAGGTTTTCCCT‑TA‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                       <  1:894923/151‑1 (MQ=255)                                                                               gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         >  1:383499/1‑149 (MQ=255)                                                                                tgccggtttttttttgGGCGATGCCCGCAT‑CCCCCAAGAAAAATAAATAGGTTTTCCCTTA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                         <  2:1066402/137‑1 (MQ=255)                                                                                  aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                    >  2:615515/1‑151 (MQ=255)                                                                                    gggCTTTTTTTTGGCGTTGCCCGAATACCCCAAGTTAAAGAAATAAGTTTTCCCT‑TA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgca                                                                  <  2:13783/151‑1 (MQ=255)                                                                                    gggCTGTTTTTGGGCGATCCCCGCATACCCCCAGTAAAAGTAATAGGGTTTCCCC‑TA‑ATAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAggg                                                                     <  1:453131/148‑1 (MQ=255)                                                                                     ggttgttttttGGCGATGCCCGCAAAGCCCGAGGTAAATTAATATGTTTTCCCT‑TA‑AAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                                 <  2:960436/145‑1 (MQ=255)                                                                                     ggttgtttttgGGGGATGCCCGCATACCCCAGGGAAAAGAAATAGGGTTTCCCC‑TA‑AAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                                 <  1:817363/145‑1 (MQ=255)                                                                                         gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                              >  1:778532/1‑150 (MQ=255)                                                                                         gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                              >  1:250241/1‑150 (MQ=255)                                                                                           ttttgGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGTTTTCCCC‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaaa                                                             <  1:238165/149‑1 (MQ=255)                                                                                               tGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                        >  1:280753/1‑150 (MQ=255)                                                                                                   gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                    >  1:110477/1‑150 (MQ=255)                                                                                                        ccacATACCCCAAGTAAAAGTAATAATTTTTCCCT‑TA‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                 <  2:67164/145‑1 (MQ=255)                                                                                                        ccGCATACCCCCAGGAAAAGTAATAAGGTTTCCCC‑TG‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                                  <  1:191431/147‑1 (MQ=255)                                                                                                         cgaaTAGCCCAAGAAAAATAAATAAGGTTTCCCC‑TA‑ATAAA‑TT‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                               <  1:818158/146‑1 (MQ=255)                                                                                                            agaccccggggTAAAGTAATAAGTTTTCCCC‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                           <  1:62866/141‑1 (MQ=255)                                                                                                              aGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                        >  2:963313/1‑151 (MQ=255)                                                                                                                cccaAGGTAAAGAAAAAAGGTTTCCCC‑TG‑ATAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                         <  1:49151/148‑1 (MQ=255)                                                                                                                 ccaggatAAAGTAAT‑AGGTTTTCCC‑TT‑AAAAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                     <  2:683065/145‑1 (MQ=255)                                                                                                                     ggTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                    >  2:604997/1‑148 (MQ=255)                                                                                                                      gTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                >  1:28082/1‑151 (MQ=255)                                                                                                                        aaaGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCGCAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                              >  1:161860/1‑151 (MQ=255)                                                                                                                         aaGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:675325/1‑151 (MQ=255)                                                                                                                         aaGAAATAGGGTTCCCCC‑GA‑AAAAA‑TT‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             <  1:446581/151‑1 (MQ=255)                                                                                                                           gTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGGAATTAATAAGgc                               >  2:450937/1‑147 (MQ=255)                                                                                                                            taataaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                             >  1:232203/1‑148 (MQ=255)                                                                                                                              ataaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGATGCTCTATAAGCGGGGTAATTGTGCGATGCGCAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAAGTAATAAGgcgc                             >  1:105224/1‑146 (MQ=255)                                                                                                                              aaaaggttttcctatG‑AAAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                          <  2:214893/136‑1 (MQ=255)                                                                                                                                aaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                      >  1:343905/1‑151 (MQ=255)                                                                                                                                    tttACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                  >  1:473359/1‑151 (MQ=255)                                                                                                                                    gttttcccataaTAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                      <  2:857158/138‑1 (MQ=255)                                                                                                                                     ttcccc‑TA‑AAAAA‑TC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                  <  2:1012130/147‑1 (MQ=255)                                                                                                                                      ttccc‑TGAATAAA‑TC‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                  <  1:560092/148‑1 (MQ=255)                                                                                                                                       tacccttaaTAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtt                  <  2:894701/142‑1 (MQ=255)                                                                                                                                        ccc‑TG‑AAAAT‑TT‑ACGCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                <  1:758155/149‑1 (MQ=255)                                                                                                                                        ccc‑TA‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt                <  1:557625/149‑1 (MQ=255)                                                                                                                                         cc‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGGTTTGa              >  1:786178/1‑150 (MQ=255)                                                                                                                                                  aa‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg      >  2:649259/1‑151 (MQ=255)                                                                                                                                                        accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt     >  2:702163/1‑148 (MQ=255)                                                                                                                                                        accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg  >  2:53698/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTG  >  NC_000913/4378655‑4378952

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