Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,332,791	0	A	C	46.0%	104.8 / 17.8	63	V161G (GTG→GGG)	basS	sensory histidine kinase in two‑component regulatory system with BasR
Reads supporting (aligned to +/- strand):  ref base A (32/2);  new base C (0/29);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.12e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332642‑4332941                                                                                                                                                        |                                                                                                                                                       cATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCAc                                                                                                                                                       >  1:210844/1‑151 (MQ=255) cATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCAc                                                                                                                                                       >  1:507886/1‑151 (MQ=255)      tGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCg                                                                                                                                                  >  2:136058/1‑151 (MQ=255)            ttACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTg                                                                                                                                            >  1:50824/1‑151 (MQ=255)             tACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTgg                                                                                                                                           >  2:360262/1‑151 (MQ=255)               ccGGAAAAAAATGCCTTCCCGCCCCGCCCCGGTTGCGGCGGCGGGGAGCCCCTCCCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑AAAA‑CGCCCCCCCGCCAGTg                                                                                                                                            <  2:473851/148‑1 (MQ=255)               ccGGAAAAAAATGACTGCCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCTTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCCCCCCCGCCAGTg                                                                                                                                            <  1:404974/148‑1 (MQ=255)                     gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                    >  2:97666/1‑150 (MQ=255)                     gaaAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGc                                                                                                                                    >  1:520761/1‑150 (MQ=255)                            cccgTCCGGCCCGCCCGCGTTGCAGCCGTTGGAACCCTCCCTCCATCATCTAAAAACGCCGCCCCCCTAAGGGCGCTCCAAAAATGGGATGCGTTTTCGCCAGCTGTCCCAGATG‑CAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTg                                                                                                                            <  2:195811/148‑1 (MQ=255)                             cTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTgc                                                                                                                           >  2:378799/1‑151 (MQ=255)                                tCCGGCACGCGCCAGTTGCAGGAGCTGGGAGACTCTCTCCATCATCTGATCAAGCCCTGCCACTAACGGTGCTTCATCAATGTGATGCGTTTTCGCCAACAGTTCAAGATG‑CAAA‑CGCACCCCCGCCAGTGGCATTCGCAGTTCGTgc                                                                                                                           >  2:746423/1‑148 (MQ=255)                                 ccGCCACGCGCCAGTTGCAGCGGCGGGGAGACCCTCTCAATATTCTGACCAAGCCGTCCCACAAGCGGTCCTACATCAATGTGATGCTTTTTCGCCAGCAGTCCCAGATG‑AAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCg                                                                                                                              <  1:234525/144‑1 (MQ=255)                                        gcgcGAGTGGCGGCAGCGGGGAGCCGCTCCCCATCTTCTGATCAGCCGGTCCCTCTAAGGGTCCAACATCAATGTGAGGCTTTTCCGCCAGCAGTCCGAGA‑G‑CAAAGCGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCa                                                                                                                  <  1:306821/149‑1 (MQ=255)                                        gcgcCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGc                                                                                                                >  1:137885/1‑151 (MQ=255)                                         cgcCAGTGGCAGCAGCGGGGAGCCGCTCCCCATCATCTGATCAAGCCGTGCCACTACCGGTGCTACATAAATGTGATGCTTTTTCGCCAGCATTTCCAGATG‑CAAC‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAg                                                                                                                 <  1:560027/149‑1 (MQ=255)                                           ccAGTTGCAGCGGTGGGGAGCCCCTCTCCACCATTTGATCAAGCCGTGCCATTAAGGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTCCCAGAGG‑CAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCa                                                                                                                  <  1:125549/146‑1 (MQ=255)                                            cAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGt                                                                                                             >  2:612752/1‑150 (MQ=255)                                              gTTCCACCAGGTGGGAGCCGCTCCCCATCATCTAACCAAGCCGCCCCACTAACGGTGCAACATCAATGTGTTGCGTTTCCGCCAGCAGTTCCAGATG‑CAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaaa                                                                                                          <  2:574251/151‑1 (MQ=255)                                               tcgcagcgggggggAGCCCCTCCCCATCTTCTGACCACGCCGTGCCACAACCGGTGCTACATCAATGTGATGCTTTTCCCCCAGCATTTCCAGATG‑CAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTaa                                                                                                           <  2:205662/138‑1 (MQ=255)                                                tGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaac                                                                                                         >  1:787062/1‑150 (MQ=255)                                                      ggcTGGGAGCCCCTCTCCTTCTTCTGACCAAGCGGGCCCATAACCGGTGCTACAACAGTGTGATGCGTTTCCGCCAGCTGTCCGAGTGG‑AAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                                  <  2:671830/150‑1 (MQ=255)                                                       gcggggAGACCCTCCCCATCATCTGATAAAGCGGTCCCAATACCGGTGCTACCTCAATGGGGTGGTTGTTCGCCAGCGGTTCGAGAAC‑GAAT‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                        <  2:797478/141‑1 (MQ=255)                                                         ggggggccGCTCCCCATCTTCGAACCAAGCCGTCCCACAACCGGTGCTACATAAATGTGATGCTTTTTCGCCAGCAGTCCCAGATG‑AAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTc                                                                                                 <  2:402143/142‑1 (MQ=255)                                                          gggggcccccctcCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACACAAATGTGATGCTTTTTCGCCAGCAGTTCCAGATG‑AAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAaaca                                                                                                        <  1:271850/131‑1 (MQ=255)                                                          ggggcgccccccccaccaTCTAAACACGCCGTCCCACTACCGGTGCTACAAAAAGGGGATGCTTTTCCGCCAACAGTTCCGAATG‑CAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCttt                                                                                                  <  2:899878/135‑1 (MQ=255)                                                            gagcccctctCCATCATCTGATCAAGCCGTCCCACTAACGGTGCTACATCAATGTGATGCTTTTCCGCCAGCATTTCCAGATG‑CAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTAt                                                                                            <  1:45315/145‑1 (MQ=255)                                                                 gCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGc                                                                                       >  1:297802/1‑151 (MQ=255)                                                                  ctctCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCg                                                                                      >  2:467253/1‑151 (MQ=255)                                                                       catcatCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTg                                                                                 >  2:712863/1‑151 (MQ=255)                                                                        atcatcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTgg                                                                                >  1:554319/1‑151 (MQ=255)                                                                        atcatcTGATCAAGCCGTACCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTgg                                                                                >  2:570197/1‑151 (MQ=255)                                                                            tcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                              >  2:387255/1‑149 (MQ=255)                                                                            tcTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTc                                                                              >  2:494465/1‑149 (MQ=255)                                                                               gATCAAGCCGTGCCACTAGGGGTGCTACATCAATGTGATGCTTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                          <  1:623359/150‑1 (MQ=255)                                                                                accacgccgccccccTAGCGGTGCTACATCAATGTGAGGCTTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAgg                                                                           <  2:552644/134‑1 (MQ=255)                                                                                accacgCCGTCCCACTAAGGGTTCTACATCAATGTGTTGTTTTTCCGCCAGCAGTTCCAGATG‑AAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGc                                                                          <  2:697325/144‑1 (MQ=255)                                                                                   acgCCGGGCCATAACCGGTGCTACATCAAGGTGTTCCGTTTTCGCCAGCGGTACCAGATC‑CAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACt                                                                      <  1:642687/148‑1 (MQ=255)                                                                                      ccGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                  >  1:448592/1‑151 (MQ=255)                                                                                            cACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGAcc                                                                  >  1:780808/1‑145 (MQ=255)                                                                                             acgagaggTGCTACATCAAGGTGAGGCGTTTCCGCCAGCCGTCCAAGAGG‑AAAA‑CACACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATc                                                             <  2:662981/143‑1 (MQ=255)                                                                                                 aCGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTa                                                        >  1:333912/1‑150 (MQ=255)                                                                                                   ggTGCACCATCAATGTGGTG‑GTTTTGCCCCGGCGGTCCAGAAGCCAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGttt                                                         <  2:413031/147‑1 (MQ=255)                                                                                                      gCTAAATAAAGGGGATGGTTTTTCCCCAGCATTCCCAGGTG‑CAAC‑GGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGtt                                                          <  1:603700/143‑1 (MQ=255)                                                                                                        tACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                                 >  1:241408/1‑150 (MQ=255)                                                                                                           aTCAATGTGATGGTTTTCCCCCAGCAGTCCCAGACG‑CAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGa                                                 <  1:741131/147‑1 (MQ=255)                                                                                                           aTCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAacc                                             >  2:689974/1‑151 (MQ=255)                                                                                                              aattggATGCTTTTTCCCCAGGCGTTCCAAAAG‑AAAA‑CCCCCCCCCGCCCGTTGCGTTCGCAGTTCGGTGGGGACCTCAGGGGTTAACCACCTTTCCTTTTCCAGCGTTCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGaaa                                               <  1:844787/141‑1 (MQ=255)                                                                                                                 gtgATGTTTTTCCGCCAGCAGTCCCGGAGG‑AAAC‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGc                                        <  1:219245/150‑1 (MQ=255)                                                                                                                         ttttCGCCACCAGTTCCAGAGG‑AAAA‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAttc                                <  1:849272/150‑2 (MQ=255)                                                                                                                         tttcCGCCAGCATTCCCAGATG‑CAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTc                                     <  2:101416/145‑1 (MQ=255)                                                                                                                          tttCGCCAGCAGTTCCAGATG‑AAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGAtt                                 <  1:891286/148‑1 (MQ=255)                                                                                                                          ttccGCCAGCATTTCAAGATG‑CAAA‑CCCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                               <  2:447167/147‑1 (MQ=255)                                                                                                                           ttCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCg                              >  1:27927/1‑150 (MQ=255)                                                                                                                             cGCCACCTGTCCAAGATC‑CAAC‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTc                               <  1:57312/147‑1 (MQ=255)                                                                                                                              gCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAggg                          >  1:379348/1‑151 (MQ=255)                                                                                                                               ccAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGt                         >  1:93541/1‑151 (MQ=255)                                                                                                                                 agcagTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTgg                       >  2:786288/1‑151 (MQ=255)                                                                                                                                    agTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                    >  1:727603/1‑151 (MQ=255)                                                                                                                                    agTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGgcgc                    >  1:56606/1‑151 (MQ=255)                                                                                                                                             tG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATgg           >  1:467537/1‑151 (MQ=255)                                                                                                                                              tgaaaa‑CGCCCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGa               <  1:56248/144‑1 (MQ=255)                                                                                                                                                      gCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGCGGGTGGCGCTGTGAATGGCAATggg    >  1:391896/1‑151 (MQ=255)                                                                                                                                                        aCCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCg  >  2:715416/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       CATGTTGATAATTACCGGAAGAAAATGACTGTCCGGCACGCGCCAGTTGCAGCAGCTGGGAGACGCTCTCCATCATCTGATCAAGCCGTGCCACTAACGGTGCTACATCAATGTGATGCGTTTTCGCCAGCAGTTCCAGATG‑CAAA‑CGCACCCCCGCCAGTGGCGTTCGCAGTTCGTGCGCGACGTCAGCGGTAAACAACCTTTCGTTATCCAGCGTGCTGGTCAGGCGACTGACCAGATCGTTTAACGCCGAAACCACCGCTTCGATTTCGAGGGTGGCGCTGTGAATGGCAATGGGCG  >  NC_000913/4332642‑4332941

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