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

ATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984294‑2984591                                                                                                                                                     |                                                                                                                                                       atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACAccc                                                                                                                                                     >  2:297082/1‑151 (MQ=255) atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACAccc                                                                                                                                                     <  2:419631/151‑1 (MQ=255)  tatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCt                                                                                                                                                    >  1:605130/1‑151 (MQ=255)     gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGc                                                                                                                                                 >  2:958235/1‑151 (MQ=255)      tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGcc                                                                                                                                                >  2:329756/1‑151 (MQ=255)       cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGccc                                                                                                                                               >  2:159338/1‑151 (MQ=255)         tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGccccc                                                                                                                                             >  1:53359/1‑151 (MQ=255)          gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGccccc                                                                                                                                             <  2:626623/150‑1 (MQ=255)             aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCccgccg                                                                                                                                         <  2:698844/151‑1 (MQ=255)              aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCccgccg                                                                                                                                         >  1:209119/1‑150 (MQ=255)                 gAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATAc                                                                                                                                     >  1:324666/1‑151 (MQ=255)                          gggAGCCAAGAAAATGTATGTCGATTTCCCACCCATCATTATACAACAGATAAATGCCGGCAGATATTTATGGGAAAGGAGAAAATGCAGATCTTTTCGTCACGTTCATGGGTCAA‑CCCCCCCCCCTGCCCCCCGCCGATACAAAGCGtt                                                                                                                             <  1:720715/150‑1 (MQ=255)                               cccAAAAAATGTATGTCGCCTTGCCACCCATTATTATTCAACAGATTAATGCCGGCAAATTTTTACCGGTAAGGAGAAAATGCAGAGCTATTCCTCCCGTTCAATGGTAAA‑CCCCACACCCTGCCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                       <  1:266025/151‑1 (MQ=255)                               ccaaaaaaaTGTATGTCGCCTTGCCAGCCATCATTATTAAACAAATAAATGACGGAAGATTTTTATCGGTAGGGAGAAAATGCAGAGTTTTTCGTCACGTTCAATGGTCAA‑TCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                         <  2:324679/144‑1 (MQ=255)                                 cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                     >  2:750618/1‑151 (MQ=255)                                  aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                    >  1:384767/1‑151 (MQ=255)                                   gAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                   >  2:284783/1‑151 (MQ=255)                                    aaaaTGTATGTCGACTTGCCAGCCATCATTATTAAAAAAATTAATGACGGAAGATTTTTATCGGAAGGGAGAAAATGCAGATCTTTTCGTCCCGTTAAATGGTCAA‑TCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                   <  1:318570/150‑1 (MQ=255)                                        tgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGTGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCgg                                                                                                              >  1:450438/1‑151 (MQ=255)                                          tatgtCTATTTGCCACCCATTATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGCAGACTTTTTCGTCCCTTTCAATGGCCAA‑TCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGct                                                                                                               <  1:938857/148‑2 (MQ=255)                                          tatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                            >  1:759333/1‑151 (MQ=255)                                               cGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                       >  2:87860/1‑151 (MQ=255)                                                gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                      >  2:823642/1‑151 (MQ=255)                                                   ttGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAAGGCAGATTTTTTCGTCACGTTAAATGGTCAA‑TCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                    <  2:617670/150‑1 (MQ=255)                                                          ccaTAATTATTCAACAGATTAATGCGGGCAGATATTTACGGGAAAGGAGAAAATGCAGAGCTTTTCGTCATTTTAATTGGTCAA‑TGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                   <  2:475864/144‑1 (MQ=255)                                                                ttatcaaaCAAATTAAGGAGGGCAAATTTTTATGGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAGTGGTAAA‑CGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  1:322189/141‑1 (MQ=255)                                                                   ttCAACAGATTAATGACGGAAAATATTTATCGGAAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAAGGGTAAA‑GCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  2:441820/149‑1 (MQ=255)                                                                      aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                   >  1:663089/1‑148 (MQ=255)                                                                      aaCAAATTAAGGCAGGCATATTTTTCGCGGAAAGGAGAAAGTGCAGGGTTTTTCGTCACTTAAAGTGGTCAG‑TGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                  <  1:297069/149‑1 (MQ=255)                                                                          gATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                             >  2:624811/1‑150 (MQ=255)                                                                           aTTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAc                                                                           >  1:404070/1‑151 (MQ=255)                                                                             tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          >  1:589971/1‑150 (MQ=255)                                                                                     gCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                 >  1:611970/1‑151 (MQ=255)                                                                                          tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            >  2:70295/1‑151 (MQ=255)                                                                                             ttATCGTTAAGGAGAAAAGGCAGAGTTTTTCGTCACTTTAAAGGGTAAA‑TGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCgg                                                          <  2:611994/150‑1 (MQ=255)                                                                                                  ggTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                    >  1:752349/1‑151 (MQ=255)                                                                                                  ggTAAGGAGAAAATGTAGAGCTATTCGTCACGTCAAATGGTCAA‑TCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                     <  2:410158/150‑1 (MQ=255)                                                                                                    tAAGGAGAAAAGGCAGAGCTTTTCGCAACGTTTAAGGGTCAC‑TGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTGCGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAg                                                       <  1:592660/146‑1 (MQ=255)                                                                                                     agggagaaAATGCAGAGTTTTTCGTCACGTTAAATGGTCAA‑TCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                     <  1:459947/145‑1 (MQ=255)                                                                                                     aaGGAGGAAATGCAGACTTTTTCGTCTCTTTAAGGGCTCAA‑TGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                        <  2:96908/144‑1 (MQ=255)                                                                                                      aggaggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                >  1:702692/1‑151 (MQ=255)                                                                                                       ggagaaAAGGCAGAGCTATTCGTAAC‑TTAAATTGGCAAGCGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                               <  1:230329/151‑1 (MQ=255)                                                                                                        gaggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                              >  2:367708/1‑151 (MQ=255)                                                                                                            aaaTGCAGATTTTTTCGTCACGTTAAAGGGTCAC‑TCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                             <  2:643045/148‑1 (MQ=255)                                                                                                            aaaTGAAGGGTTTTTCGTCACGTACAGGGATAAA‑TCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                          <  1:79011/151‑1 (MQ=255)                                                                                                            aaaGGCAGAGCTTTTCGTCCCGTTAAATGGTAAA‑TCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                            <  2:697546/149‑1 (MQ=255)                                                                                                                         ttcgtaattttaagtGGACAA‑TGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                  <  2:483346/132‑1 (MQ=255)                                                                                                                         ttCGTCACGTTCAATGGTCAA‑TCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                               <  1:968428/149‑1 (MQ=255)                                                                                                                          tCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            >  2:911813/1‑151 (MQ=255)                                                                                                                                  ttCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  2:421573/1‑150 (MQ=255)                                                                                                                                  ttCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  2:839468/1‑150 (MQ=255)                                                                                                                                    cAATGGTAAA‑CGCCCCCCCATGGCCACCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     <  2:450451/148‑1 (MQ=255)                                                                                                                                     aaggttCAA‑TCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcat                   <  2:846969/144‑1 (MQ=255)                                                                                                                                       ttttCAA‑TCCCCCACCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAATCCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattg                 <  2:629544/146‑1 (MQ=255)                                                                                                                                           cAA‑TCCCACCCCCGGCCCCCCGCCGATACAAAGCGTCGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           <  1:329742/151‑1 (MQ=255)                                                                                                                                             a‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCattgacattgac         >  1:817092/1‑151 (MQ=255)                                                                                                                                               tGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt            <  1:800535/147‑1 (MQ=255)                                                                                                                                               tGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAcc        >  2:43291/1‑151 (MQ=255)                                                                                                                                                 ccACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg      >  1:170581/1‑151 (MQ=255)                                                                                                                                                  cacaCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCg      >  1:518378/1‑150 (MQ=255)                                                                                                                                                    cacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACg   >  1:950171/1‑151 (MQ=255)                                                                                                                                                    cacCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACg   >  2:584098/1‑151 (MQ=255)                                                                                                                                                     acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:758141/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAA‑TGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984294‑2984591

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