Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,598,973	0	A	T	41.9%	98.2 / 18.9	62	M229K (ATG→AAG)	livJ	branched‑chain amino acid ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (31/5);  new base T (0/26);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.62e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAACCTGATCGTAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTTTCTTCAGGCCGTCCTGCACC  >  NC_000913/3598824‑3599119                                                                                                                                                      |                                                                                                                                                    gAACCTGATCGTAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCAt                                                                                                                                                    <  1:289669/151‑1 (MQ=255)  aaCCTGATCGTAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCAtt                                                                                                                                                   >  2:1077665/1‑151 (MQ=255)       gATCGTAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCgg                                                                                                                                             >  1:791943/1‑151 (MQ=255)           gTAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGa                                                                                                                                         <  2:537160/151‑1 (MQ=255)            tAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGAt                                                                                                                                        >  2:590089/1‑151 (MQ=255)                          ggTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCgtagta                                                                                                                           >  2:911091/1‑150 (MQ=255)                            tCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGAATTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAAc                                                                                                                        >  2:409955/1‑151 (MQ=255)                             cccccccccccctccctCTGATTCCCCCGCATTTTAAGACAGCAAAACTTTGCCCCCCCCTCCCGCCCCAATAAACGGTTTTTTCCCCCCTCCCGCGCGTCCGTGACGCGGGTTTTGCCCCTTTC‑CGGGGTGATAACCGCCGTAGTaaa                                                                                                                         <  2:1009589/140‑1 (MQ=255)                              aCCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGAGATAACCGCCGTAGTAAACg                                                                                                                       >  2:743721/1‑150 (MQ=255)                                        cTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATAtt                                                                                                            >  2:542184/1‑151 (MQ=255)                                           ccGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATtctc                                                                                                         >  1:141588/1‑151 (MQ=255)                                           ccGCTGATCCGCCCGCAATGTTAGACAGCAAAACTTTACCCACCCCTCCCGCCCCAAAAAACTGATTTTCCACCCCGGCCGCGCGTGCCTGACGAAGGTTTTCCCCCTTTC‑CGGGGTGATAACCGCCGTAGTAAACGAAGTCGATATtctc                                                                                                         <  1:455184/151‑1 (MQ=255)                                            cGCTGATTCCCCCGAAATGTTAGACAGCAAAGCGTTACCCACCCCTCCCGCCCCAATAAACTGATTTTTCACCCCGGCCGCGCGTGCCTGACGCAGGTTTTGCCCCTTTC‑CGGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTctct                                                                                                        <  1:804987/151‑1 (MQ=21)                                              cTGATCCGCCCGCAGTGTTAGACAGCGAAACGTTACCCACCCCTCCCGCCCCAATAACCGGATTTTTCAGCCCGGCCCCGCTGGCCTAACGCAGTATTTGCCCCTTTT‑CCGGGGGATAACCGCCGTAGTAAACGAAGTCGATATTCTCtt                                                                                                       <  2:410933/150‑1 (MQ=21)                                               tGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCtttt                                                                                                     >  1:714740/1‑151 (MQ=255)                                                gttcccccgacattGTTAGACAGCAAACGGTTAGCCCCCCCTCCCGCCCCCATAACCTGATTTTCCACCCCTCCCGCGCGTGCCTGACGAAGGTTTTCCCCCTTTT‑CGGGGTAAAAACCGCCGTAGTAAACGAAGTCGATATTCTCttt                                                                                                      <  2:540092/136‑1 (MQ=255)                                                         gaaaTTTTAGACGGAAAAACTTTAGCCCCACCTTCCCCCCCAATAAACTGATTTTTCAACCCCGCCGCGCGTGCATGGCGAAGGTTTCCCCCCTTTT‑CCGGTTAAAAACCGCCGTAGTAAACGAAGTCGATATTCTCttt                                                                                                      <  1:214759/138‑1 (MQ=16)                                                         gCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGCGATAACCGCCGTAGTAAACGCAGTCGATATTCTCTTTTTTCAGAcgcg                                                                                           >  2:808131/1‑151 (MQ=255)                                                           attGTAAGACAGCAAACCGTTACCCCCCCCTCCCGGCCCAATAACCTGATTTTCCACCCCTGCCGCGCGTGCCTGACGCAGGATTTCCCCCTTTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGAcgc                                                                                            <  1:971145/146‑1 (MQ=21)                                                             tGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGccac                                                                                       >  1:183970/1‑151 (MQ=255)                                                                tAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGccacc                                                                                      >  1:661148/1‑149 (MQ=255)                                                                   aCAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGt                                                                                 >  2:1024668/1‑151 (MQ=255)                                                                      gcaaaccGTGAGCCACCCCTCCCGCCCCAATAAATTGATTTTCCACCCCTGCCGCGCGTGCTTGCCGTAGGTTTTCCCCCTTTT‑CGGGGTAAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAg                                                                                    <  2:905352/139‑1 (MQ=19)                                                                         aaaCGTTCCCCACCCCTCCCGCCCCCATAAACTGATTTTTCAGCCCTGCCGCGCGTGCTTGACGCAGGATTTGCCCCATTT‑CGGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGc                                                                                   <  1:547468/143‑1 (MQ=25)                                                                         aaaCGTTACCCACCCCTCCCGCCCCCAAAAACTGATTTTTCACCCCTGCCGCGCGTGCCTGCCGCAGGTTTTCCCCCTTTT‑CGGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGt                                                                                 <  1:682449/145‑1 (MQ=25)                                                                          accgttcccccccccTCCCGCCCCCAAAACCTGGTTTTCCACCCCCCCCGCGCTGGCGTGACGAAGATTTTGCCCCTTTC‑CGGGGTAAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCg                                                                                  <  2:395605/130‑1 (MQ=12)                                                                          accGTTACCCCCCCCTCCCGCCCCAATAACCTGATTTTCCGCCCCTGCCGCGCGTCCCTGACGCAGGATTGGCCCCATTT‑CGGGGTGAAACCCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCg                                                                                  <  2:371455/141‑1 (MQ=21)                                                                             gTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCt                                                                       >  1:662340/1‑151 (MQ=255)                                                                              taaccccccccTTCCGGCCCCATAACTTGTTTTTTCCGCCCTCCCCCGCGTCCCGGACGAAGGTTTTCCCCCTTTT‑CGGGTTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGaaa                                                                          <  2:169452/138‑1 (MQ=21)                                                                                  ccccccctcccccccccAAAAACTGATTTTTCACCCCCGCCGCGCGTGCCTGACGCAGGATTTCCCCATTTT‑CGGGGTAAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCt                                                                       <  1:952469/133‑1 (MQ=21)                                                                                    acacCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTccc                                                                >  2:269405/1‑151 (MQ=255)                                                                                      acCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCcg                                                              >  1:548069/1‑151 (MQ=255)                                                                                         ttCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCcggc                                                            >  2:681215/1‑150 (MQ=255)                                                                                           ccGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCcggcgg                                                          >  2:850191/1‑150 (MQ=255)                                                                                           ccGCCCCCATAAACTGTTTTTCCACCCCTGCCGCGCGTCCCTGACGCAGGTTTTGCCCCTTTT‑CGGGGTGAAACCCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCcggc                                                            <  2:558404/148‑1 (MQ=25)                                                                                            cgcccccAAAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCGGGATTTCCCCCTTTC‑CGGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTg                                                        <  1:713100/148‑1 (MQ=25)                                                                                             ggCCCCAAAAAGTGATTTTCCACCCCGGCCGCGCGTGCCGGACGGGGGTTTTGCCCCTTTT‑CGGGGTGAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTg                                                        <  1:555023/150‑1 (MQ=21)                                                                                              gCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGAt                                                      >  1:753693/1‑151 (MQ=255)                                                                                                   aaaaccggAGTTTTCACCCCGGCCGCGCGTGCCGGACGCAGGATTTGCCCCTTTC‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCAt                                                 <  1:662828/144‑1 (MQ=25)                                                                                                   aTAACCTGATTTTCCGGCCCTGCCGGGCGTGCCTGACGCGGGTTTTCCCCCTTTT‑CGGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATg                                                     <  2:161458/147‑1 (MQ=25)                                                                                                    tAAACTGAGTTTTCAGCCCTGCCGCGCGCGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGCCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATc                                                >  1:374135/1‑151 (MQ=255)                                                                                                      accTGAGTTTTCACCCCTGCCGCGCGTGCCTCACGAAGGTTTTGCCCCTTTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATc                                                <  2:247429/147‑1 (MQ=37)                                                                                                         tGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCaa                                              >  1:196993/1‑148 (MQ=255)                                                                                                          ggttttcCACCCCCGCCGCGCTTGCCTGACGCAGGTTTTCCCCCTTTT‑CGGGGTAATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCaaa                                             <  1:505294/145‑1 (MQ=21)                                                                                                              tttCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAAcacc                                      >  1:1049472/1‑151 (MQ=255)                                                                                                              ttccaccccGGCCGCGCGGCCCTGCCGCAGATTTTCCCCCATTC‑CGGGGTGAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAacac                                       <  1:328442/144‑1 (MQ=21)                                                                                                               ttCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAAcacc                                      >  2:449379/1‑150 (MQ=255)                                                                                                                   gCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGtt                                 >  2:487072/1‑151 (MQ=255)                                                                                                                    cccGGCCGCGCGTGCCTGACGCGGGTTTTCCCCCTTTC‑GCGGGTGAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGttt                                <  2:855411/151‑1 (MQ=25)                                                                                                                      cTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGc                              >  1:329951/1‑151 (MQ=255)                                                                                                                         ccGCGCGTGCCTGACGCAGGATTTGCCCCTTTC‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCAt                            <  1:42638/150‑1 (MQ=255)                                                                                                                          cgcgcgTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCAtt                           >  1:219631/1‑150 (MQ=255)                                                                                                                            cgcgTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTg                          >  1:518978/1‑149 (MQ=255)                                                                                                                            cgcgTCCCGGACGCAGGTTTTGCCCCTTTC‑CGGGGTAATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCa                             <  2:950127/146‑1 (MQ=25)                                                                                                                               ggccgtaccGCAGGCTTTCCCCCTTTA‑CCGGAAGAGAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCAt                            <  2:746704/136‑1 (MQ=25)                                                                                                                                  ccTGCCGGAGGTTTTGCCCCTTTT‑CCGGGTAAAAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTttct                   <  1:664635/150‑1 (MQ=25)                                                                                                                                   cTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTttcttc                 >  2:25471/1‑151 (MQ=255)                                                                                                                                               tttccccTTTTCCCGGTTAATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTTTCTTCAGGcc            <  2:1075918/142‑1 (MQ=25)                                                                                                                                               ttGCCCCTTTT‑CGGGGTAAAAGCCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTTTCTTCAGGCCGTcc        <  2:707009/148‑1 (MQ=255)                                                                                                                                               ttGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTTTCTTCAGGCCGTCCTGc     >  2:917458/1‑151 (MQ=255)                                                                                                                                                  ccccATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTTTCTTCAGGCCGTCCTGCAcc  >  1:856480/1‑151 (MQ=255)                                                                                                                                                  ccccATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCCCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGTAGCCCTCAAAGAACAGCACGTTTGCATTGCCTTTCTTCAGTCCGTCCTGCAcc  >  2:110724/1‑151 (MQ=37)                                                                                                                                                      |                                                                                                                                                    GAACCTGATCGTAGTTCTTCGGCTTGGTCACCAGCAGCCCTTCCGCTGATTCGCCCGCAATGTTAGACAGCGAAACGTTAGCCACACCTTCCGGCCCCATAAACTGAGTTTTCAGCCCTGCCGCGCGTGCCTGACGCAGGATTTGCCCCATTT‑CCGGGTGATAACCGCCGTAGTAAACGAAGTCGATATTCTCTTTTTTCAGACGCGCCACCAGCGTTGAGAAATCTTTTTCCCCGGCGGTGATGCCATCAAAGAACACCACGTTTGCATTGCCTTTCTTCAGGCCGTCCTGCACC  >  NC_000913/3598824‑3599119

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