Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	56.7%	55.9 / 91.2	60	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (20/6);  new base C (0/34);  total (20/40)
Fisher's exact test for biased strand distribution p-value = 5.49e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-11

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591                                                                                                                                                          |                                                                                                                                                       gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCacac                                                                                                                                                       <  2:456860/151‑1 (MQ=255)   atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACAccc                                                                                                                                                     <  1:111860/151‑1 (MQ=255)     atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTg                                                                                                                                                   <  2:350799/151‑1 (MQ=255)      tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTg                                                                                                                                                   >  2:938378/1‑150 (MQ=255)       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                                                            >  1:650075/1‑151 (MQ=255)                                                                                              tatTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            >  2:859028/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                                >  1:755847/1‑150 (MQ=255)                                                                                                                            tATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                >  2:28313/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)                                                                                                                                                    g‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc       >  1:205943/1‑151 (MQ=255)                                                                                                                                                       caccCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa     <  1:811337/151‑1 (MQ=255)                                                                                                                                                          acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  1:504882/1‑151 (MQ=255)                                                                                                                                                          |                                                                                                                                                       GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCC‑AGCCATCATT‑ATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGC‑AGAGCTATTCGTCACGTTCAATGGTCAATG‑CCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591

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