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

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948                                                                                                                                                           |                                                                                                                                                   ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑a                                                                                                                                                    <  1:806410/150‑1 (MQ=255)      cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgcca                                                                                                                                              <  2:867671/151‑1 (MQ=255)      cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgc                                                                                                                                                <  2:448839/149‑1 (MQ=255)       gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccac                                                                                                                                             >  1:1001051/1‑151 (MQ=255)       gAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑CT‑accgccac                                                                                                                                             >  2:1273242/1‑151 (MQ=255)         tcAATAGCGTCGCCACCAAGCACGATAGAAAATGGTTTGCCGACAGAACCTAAGTCAAACAACGTTTGCTGGGTTACGGGTTGATTTATTGCGCGGTCCGCATATCCGCAAGTTAAGT‑ATTAAGCTCTACCC‑TG‑ATAAA‑TT‑atagcagc                                                                                                                                             >  1:77328/2‑142 (MQ=255)          cAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccaccgc                                                                                                                                          >  2:1127126/1‑151 (MQ=255)           aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccaccgcc                                                                                                                                         >  2:1078/1‑151 (MQ=255)            aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑accgccaccgcc                                                                                                                                         >  1:907200/1‑150 (MQ=255)                   cGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCgg                                                                                                                                 >  2:777367/1‑151 (MQ=255)                      caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGAt                                                                                                                              >  2:770063/1‑151 (MQ=255)                      caccaAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGAt                                                                                                                              >  2:970991/1‑151 (MQ=255)                          aggcccGCCAGTAAATTTTTTGCTGCCCAACCCTAACAAAAACACCTTTTGTTGTGTGACGGGCTGTTTTTTGGCGATGCCCGCACCCCCCCAGGAAAAGT‑AAAAGGTTTTCCCC‑TA‑ATAAA‑TC‑CCCGCCACCGCCATACCCGGGAt                                                                                                                              <  1:158302/142‑1 (MQ=255)                          aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCttt                                                                                                                          >  1:336550/1‑151 (MQ=255)                           accccccctgTAAATTTTTTTCTGCCCAACCCCAACACAAACAACGTTTGCTGTGTGACGGTTTTTTTTTGGGCGATGCCCCCATCCCCCCGGGAAAATT‑AAAAGGTTTTCCCC‑TA‑AAAAA‑TC‑CCCCCCACCGCCATACCCGGGATCtt                                                                                                                           <  2:1035162/140‑1 (MQ=255)                              cccccAGAAATTTTTTTCCGGCCCGACCCTACCCCAAACAATGTTTGCTGTGGGAGGGTTTTTTTTTTGGCTATCCCCACACCCCCCCAGGAAAAAT‑AAAAGGTTTTCCCC‑TG‑AAAAT‑TC‑CCCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                      <  1:995467/148‑1 (MQ=255)                               cGCCAGTAATTTTTTTGCTGACCGAACCAAACAAAAACACCGTTTGCTGTGTGACGGGTTGTTTTTTGGCGATGTCCGAATACCCCCAGGTAAAGT‑AATAGGTTTTCCCT‑TA‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                       <  2:31858/149‑1 (MQ=255)                                gccggaaaaTTTTTTGCTGCCCGAACATACCACAACAACCGTTTGCGGTGTGACGGGCTGTTTTTGGGCGTGGCCCGCATGGCCACAGGAAAAGT‑AAAAAGGTTTCCCC‑AA‑ATAAA‑TT‑CCCGCCCCCGCCATACCCGGGATCTTttgttg                                                                                                                     <  2:180385/144‑1 (MQ=255)                                 ccAGTAATTTTTTTGCTGCCCAACCCTAACCCAAACACCGTTTGCTGTGGGCGGGGCTGTTTTTTGGGGTTGCCCGCATACCCCAGGGTAAAGT‑AAAAAGTTTTCCCC‑TG‑AAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTttgt                                                                                                                       <  2:112799/147‑1 (MQ=255)                                     tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                               >  1:779949/1‑151 (MQ=255)                                      aattGTTTTGCTGCCCGACCCTAACAAAAACAACGTTTGCTGTGGGAGGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGGAAAAGT‑AATAAGTTTTCCCT‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 <  2:72665/145‑1 (MQ=255)                                      aaaTTTTTTGCTGCCCGACCCAACCTAAACCACCGTTTGCTGTGTGACGGGTTGTTTTTTGGCGTTGCCCGCATAGCCCCAGGAAAATT‑AATAGGTTTTCCCT‑TA‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGctct                                                                                                                 <  2:517817/148‑1 (MQ=255)                                           ttttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgg                                                                                                         >  2:1171458/1‑151 (MQ=255)                                            tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                        >  1:640523/1‑151 (MQ=255)                                             ttGCTGACCAACCCAAACAAAAACACCTTTGGTTGTGTGACGGGTTGTTTTTGGGCGATGTCCGCATACCCCAAGTAAAAGT‑AAAAAGTTTTCCCC‑TG‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                        <  2:474461/150‑1 (MQ=255)                                                 tgcccGACCCTAACTAAAACACCTTTTGCTGTGTGCGGGGCTGTTTTTTGGCGATGCCCGCATAGCCCCAGGAAAAGA‑AATAGGGTTTCCCC‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTa                                                                                                     <  2:1166144/146‑1 (MQ=255)                                                  gcccGAACCTTACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGCATTAACT‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGGTCTTTTGTTGGTCTATAAGTGCGGTAAtt                                                                                                  >  1:436366/3‑151 (MQ=255)                                                    ccGAACCTACCTCAAACACCGTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGCCCGCAAAGCCCGAGAAAAAGT‑AATAGGGTTTCCCC‑TA‑AAAAT‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                    <  2:240924/147‑1 (MQ=255)                                                       aaCCTACCCAAAACACCGTTGGCTGTGTGACGGGCTGTTTTTTGGCGATTCCCGCATACCCCCAGGAAAAGA‑AAAAAGTTTTCCCC‑TA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             <  2:1257743/151‑1 (MQ=255)                                                       aaCCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                             >  1:356321/1‑151 (MQ=255)                                                          cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                          >  2:14037/1‑151 (MQ=255)                                                          cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                          >  2:1055806/1‑151 (MQ=255)                                                           tAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                         >  2:787265/1‑151 (MQ=255)                                                            aaCTAAAACAGCTTTTGTTGTGGGAGGGTCTGTTTTTGGGCGATGCCCACAAACCCCCAGATAAAGT‑AATAGGTTTTCCCT‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAttc                                                                                         <  1:395603/150‑3 (MQ=255)                                                                cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                    >  1:210121/1‑151 (MQ=255)                                                                 aaaaatcGTTTGTTGTGTGATGGGTTTTTTTTTGGTGATTCGCGCAAAGCCCGAGAAAAAGT‑AGT‑AGGTTTCCCC‑TGTAAAAA‑TT‑GCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                   <  1:138386/145‑1 (MQ=255)                                                                 aaaCACCGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGCCCGCATACCCCATGTTAAATT‑AAAAAGTTTTACCT‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAata                                                                                   <  2:44765/151‑1 (MQ=255)                                                                 aaaCAACTTTTGTTGTGTGGGGGGCTTTTTTTGGGCAATGCCCGCATACCCCCAGGAAAAGT‑AATAGGTTTTCCCC‑TG‑AAAAT‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                      <  1:527124/148‑1 (MQ=255)                                                                  aacaacGTTTGCTGTGTGACGGGCTGTTTTTTGGCGATGTCCGCATAGCCCCAGAAAAAGT‑AATAGGTTTTCCCC‑TG‑AAAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                  <  2:597233/151‑1 (MQ=255)                                                                  aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                  >  2:252994/1‑151 (MQ=255)                                                                     aacTTTTGTTGTGGGACGGGTTTTTTTTGGGCGATGCCCGAATAGCCCCAGTAAAAGA‑AATAGGTTTTCCCC‑TA‑AAAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                 <  2:849015/149‑1 (MQ=255)                                                                     aacTTTTGTTGGGTGTGGGGCTTTTTTTTGGCTTTGTCCACATAGCCCGAGTTAAATT‑AATAGGGTTTCCCC‑TA‑AAAAA‑TT‑CCCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGt                                                                               <  1:974947/151‑1 (MQ=255)                                                                     aacTTTGGCTGTGTGGCGGGCTTTTTTTTGGCGAGTCCCGCATACCCCGAGGAAAATA‑AAAAAGTTTTACCT‑TA‑AAAAT‑TT‑CCGGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATc                                                                                 <  2:522827/149‑1 (MQ=255)                                                                      acGTTTGGGGTTTTAGGGGCTGTTTTTTGGCGATGTCCCCATCCCCCCGGGAAAAGT‑AATAATGTTTCCCC‑TA‑ATAAA‑TT‑CCCCCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                <  2:425103/149‑1 (MQ=255)                                                                          ttgttgTGTGACGGTCTGCTTTTTGGCTTTGCCCCCACACCCCCAGTAAAATT‑AATAAGTTTTCCCC‑TA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                           <  1:410539/150‑1 (MQ=255)                                                                            gttgtgtgCGGGGCTTTTTTTGGGCGATGTCCCCATACCCCCAGGAAAAG‑‑AAAAAGGTTTTCCCTTA‑AAAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        <  2:485898/149‑1 (MQ=255)                                                                            gCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑GCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                         >  1:164936/1‑150 (MQ=255)                                                                             gtgtgtgACGGGCTGTTTTTTGCAAATTCCCGCCT‑GCCCCGAGTAAAATAAATAGGTTTTCCCC‑TA‑AAATT‑TT‑CCCGCTACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        <  1:694335/149‑1 (MQ=255)                                                                              tgtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                      >  1:443459/1‑151 (MQ=255)                                                                              tgtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                      >  2:857973/1‑151 (MQ=255)                                                                               gtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     >  2:345419/1‑151 (MQ=255)                                                                                 gtgAGGGGTTGTTTTTGGGCGATGCCCGCATACCCCCGGGAAAAGT‑AAT‑AGGTTTTCCC‑TT‑AAAAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                        <  2:887802/146‑1 (MQ=255)                                                                                 gggcggggTTGTTTTTGGGCGATGCCCGCATCCCCCAAGTAAAAGA‑AATAAGTTTTCCCC‑TG‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                     <  1:18240/144‑1 (MQ=255)                                                                                      gggCTTTTTTTTGGCGATGTCCGCATACCCCCAGAAAAAG‑‑AAAAAAGTTTTCCC‑TT‑AAAAA‑TTCCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                  <  1:497771/147‑1 (MQ=255)                                                                                       ggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                               >  1:757035/1‑149 (MQ=255)                                                                                          tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                           >  2:730961/1‑150 (MQ=255)                                                                                             ttttgGGCGATGCCCGCATACCCCCAGGTAAAGT‑AATAGGTTTCCCCC‑TG‑ATAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                        <  2:1156419/150‑1 (MQ=255)                                                                                                  ggCGATGCCCGCAAACCCCGAGGTAAAGT‑AAGATGTTTCCCCC‑TA‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                    <  2:152725/149‑1 (MQ=255)                                                                                                     gATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  2:878774/1‑150 (MQ=255)                                                                                                     gATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                >  2:945773/1‑150 (MQ=255)                                                                                                       tgcccGCATGGCCCAAGGAAAAGT‑AATAAGTTTTCCCC‑TA‑ATAAA‑TT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                <  2:336604/145‑1 (MQ=255)                                                                                                          ccGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                          >  2:378639/1‑151 (MQ=255)                                                                                                          ccGCATAGCCCCAGGAAAAGA‑AATAAGTTTTCCCC‑TG‑AAAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAag                                              <  2:543987/147‑1 (MQ=255)                                                                                                          ccGCATACCCCGAGGAAAAGT‑AATAAGTTTTACCC‑TA‑ATAAA‑TC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  1:769949/148‑1 (MQ=255)                                                                                                            gcatccccccgttaaaaGT‑AAGAAGTTTCCCCA‑TG‑ATAAT‑TC‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                             <  2:614326/132‑1 (MQ=255)                                                                                                             cATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                        >  1:471216/1‑150 (MQ=255)                                                                                                              aTAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                       >  1:199411/1‑150 (MQ=255)                                                                                                                acccccAGTAAAAGT‑AATAAGGTTTCCCC‑TA‑AAAAA‑TC‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                     <  1:966300/148‑1 (MQ=255)                                                                                                                  ccccAGGAAAAGT‑AATAAGGTTTACCC‑TA‑AAAAA‑TG‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                  <  2:936412/151‑1 (MQ=255)                                                                                                                      cggTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               >  2:559840/2‑150 (MQ=255)                                                                                                                      aGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                               >  1:1016311/1‑150 (MQ=255)                                                                                                                        gTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                            >  2:811088/1‑151 (MQ=255)                                                                                                                         tAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAg                             >  1:1150673/1‑149 (MQ=255)                                                                                                                           aaGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                         >  2:693406/1‑151 (MQ=255)                                                                                                                              t‑aataaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                         >  2:1195993/1‑148 (MQ=255)                                                                                                                              t‑aataaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        >  1:680682/1‑149 (MQ=255)                                                                                                                              t‑aataaGGCTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGAGCGACGCACAATATCGTTGATTTGTTGAGGGCCAGCAAATGTGGAGCAAGAGGCTGTAATTAATAAGgc                           >  1:384623/1‑146 (MQ=255)                                                                                                                                aataaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                       >  1:881446/1‑149 (MQ=255)                                                                                                                                aataaGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                       >  2:571500/1‑149 (MQ=255)                                                                                                                                 ataaGTTTTCCCC‑TG‑ATAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                    <  2:435505/151‑1 (MQ=255)                                                                                                                                 aaaaGGTTTTCCCTTA‑AAAAA‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                        <  1:865223/146‑1 (MQ=255)                                                                                                                                        ttACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt             >  2:410211/1‑151 (MQ=255)                                                                                                                                         tACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAAGCGGTAATTAATAAGGCGCAGAGcgtcgt               >  1:292906/1‑148 (MQ=255)                                                                                                                                          cccc‑TA‑AAAAT‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg           <  2:329307/150‑1 (MQ=255)                                                                                                                                            cccttaaaaaa‑TT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt            <  1:1273080/142‑1 (MQ=255)                                                                                                                                               tG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc        >  2:321560/1‑150 (MQ=255)                                                                                                                                                  aTAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt      >  2:907202/1‑150 (MQ=255)                                                                                                                                                     aa‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg  >  2:1020365/1‑151 (MQ=255)                                                                                                                                                           |                                                                                                                                                   CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGT‑AATAAGGTTTACCC‑TG‑ATAAA‑TT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGG  >  NC_000913/4378653‑4378948

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