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

AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG  >  NC_000913/2984293‑2984587                                                                                                                                                      |                                                                                                                                                   aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACAcc                                                                                                                                                  >  1:945636/1‑151 (MQ=255) aaTATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACAcc                                                                                                                                                  <  2:405169/151‑1 (MQ=255)  atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACAccc                                                                                                                                                 >  2:846378/1‑151 (MQ=255)       tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGcc                                                                                                                                            >  2:553118/1‑151 (MQ=255)           gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGccccc                                                                                                                                         <  2:522571/150‑1 (MQ=255)               aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGa                                                                                                                                    <  2:997694/151‑1 (MQ=255)                         gAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGt                                                                                                                          >  1:404580/1‑151 (MQ=255)                         gAGGGAGCCAAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTACCGGTAAGGAGAAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCg                                                                                                                           <  1:161783/150‑1 (MQ=255)                          aGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGtt                                                                                                                         >  2:635502/1‑151 (MQ=255)                           gggAGCCAAAAAAATGTATGTGGATTTGCCAGCCATAATTATCAAAAAGATTAATGACGGCAGATATTTACCGGTAAGGAGAAAAGGCAGAGCTTTTCTTCACGTTCATTGGTCAAT‑GCCCCCCCCGGGCCCCCGCCGATACAAAGCGtt                                                                                                                         <  2:154995/150‑1 (MQ=255)                                cccAAAAAATGTATGTCGATTTGCCAGCCATCATTATTAAACAGTTTAAGGCCGGAAGATTTTTATCGGTAAGGGGGAAATGCAGAGTTTTTCGTCATTTTCATGGGTCAAT‑CCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                     <  1:287919/149‑1 (MQ=255)                                    aaaaaTTTATGTCGACTTCCAAGCCATCTTTATTCAAAAGATAAATGACGGCAGATTTTTATCGGTAGGGAGGAAATGCAGCGCTTTTCGTCCCTTTCATTGCTAAAT‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                               <  2:752324/150‑1 (MQ=255)                                     aaaaTTTATGTCGTCTTGCATGCCATATTTATTCAACAAATTAAGGACGGCATATATTTATGGGTAAGGGGAAAATGCAGAGCTTTCCGTCATTTTCAATGTCAAAT‑CCCCCACCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                              <  2:837235/151‑1 (MQ=255)                                       aaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                            >  2:59845/1‑151 (MQ=255)                                         tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                          >  1:103973/1‑151 (MQ=255)                                          gtaggacGACTTGCCACCCATCATTATTCAACAGATTAATGACGGCAGATATTTACCGGTAGGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGCCAAT‑GCCTCCCCCTTCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                         <  2:235583/145‑1 (MQ=255)                                             tgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                                       >  1:353688/1‑150 (MQ=255)                                             tgtCGACTTGCCAGCCATAATTATTCAACAGATAAATGACGGCAGATATTTACCGGAAAGGAGAAAATGCAGAGCTTTTCTTCAGTTTCATTGGTCATT‑GCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                            <  2:2632/145‑1 (MQ=255)                                                cgattgcccagccaTCATTATTCAACAGATTAATGACGGCATTTATTTAGGGATAGGGAAAAATTGCAGAGCTTTTCGTCCTGTTCAATGGTCAAC‑CCCCCCCCCTCACCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                    <  2:615497/143‑1 (MQ=255)                                                cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGTTTTTTATGGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAAT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                     <  2:698264/149‑1 (MQ=255)                                                cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                   >  2:591257/1‑151 (MQ=255)                                                   cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAGGGAGAAAATGCAGAGCTATTCGTCACTTTCAATGGTCAAT‑GCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                    <  1:985063/147‑1 (MQ=255)                                                    ttGCCAGCCATCATTATTCAACAGATTAATGCCGGCAGATATTTACCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTAAAT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                 <  2:337179/149‑1 (MQ=255)                                                         agccaTCATTATTCAACAAATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCGTTAAATGGTCAAT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                          <  1:495490/151‑1 (MQ=255)                                                           ccaTCATTATTCAACAGATTAATGCCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCCCGTTAAATGGTCAAT‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                        <  1:863654/151‑1 (MQ=255)                                                             atcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                      >  2:593883/1‑151 (MQ=255)                                                              ttattattAAACAAATAAATGCGGGAAGATATTTAGCGAAAGGGAGGAAATGGAGAATTTTTCGTCGCGTTAAGGGTTCAAC‑GCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                       <  1:657134/147‑1 (MQ=255)                                                                 ttattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTAAATGGTCAAT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                    <  1:963273/149‑1 (MQ=255)                                                                 ttattCAACAGATTAAGGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGAGTTATTCGTCATGTTCAATGGACAAT‑CCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                      <  2:872660/147‑1 (MQ=255)                                                                  tattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                  >  2:622942/1‑150 (MQ=255)                                                                   attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                >  1:561477/1‑151 (MQ=255)                                                                    ttCAACAGATTAATGACGGAAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAAT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                 <  1:648360/149‑1 (MQ=255)                                                                      cAATAGATAAATGGCGGCAGATATTTATCGGTAGGGAGAAAATGCAGAGTTTTTCGTCCCTTAAAATGGTCAAC‑CCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                              <  2:548160/150‑1 (MQ=255)                                                                            aTTAAGGACGGTAGTTATTTATCGGGAAGGAGGAAAGGAAGATTTATTCGTCCCGTTAAATGCTCAGT‑CCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                       <  1:27115/151‑1 (MQ=255)                                                                                   aCGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                >  1:618654/1‑151 (MQ=255)                                                                                   aCGGAAAATATTCATCGGGAAGGCAGAAATGGAGGTTTATTCCTCGGATGCAATGGTCATG‑GCCCCCCCGCGGCCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                <  1:598795/151‑1 (MQ=255)                                                                                    cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                               >  2:174502/1‑151 (MQ=255)                                                                                    cGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                               >  2:109683/1‑151 (MQ=255)                                                                                              ttATCGGTAAGGAGAAAATGCAGAGCTATTCTTCCCGTTCAATGGTCAAT‑GCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                      <  1:489123/150‑1 (MQ=255)                                                                                                   ggaagggaggaAATGCAGAGCTTTTCGTCCCTTTCAATGGTAAAT‑GCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                  <  2:887554/144‑1 (MQ=255)                                                                                                   ggTAGGGAGAAAATGCAGGCGTTTTGGTCA‑GTTTAAGGGTCAATCGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                       <  2:870857/144‑1 (MQ=255)                                                                                                    gTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                               >  1:928908/1‑151 (MQ=255)                                                                                                     tAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                              >  2:875928/1‑151 (MQ=255)                                                                                                      aaGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                              >  2:924883/1‑150 (MQ=255)                                                                                                        ggaggaAATGCAGAGTTATTCGTCCCGTTCAATGTTCAAT‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCt                                              <  1:848908/148‑1 (MQ=255)                                                                                                         gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                          >  1:520016/1‑151 (MQ=255)                                                                                                            gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                       >  1:687647/1‑151 (MQ=255)                                                                                                                        tATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                            >  2:900226/1‑150 (MQ=255)                                                                                                                        tATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                            >  1:80255/1‑150 (MQ=255)                                                                                                                           tCGGCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                        >  1:514018/1‑151 (MQ=255)                                                                                                                               cACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGcc                    >  2:478271/1‑151 (MQ=255)                                                                                                                                   ttCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                 >  2:997462/1‑150 (MQ=255)                                                                                                                                      aaTGGTCAAT‑GCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt              <  1:575844/150‑1 (MQ=255)                                                                                                                                      aaTGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg             >  2:335176/1‑151 (MQ=255)                                                                                                                                       aTGGTCAAT‑CCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt              <  2:687585/149‑1 (MQ=255)                                                                                                                                        tGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac           >  1:14407/1‑151 (MQ=255)                                                                                                                                          gTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat         >  1:500214/1‑151 (MQ=255)                                                                                                                                             aaT‑GCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt        <  1:739983/149‑1 (MQ=255)                                                                                                                                             aaT‑GCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg       <  1:378545/150‑1 (MQ=255)                                                                                                                                               t‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc    >  1:207781/1‑151 (MQ=255)                                                                                                                                                  ccACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg  <  2:140289/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                   AATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAAT‑GCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCG  >  NC_000913/2984293‑2984587

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