Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,447	0	G	C	46.8%	72.9 / 23.4	62	G382G (GGC→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (21/12);  new base C (1/28);  total (22/40)
Fisher's exact test for biased strand distribution p-value = 4.64e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCAT  >  NC_000913/2984298‑2984595                                                                                                                                                          |                                                                                                                                                     gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTgg                                                                                                                                                      <  1:300832/150‑1 (MQ=255)  tCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGcc                                                                                                                                                    >  2:68046/1‑151 (MQ=255)      gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGccccc                                                                                                                                                 <  2:898847/150‑1 (MQ=255)      gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGccccc                                                                                                                                                 <  1:668270/150‑1 (MQ=255)      gagaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCccg                                                                                                                                                <  2:512397/151‑1 (MQ=255)         aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCccgccg                                                                                                                                             >  2:60060/1‑151 (MQ=255)                      gggAGCCAAAAAAATGTATGTCGACTTGCC‑ACCCATCTTT‑ATTCAACAGATTAATGACGGCAGATTTTTATCGGTAAGGAGGAAATGC‑AGAGCTTTTCGCCACGTTAAAGGGTAAATG‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTg                                                                                                                                <  1:300087/151‑1 (MQ=255)                               gAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                       >  2:79667/1‑151 (MQ=255)                                 aaaTGTATGCCGA‑TTTCC‑AACCAATATTAATTCAAAAGTTTAAGGTGGGCAGATTTTTTTCGGAAGGAAGGAAAGCC‑AGGGCTATTCTCCAGGTTAAATGTTAACTG‑CCACCCCCCGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTg                                                                                                                     <  1:390724/151‑1 (MQ=255)                                      tatgtCGACTTGCC‑AGCCATCATT‑ATCAAACAGATTAATGCCGGAAGATATTTATCGGTAAGGAGGAAATGC‑AGAGTTATTCGTCACTTTAAATGGTCAACG‑CCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                                 <  2:297631/150‑1 (MQ=255)                                       atgtCGA‑TTGCCAAGCCATCATT‑ATCAAAAAGATAAATGCCGGAGAATATTCAAGGGAAAGGGGAAAAGGC‑AGGTCTTTTCGCCACGTTCAATGGTAACTG‑CCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                                <  2:313796/150‑1 (MQ=255)                                        tgtCGCCTTGCC‑AGCCATTATT‑ATTCAAAAGATTAAGGAGGGAAGATATTTATCGGTAAGGAAGAAATGC‑AGGGGTTTTCGTCGCTTTCAAGGGACAGTC‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCa                                                                                                               <  2:355354/150‑1 (MQ=255)                                        tgtCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                              >  1:287727/1‑151 (MQ=255)                                         gccgccTTGCC‑CGCCATCTTT‑ATCCAAAAAATTAATGCCGGAAAATTTTTATCGGTAAGGAGAAAATGC‑AGACTTTTTCGTCCCGTTAAGTGGTAACTG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                              <  1:256707/145‑1 (MQ=255)                                            gACTGGCC‑ACCCATCATT‑ATTCAAAAGATTAATGACGGCAGATATTTATCGGTAGGGAGAAAATGC‑AGATCTTTTCGTCACGTTAAATGGTCAACC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                          <  1:772322/151‑1 (MQ=255)                                               ttGCC‑AGCAATCATT‑ATTAAACAGATTAATGCCGGCATTTATTAATCGGTAGGGAGGAAATCA‑AGAGCTTTTCGCCCCGTTAATTGGT‑AATGCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                                 <  1:656917/141‑1 (MQ=255)                                                tccc‑acccaTATTT‑ATCCAAAAAATTAATGACGAAAGATTTTTATGGGTAAGGAGAAAATGC‑AGAGTTTTTCGTCACGTTCAATGGTCAATC‑CCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                       <  1:802051/148‑1 (MQ=255)                                                        catcatT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                               >  1:32894/1‑151 (MQ=255)                                                            att‑attCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                           >  2:811399/1‑151 (MQ=255)                                                             tt‑attCAACAAATTAATGACGGCAGATATTTATGGGAAGGGAGAAAATGA‑AGAGCTTTTCGTCACGTTAAATGGTCAATG‑CCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                                  <  2:556372/143‑1 (MQ=255)                                                                attCAACAGATTAATGCCGGCAAATATTTA‑CGGTTAGGGAGAAAAGCAAGAGCTTTTCGTCCCGTTAAAGGGTCAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                          <  1:60453/149‑1 (MQ=255)                                                                    aagaGATAAATGCCGCAAGATATTTTGGGGTAGGGAGAAAAGGC‑AAAGCTTTCCGTCACCTTCAATGG‑CAACGACCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                      <  2:346558/146‑1 (MQ=255)                                                                    aaagatttAACGCCGGCACATATTTATCGGTAAGGAGAAAACCC‑AGCTCTGTTCGTCCCTTTAAGCGCTAAAGG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                                    <  1:680287/145‑1 (MQ=255)                                                                    aaaaaaTTAATGACGGAAGATATTTACCGGTAAGGAGAAAAGGA‑AGGGTTTTTCGTCACTTTAAAGGTTCAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCaga                                                                                    <  1:959112/146‑1 (MQ=255)                                                                    aaCAGATTAATGACGGCATATTTTTATCGGTAAGGAGAAAATGC‑AGAGCTTTTCGTCCCGTTAAATGGTCAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                      <  1:472269/149‑1 (MQ=255)                                                                    aaCAGATTAATGACGGCAGATATTTAGCGGAAAGGAGGAAATGC‑AGATCTTTTCGTCACGTTCAAGGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                        <  2:388802/147‑1 (MQ=255)                                                                       aGATTAATGACGGCAGATATTTATCGGTAAGGAGAAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                                  <  2:280025/150‑1 (MQ=255)                                                                          ttAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                              >  2:225657/1‑151 (MQ=255)                                                                            aaTGACGGCAAATTTTTATCGGAAAGGAGAAAACGC‑AGAGTTTTTCGTCACTTTAAATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                             <  1:292793/150‑1 (MQ=255)                                                                             aTGAGGGCAAATTTTAATCGGTAAGGAGAAAATGC‑AGACCTTTTCGTCACTTTAAATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                             <  2:726886/149‑1 (MQ=255)                                                                                aCGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAtt                                                                        >  2:16959/1‑151 (MQ=255)                                                                                  ggCAGTTATTTATCGTTAAGGAGAAAAGGC‑AGACTTTTTCGCCCCTTTCAATGGTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                       <  2:830603/150‑1 (MQ=255)                                                                                  ggCAGATATTTATCGGAAGGGAGAAAAGGC‑AGGCCTTTTCGAAACGTTAAGTGGACAAGG‑CCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                          <  1:890715/147‑1 (MQ=255)                                                                                     aGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                    >  1:943224/1‑150 (MQ=255)                                                                                        tttttATCGGAAAGGAGGAAATGC‑AGAGCTTTTCGTCACGTTAAATGGTCAATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTG‑CAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                                 <  2:684322/147‑1 (MQ=255)                                                                                        tatTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                                >  2:859028/1‑151 (MQ=255)                                                                                        tatTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                                >  1:650075/1‑151 (MQ=255)                                                                                          tttATCGGTAAGGAGAAAATGC‑AGAGCTTTTCGTCACGTTAAATGGTCAATC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                                <  2:195310/149‑1 (MQ=255)                                                                                          tttATCGGAAAGGAGGAAATGC‑AGAGCTTTTCGTCACGTTCAGTGGTCAATG‑CCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                       <  1:503766/142‑1 (MQ=255)                                                                                                 gTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                       >  2:874727/1‑151 (MQ=255)                                                                                                    aggaggAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                    >  2:110801/1‑151 (MQ=255)                                                                                                     gaagaataTGC‑AGAGCTTTTCGGCACGCTTAGTGCCAAACC‑CCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                    <  1:462213/143‑1 (MQ=255)                                                                                                          aaaGGC‑AGAGCTTTTCTCCACGTTAAATGGTCAACC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                              <  1:428719/151‑1 (MQ=255)                                                                                                            aTGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAc                                            >  2:122200/1‑151 (MQ=255)                                                                                                             tGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGAcc                                           >  2:394552/1‑151 (MQ=255)                                                                                                                      tATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                    >  2:28313/1‑150 (MQ=255)                                                                                                                      tATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                    >  1:755847/1‑150 (MQ=255)                                                                                                                         tCGTCACTTTCAATGGTCAATG‑CCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                                 <  2:849568/150‑1 (MQ=255)                                                                                                                          cGTCACGTTCAATGGTCAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTCGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                                 <  2:120835/149‑1 (MQ=255)                                                                                                                           gccccGTTAAATGGCAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                               <  1:87813/146‑1 (MQ=255)                                                                                                                             cacttacAATGGCAACTG‑CCCCCCCTTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                                 <  1:493535/140‑1 (MQ=255)                                                                                                                                tttaaaTGGTCAATG‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcgc                               <  1:470249/141‑1 (MQ=255)                                                                                                                                 ttCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                         >  2:180812/1‑150 (MQ=255)                                                                                                                                       ggTAAATC‑CCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt                      <  1:700475/147‑1 (MQ=255)                                                                                                                                        gTCAATG‑CCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacat                 <  2:144710/151‑1 (MQ=255)                                                                                                                                            agtcccACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt                <  1:412279/145‑1 (MQ=255)                                                                                                                                            acccccccccccTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattga                    <  1:60052/136‑1 (MQ=255)                                                                                                                                              g‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc           >  1:205943/1‑151 (MQ=255)                                                                                                                                                 cccccccGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa         <  2:607466/147‑1 (MQ=255)                                                                                                                                                 caccCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa         <  1:811337/151‑1 (MQ=255)                                                                                                                                                    acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc      >  1:504882/1‑151 (MQ=255)                                                                                                                                                      ccTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTc    >  1:423820/1‑151 (MQ=255)                                                                                                                                                        tGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCAt  >  1:387679/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                     GTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGCTCAT  >  NC_000913/2984298‑2984595

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