Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,568,512	0	C	A	21.3%	128.5 / 14.4	47	R3L (CGT→CTT)	mdtM	multidrug efflux pump/bile salt:H(+) antiporter/Na(+):H(+) antiporter/K(+):H(+) antiporter
Reads supporting (aligned to +/- strand):  ref base C (8/29);  new base A (10/0);  total (18/29)
Fisher's exact test for biased strand distribution p-value = 8.45e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACAC  >  NC_000913/4568379‑4568648                                                                                                                                      |                                                                                                                                          gACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCa                                                                                                                                     <  1:472868/140‑1 (MQ=255) gACATCGGCAATAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCa                                                                                                                                     >  1:888598/1‑140 (MQ=255)   cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTGAAAAAAAGttggcaa                                                                                                                                   >  2:453533/1‑134 (MQ=255)   cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATa                                                                                                                                   <  1:497930/140‑1 (MQ=255)   cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATa                                                                                                                                   >  1:344293/1‑140 (MQ=255)   cATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATa                                                                                                                                   >  2:497930/1‑140 (MQ=255)        gCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACAAGGCATaaaaaa                                                                                                                              >  2:897636/1‑140 (MQ=255)            tAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAaacaac                                                                                                                          >  1:786655/1‑140 (MQ=255)                   aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAATGGGATAAAAAAAAAC‑ACCtttt                                                                                                                  >  2:943842/1‑138 (MQ=255)                      tACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCg                                                                                                               <  2:501273/140‑1 (MQ=255)                      tACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCg                                                                                                               >  1:73868/1‑140 (MQ=255)                       aCCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCgg                                                                                                              <  2:373288/140‑1 (MQ=255)                        ccACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAATGGCATAAAAAAAAAA‑TCCTTGTCCggg                                                                                                             >  2:352476/1‑139 (MQ=255)                           caTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAAC‑TCCTTGTCCGGAGaa                                                                                                          >  1:779264/1‑140 (MQ=255)                                aTGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAAA‑TCCTTTTCCGGAGAAccccc                                                                                                     >  2:443280/1‑139 (MQ=255)                                  gATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAAAAAACAAC‑TCCTTGTCCGGAGAACCCCAAc                                                                                                   >  1:154640/1‑140 (MQ=255)                                   aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGGATAACAAAAAAC‑ACCTTGTCCGGGGAACCCCAACg                                                                                                  >  2:785071/1‑140 (MQ=255)                                     cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAAC‑TCCTTGTCCAGAGAAACCCAACggg                                                                                                >  1:69598/1‑140 (MQ=255)                                     cccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACAGGGCACAACAAAAAAC‑TCCTTGTTCGGAGAAACCCAACggg                                                                                                >  2:870482/1‑140 (MQ=255)                                      ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACgggg                                                                                               >  1:738008/1‑140 (MQ=255)                                      ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAAAAAAAAAC‑TCCTTTTCCGGAGAACCCCAACgggg                                                                                               >  2:100116/1‑140 (MQ=255)                                      ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAACAAACAAC‑ACCTTGTTCGGAGAACCCCAACgggg                                                                                               >  1:269993/1‑140 (MQ=255)                                      ccAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGATTAAAAACAAC‑CCCCTTTCCGGAGAACCCCCACgggg                                                                                               >  2:175807/1‑140 (MQ=255)                                        aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAC‑TCCTTGTCCGGAGAACCCCAACggggga                                                                                             >  2:560683/1‑138 (MQ=255)                                          gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGGGGCATAACAAAAAAC‑ACCCTTTCCGGGGAAACCCAACGGGGaaaa                                                                                           >  2:682062/1‑140 (MQ=255)                                          gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGGATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAAAGGGGgaaa                                                                                           >  1:420943/1‑140 (MQ=255)                                          gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATAAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGaaaa                                                                                           >  1:735129/1‑140 (MQ=255)                                            tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                         <  1:943842/140‑1 (MQ=255)                                            tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                         <  2:94847/140‑1 (MQ=255)                                            tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                         <  1:904978/140‑1 (MQ=255)                                            tGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                         <  1:870482/140‑1 (MQ=255)                                              gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGGAAAAAACGTGGCATAAAAAAAAAT‑TCCTCTGCCGGGGAAACCCAACGGGGAAAACGaa                                                                                       >  2:753905/1‑140 (MQ=255)                                               aTCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGAGGGATAACAAAAAACACCCTTGT‑CGGAGAAACCCCACGGGGgaaaaaaaa                                                                                      >  2:24132/1‑135 (MQ=255)                                                   gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAAA‑CCCCTTTCCGGGGAACCCCCACGGGGGAAAAAAAAATCa                                                                                  >  2:898550/1‑140 (MQ=255)                                                    aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCaa                                                                                 >  1:550475/1‑140 (MQ=255)                                                     tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAAt                                                                                <  1:443280/140‑1 (MQ=255)                                                              aGATATGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAAAAAACAAA‑TCCTTGTCCGGAGAACCCCCACGGGGGAAAAGAAAATCAATACCGTCggg                                                                       >  2:68274/1‑139 (MQ=255)                                                              aGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAACAAAAAAA‑ACCTTTTCCCGGGAACCCCCACGGGGGAAAAGAAAATCAAAACCGTCggg                                                                       >  2:788458/1‑139 (MQ=255)                                                              aGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCTTAAAAAACAAC‑TCCTTGTCCGGAGAAACCCAACGGGGGAAAAGAAATTCAAAACCGTCggg                                                                       >  2:387648/1‑139 (MQ=255)                                                                     cagcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACg                                                                <  2:344293/140‑1 (MQ=255)                                                                       gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACATGGCATAAAAAAAAAC‑TCCTTGTCCGGAGAACCCCCAAGGGGAAAAAAAAAATTAAAAACCTCGGGAATGACGaa                                                              >  2:405069/1‑140 (MQ=255)                                                                        caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa                                                             <  2:179990/140‑1 (MQ=255)                                                                        caAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGaaa                                                             <  1:175807/140‑1 (MQ=255)                                                                                  cAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAAAAAC‑TCCTTGTCCGGAGAAACCCAACGGGGAAAAAGAAAATCAATACCGTCGGGAATGACGAAAGGCAGGGgag                                                   >  2:742588/1‑140 (MQ=255)                                                                                   aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                                  <  1:352476/140‑1 (MQ=255)                                                                                   aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTgaga                                                  <  1:897636/140‑1 (MQ=255)                                                                                       tCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAAA‑TCCTTGTCCGGAGAACCCCCACGGGGGAAAAGAAAATCAAAACCGTCGGGAATGACGAAAGGCAGGGGAGAAggg                                              >  2:204611/1‑140 (MQ=255)                                                                                           cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAAc                                          <  2:786655/140‑1 (MQ=255)                                                                                             ccATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACga                                        <  1:682062/140‑1 (MQ=255)                                                                                                tCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                                     <  2:269993/140‑1 (MQ=255)                                                                                                tCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGAc                                     <  1:100116/140‑1 (MQ=255)                                                                                                     aaaaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaa                                <  1:24132/140‑1 (MQ=255)                                                                                                       aaaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAAt                              <  2:328282/140‑1 (MQ=255)                                                                                                         aCAGCGTGGCGGCATGGCGGGTAAAAAAACGTAGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGt                            <  1:405069/140‑1 (MQ=255)                                                                                                          cAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAAC‑TCCTTTTCCGGGGAACCCCCACGGGGGAAAAGAAATTTAATACCGTCGGGAATGACGAAAGGGAGGGGGGAAGGGCAACGAGACTGCAAATGTa                           >  2:142942/1‑140 (MQ=255)                                                                                                            gCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAC‑TCCTTGTCCGGGGAACCCCAACGGGGGAAAAGAAATTCAATACCGTCGGGAATGACGAAAAGCAGGGGGGAAGGGGAACGAGAATGCAAATGTACt                         >  2:485110/1‑140 (MQ=255)                                                                                                                   ggcATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTc                  <  1:387648/140‑1 (MQ=255)                                                                                                                    gcATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCa                 <  1:753905/140‑1 (MQ=255)                                                                                                                    gcATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCa                 <  2:550475/140‑1 (MQ=255)                                                                                                                        ggCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa             <  1:742588/140‑1 (MQ=255)                                                                                                                        ggCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa             <  1:142942/140‑1 (MQ=255)                                                                                                                        ggCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCAGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTa             <  2:69598/140‑1 (MQ=255)                                                                                                                                aaaaaaCGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACca     <  2:73868/140‑1 (MQ=255)                                                                                                                                   aaaCGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCacac  <  2:748707/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                          GACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAAC‑TCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACAC  >  NC_000913/4568379‑4568648

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