Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	47.5%	96.7 / 13.2	60	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/4);  new base C (0/28);  total (27/33)
Fisher's exact test for biased strand distribution p-value = 8.93e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTG  >  NC_000913/2984292‑2984582                                                                                                                                                       |                                                                                                                                              gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCacac                                                                                                                                              >  1:38798/1‑151 (MQ=255)  aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACAcc                                                                                                                                             <  1:464136/151‑1 (MQ=255)   atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACAccc                                                                                                                                            <  1:424029/151‑1 (MQ=255)       gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGTCACACCCTgg                                                                                                                                         <  1:236090/150‑1 (MQ=255)              gaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCccgcc                                                                                                                                 >  2:200682/1‑151 (MQ=255)                          gAGGGAGCCAAAAAAATGTATGTCGCCTTGCCAGCCATCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGT‑AAATGCCCCCCCCTGGCCCCCGCCGATACAAAGc                                                                                                                       <  1:12965/149‑1 (MQ=255)                            gggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTg                                                                                                                   >  1:493924/1‑151 (MQ=255)                                    aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                           >  2:716077/1‑151 (MQ=255)                                      aaaaTGTATGTCGATTTGCCAGCCATCATTATCCAACAGATAAAGGACGGCAAATATTTATCGGTAGGGAGGAAATGCAGAGTTTTTCGTCACGTTCATTGTT‑CATTCCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                         <  1:395437/151‑1 (MQ=255)                                      aaaaTGTATGTCGACTTGCCAGCCATAATTATTCAACAAATTAAGGCCGGCAGTTATTAATGGGTAAGGAGAAAATGCAGAGCTATTCGTCACTTTAAATGGT‑CAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                         <  1:711287/151‑1 (MQ=255)                                       aaaTGTATGTCGTCTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAAATATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCCTCACGTTAAAGGGT‑AAATGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                        <  1:512008/151‑1 (MQ=255)                                       aaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGT‑CAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                        <  2:539064/151‑1 (MQ=255)                                                 cGACTTGCCAGCCATCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCCCGTTCAATGGT‑CAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                               <  2:623379/150‑1 (MQ=255)                                                 cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                              >  2:340620/1‑151 (MQ=255)                                                  gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                             >  2:87403/1‑151 (MQ=255)                                                   aCTTGCCACCACTCTTTATAAAACAATTAAAGGATGGCAGTTTTTTATGGGAAAGGAAGAAATGCAGAGCTTTTCGTCGCTTTCAGGGGT‑CAATCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                                  <  1:817449/145‑1 (MQ=255)                                                      tGCCAGCCATCATTATTCAAAAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTAAATGGT‑CAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                           <  2:405066/149‑1 (MQ=255)                                                         cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                      >  2:817492/1‑151 (MQ=255)                                                                catTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                               >  2:436394/1‑151 (MQ=255)                                                                  ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                              >  2:571910/1‑150 (MQ=255)                                                                   tattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                             >  1:52517/1‑150 (MQ=255)                                                                        aaCAGATTAATGACGGCAGATTTTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACTTTAAATGGT‑CAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                          <  2:591551/148‑1 (MQ=255)                                                                        aaCAAATTAATGGCGGCAGATATTTCTCGGAAAGGAGGAAATGAAGAGTTTTTCGTC‑CGTTTAAAGGTCCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                       <  2:755539/151‑1 (MQ=255)                                                                                aaTGACGGCAGATTTTTATCGGAAAGGAGGAAATGCAGAGTTATTCGTCACGTTCAATGGT‑AAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                    <  2:868379/146‑1 (MQ=255)                                                                                aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAATGGT‑CAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                <  1:708457/150‑1 (MQ=255)                                                                                aaTGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                               <  2:151992/151‑1 (MQ=255)                                                                                     cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                          >  1:290275/1‑151 (MQ=255)                                                                                      ggaaGATTTTTATCGGTAAGGAGAAAATGCAGGGCTATTCGTCCCTTTAAATGGT‑CAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                          <  2:314462/147‑1 (MQ=255)                                                                                      ggCAGATATTTATGGGAAAGGAGAAAATGCAGAGCTATTCGTCACTTTCAATGGA‑CAAGGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                             <  1:282754/147‑1 (MQ=255)                                                                                      ggCAGATATTTATCGGTAAGGAGGAAATGCAGAGTTATTCGTCACGTTAAATGGT‑CAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                         <  2:75996/151‑1 (MQ=255)                                                                                         atattttttTCGGAAAGGAAAAAATGAAGCGTTATTTGCCGCTTACAAGGAT‑CAATGCCCCGCCCTTTCCCCCGGCGATACAAAGCGTTGCCAGTCCTTTGCTGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                              <  1:378012/138‑1 (MQ=255)                                                                                          gATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCa                                                      >  1:194833/1‑150 (MQ=255)                                                                                              tttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                 >  1:889122/1‑151 (MQ=255)                                                                                              tttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                 >  1:140771/1‑151 (MQ=255)                                                                                              tttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                 >  1:105700/1‑151 (MQ=255)                                                                                               ttATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                 >  2:742705/1‑150 (MQ=255)                                                                                                  tCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                             >  2:896460/1‑151 (MQ=255)                                                                                                    ggaaaGGAGAAAATGCAGAGTTATTCGCCCCGTAAAAGGGT‑AAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                            <  1:644509/147‑1 (MQ=255)                                                                                                        aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                       >  2:391629/1‑151 (MQ=255)                                                                                                         ggagaaAATGCAGAGCTATTCGTCACTTTCAAGGGT‑CAAGGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                      <  1:266243/151‑1 (MQ=255)                                                                                                         ggagaaAATGCAGAGCTATTCGTCACGTTAAATGGT‑CAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                         <  1:911504/148‑1 (MQ=255)                                                                                                           aggaAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                    >  2:103627/1‑151 (MQ=255)                                                                                                           aggaAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                    >  2:306928/1‑151 (MQ=255)                                                                                                                   cAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                             >  1:70400/1‑150 (MQ=255)                                                                                                                    agagCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                             >  1:162636/1‑149 (MQ=255)                                                                                                                        cTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                       >  1:817699/1‑151 (MQ=255)                                                                                                                         tATTCGTCACGTTAAATGGC‑CAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                        <  1:932340/149‑1 (MQ=255)                                                                                                                         tATTCGCCACGTTCAATGGT‑AAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGt                          <  2:144576/147‑1 (MQ=255)                                                                                                                          aTTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                     >  2:605678/1‑151 (MQ=255)                                                                                                                           ttCGTCCCGTTCAATGGT‑AAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                      <  2:665527/149‑1 (MQ=255)                                                                                                                           ttCGTCACTTTAAATGGT‑AAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                      <  2:930635/149‑1 (MQ=255)                                                                                                                           ttCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGGGATc                    >  2:478692/1‑151 (MQ=255)                                                                                                                             cGTCACGTTAAATGGT‑CAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                   <  1:404467/150‑1 (MQ=255)                                                                                                                              gTCACGTTCAATGGT‑CAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                  <  1:493368/150‑1 (MQ=255)                                                                                                                                       aagggT‑CAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat          <  2:44644/146‑1 (MQ=255)                                                                                                                                       aaTGGT‑CAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTAACCGAGTGCGATCGCGCCACcatt         <  2:712288/150‑1 (MQ=255)                                                                                                                                       aaTGGT‑AAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt         <  2:945558/150‑1 (MQ=255)                                                                                                                                         tGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGAGCGATCGCGCCACcattgac      >  1:377117/1‑151 (MQ=255)                                                                                                                                          ggT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac      >  2:580693/1‑150 (MQ=255)                                                                                                                                               aaTGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg  <  1:72603/150‑1 (MQ=255)                                                                                                                                               aaTGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg  >  2:562200/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                              GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGT‑CAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTG  >  NC_000913/2984292‑2984582

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