Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,568,512	0	C	A	39.3%	98.4 / 19.7	61	R3L (CGT→CTT)	mdtM	multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (9/28);  new base A (24/0);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 7.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAAC  >  NC_000913/4568366‑4568659                                                                                                                                                   |                                                                                                                                                    cagGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg                                                                                                                                                 <  1:613452/151‑1 (MQ=255)     ggCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATa                                                                                                                                             <  2:971971/151‑1 (MQ=255)       ccAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaa                                                                                                                                            >  1:425545/1‑150 (MQ=255)       ccAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGAAAAAAAAAGAGGcgaaa                                                                                                                                            >  2:809573/1‑146 (MQ=255)        cAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaaca                                                                                                                                          <  2:517599/151‑1 (MQ=255)        cAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGAGGCCTaa                                                                                                                                            >  2:549127/1‑149 (MQ=255)           aCTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGTCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACACCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaacc                                                                                                                                       <  2:392839/151‑3 (MQ=255)              gACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAaacaac                                                                                                                                    <  1:748263/151‑1 (MQ=255)                          aaaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCgg                                                                                                                         <  1:804845/150‑1 (MQ=255)                          aaaaTCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGga                                                                                                                        <  2:334415/151‑1 (MQ=255)                                aCGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGAAAAAAAAAGTGGGATAACAAAAAACTCCTTGTCCGGAGAAccc                                                                                                                  >  2:645805/1‑151 (MQ=255)                                       acaTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAACTCCTTGTCCGGAGAAACCCAAAggg                                                                                                           >  2:909541/1‑151 (MQ=255)                                         aTTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACAAGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAAAgggg                                                                                                          >  2:460968/1‑150 (MQ=255)                                             aTGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAAc                                                                                                     <  2:387941/151‑1 (MQ=255)                                               gATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACg                                                                                                    <  1:810126/150‑1 (MQ=255)                                               gATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                                                                                   <  2:318712/151‑1 (MQ=255)                                                aTCCCAGGCTGGCTCAGACCCGCCGCCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGCCGACAAAAACTGGCATAAAAAAAAACCCCCTTGCCGGAGCACCCCCACCGGGAAAAcca                                                                                                   >  2:798768/1‑148 (MQ=255)                                                aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGaa                                                                                                  <  2:754428/151‑1 (MQ=255)                                                aTCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAAAGa                                                                                                   >  1:655280/1‑150 (MQ=255)                                                    cAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACATGGGCTAACAAAAAACACCTTGTTCGGAGAACCCCAACGGGGAAAAAGaaa                                                                                                 >  2:699432/1‑148 (MQ=255)                                                    cAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGCGGCCTAAAAAACAACCCCTTGTCCGGAGAACCCCAAAGGGGAAAACGaaa                                                                                                 >  2:150443/1‑148 (MQ=255)                                                     aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAAACAACAACCCCTTGTCCGGAGAACCCCAACGGGGAAAAAGAAAtt                                                                                               >  2:1006524/1‑149 (MQ=255)                                                     aGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGaaaat                                                                                               >  1:606510/1‑147 (MQ=255)                                                       gCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAAt                                                                                           <  1:655840/151‑1 (MQ=255)                                                          ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCCACGGGGGAAACGAAATTCAATa                                                                                          >  1:302194/1‑149 (MQ=255)                                                              cAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCg                                                                                    <  2:1792/151‑1 (MQ=255)                                                               aGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAACAGGCATAACAAAAAACCCCCTGTCCGGAGAACCCCAACGGGGGAAAAGAAAATCAAAACCCTCgg                                                                                   >  2:439605/1‑151 (MQ=255)                                                                gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACCCCCTGTCCGGGGAACCCCAAAGGGGAAAACGAAATTCAATACCGTCg                                                                                    >  2:111494/1‑149 (MQ=255)                                                                  tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCATAAAAAACAACTCCTTGTCCGGAGAACACCAACGGGGGAAAAGAAAATCAAAACCGTCgg                                                                                   >  2:73771/1‑148 (MQ=255)                                                                  tCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGAAAAAAAAAGCGGCATAAAAAAAAAACCCTTTTCCGGGGAAACCCAAAGGGGAAAAAGAAAATCAAAACAGGCggg                                                                                  >  2:168709/1‑148 (MQ=255)                                                                   ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAAAAAACAAATCCTTGTCCCGAGAACCCCAACGGGGGAAAAGAAATTTAATACcggcgg                                                                                   >  2:348280/1‑147 (MQ=255)                                                                       cGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACg                                                                           <  2:512697/151‑1 (MQ=255)                                                                        gACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGGAAAAGAAATTCAATACCGTCGGTAATGACg                                                                           >  1:784355/1‑150 (MQ=255)                                                                          cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGGAATGACGaa                                                                         >  1:96757/1‑150 (MQ=255)                                                                           aGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAg                                                                       <  2:214991/151‑1 (MQ=255)                                                                           aGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAg                                                                       <  1:933005/151‑1 (MQ=255)                                                                           aGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAg                                                                       <  1:914651/151‑1 (MQ=255)                                                                            gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAg                                                                       <  1:850808/150‑1 (MQ=255)                                                                             aTACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc                                                                     <  1:942075/151‑1 (MQ=255)                                                                                cGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCCTGTCCGGAGAACCCCCAAGGGGGAAACGAAATTCAATACCGTCGGGAAAGACGAAAGGCAgg                                                                  >  1:1002404/1‑151 (MQ=255)                                                                                      aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATTACAAAAAACTCCTTGTGCCGGGAAACCCAACAGGGAAAAAGAAAATTAAAACACTCGGTACAGACGAAAAGGAGGGgaga                                                             >  2:965497/1‑150 (MQ=255)                                                                                       aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAg                                                           <  2:333342/151‑1 (MQ=255)                                                                                          tCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGCAAAAAAAACTGGCAGAACAAACAACTCCTTGTCCGGAGAACCCCCACCGGGGAAAAGAAATTAAATAACGTCGGGCATGGCGAAAGGCAGGGGAGaa                                                            >  2:1002240/1‑147 (MQ=255)                                                                                            aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGGGGCATAAAAAAAAACTCCCTTTCCGGGGAACCCCAACGGGGAAAACGAAATTCAAAACCGCCGGGAATGACGAAAAGCCGGGGAAAAGGGCa                                                       >  2:932175/1‑150 (MQ=255)                                                                                            aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAACTCCCTTTCCCGAGGACCCCCAAGGGGGAAACGAAATTTAATACCCCCGGGAATTAAGAAAAGCAGGTGAGAAggg                                                         >  1:626207/1‑148 (MQ=255)                                                                                             tACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGGAAACGAAATTCAATAACGTCGGGAATGAAGAAAGGCAGGGGAGAAGGGCaa                                                      >  1:118899/1‑150 (MQ=255)                                                                                              aCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAAc                                                     <  1:321341/150‑1 (MQ=255)                                                                                                aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACgag                                                  <  1:764153/151‑1 (MQ=255)                                                                                                   aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAACTCCTTGTCCGGAGAACCCCAAAGGGGAAAACCAAATTTAATACCGTCGGTAATGAAGGAAGGGAGGGGGGAAGGGGAACgagaat                                               >  1:686198/1‑149 (MQ=255)                                                                                                         gCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAGCTCCGTGTCCGGAGAACGCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTg                                              <  2:13819/146‑1 (MQ=255)                                                                                                         gCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAAt                                         <  1:188848/151‑1 (MQ=255)                                                                                                         gCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAACTCCTTGTCCGGGGAACCCCAACGGGGGAAACGAAATTCAATACCGTCGGGAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCaaa                                          >  1:194579/1‑150 (MQ=255)                                                                                                           cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTTAAAAAACGTGGCCTAAAAAAAAACTCCTTGTCCGGGGAAACCCAACGGGGGAAAAGAAAATCAATACCGTCGGGAATGACGAAAAGACGGGGGGAAGGGCAACGAGAATGCAAATg                                        >  2:894359/1‑150 (MQ=255)                                                                                                           cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGt                                       <  1:426903/151‑1 (MQ=255)                                                                                                           cATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCAAAACAAACAACTCCTTGTCCGGAGAAACCCAAAGGGGAAAACGAAATTCAATAACGTCGGGAAAGACGAAAGGCAGGGGAGAAGGGCAAAGAGACCTCAAATg                                        >  1:248053/1‑150 (MQ=255)                                                                                                                 gAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGaa                                 <  1:563165/151‑1 (MQ=255)                                                                                                                  aaaaaCAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAACAAATCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTTAAAAACGTCGGGAATGACGAAAGGAAGGGGAGAAGGGCAAAGAGACTGCAAAAGTAATGaa                                 >  2:947862/1‑150 (MQ=255)                                                                                                                   aaaaCAGCGTGGCGGCATGGCGGGTAAAAAAAAGTGGCATAACAAAAAACTCCCTGTCCGGGGAACCCCCACCGGGGAAACGAAAATTCACAACGTCGGGAATGACAAAAAGCAGGGGAGAAGGGAAACCAGACTGCAAATGGAATGAAtt                               >  2:731721/1‑151 (MQ=255)                                                                                                                        aGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATg                          <  1:916590/151‑1 (MQ=255)                                                                                                                               cggcATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCaa                   <  1:348479/151‑1 (MQ=255)                                                                                                                                                 aaCGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAac  >  1:19604/1‑150 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAGGGGCCAGACTGACATCGGCATTAAAATCACGTACCACATTAATGATCCCAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGCAGGTGAGAAGGGCAACGAGACTGCAAATGTACTGAATTTCATGTAATCAACCACACTGCCTGTGAAC  >  NC_000913/4568366‑4568659

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