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

TGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984297‑2984591                                                                                                                                                 |                                                                                                                                                        tGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                   >  2:242251/1‑151 (MQ=255)  gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTgg                                                                                                                                                   <  1:515088/150‑1 (MQ=255)        agaATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccgc                                                                                                                                           >  2:920859/1‑151 (MQ=255)          aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccgccg                                                                                                                                         <  1:1215060/151‑1 (MQ=255)          aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCccgccg                                                                                                                                         <  2:469855/151‑1 (MQ=255)           aTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGa                                                                                                                                        <  2:4398/151‑1 (MQ=255)                aaaaGGAGGGACCCAAAAAAATTTATGCCGACTTCCCACCCATCATTATTCAACAGATTAATGACGGAAGATATTTACCGGAAAGGAGAAAACGCAGAGTTATTCGTCACGTTCAATGATCAATCCCACCCCCTGC‑CCCCCGCCGATACaa                                                                                                                                   <  2:600487/151‑1 (MQ=255)                            ccaaaaaaaTGTATGTCGTCTTGCCACCCATCATTATACAACAGATAAATGACGGCAGATTTTTATCGGTAAGGAGGAAATGCAGAGCTTTTCGTCACTTTAAATGGTCAATCCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                        <  1:857948/145‑1 (MQ=255)                              cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAACGTTGCCAGTcc                                                                                                                     >  1:1179332/1‑151 (MQ=255)                                 aaagTGTATGTCGGCTTCCGGGACATCATTATTAAACAAATTAATGGCGGCAGATATTTAAGGGAAGGGAGAAAAGGCAGACCTTTTCGTCACTTTAAATGGTCAGTACCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                   <  1:249189/150‑1 (MQ=255)                                 aaaaTGTATGTCGTCTGCCAAGCCTTAATTTTTAAACAGATTAATGACGGCAGATTTTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCCCGTTCAATGGTCAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  2:634815/151‑1 (MQ=255)                                    atgtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                               >  2:229058/1‑151 (MQ=255)                                     tgtttgcCGATTTGCCAGCCAACATTATTAAAAAAATTAAGGCCGGAAGATTTTTATCGGAAAGGAGAAAATGCAGACTTTTTGGTCACGTTCAATGGTCAATGCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                               <  2:1074094/150‑1 (MQ=255)                                      gtatgtCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCggg                                                                                                             >  1:771515/1‑151 (MQ=255)                                            cgccTTGCCAGCAATCTTTATTAAACAAATTAAGCCCGGCATATATTTATCGGAAAGGAGAAAACGCGGATCTATTCGCCTCGTCCAATGCTCAGTCCCCCCCCGGCC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAt                                                                                                          <  2:783736/145‑1 (MQ=255)                                             gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                      >  2:946919/1‑151 (MQ=255)                                              aCTTGCCAGCCATCATTATTAAACAGATAAATGACGGCAGATATTTATGGGTAAGGAGGAAATGCAGAGCTTTTCGTCACGTTCAGGGGACAACGCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                      <  2:630594/150‑1 (MQ=255)                                               cTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGtt                                                                                                    >  1:650682/1‑151 (MQ=255)                                                    cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               >  2:844640/1‑151 (MQ=255)                                                    cagccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               >  1:969052/1‑151 (MQ=255)                                                     atcaataattattAAACAATTTAAGGACGGCAGATATTTTTGGGTAAGGAGAAAATGCATATTTATTCGCCACTTTAAATGGTAAAGGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCa                                                                                              <  2:1152966/144‑1 (MQ=255)                                                            attattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       >  2:741687/1‑151 (MQ=255)                                                              tattCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                      >  2:71508/1‑150 (MQ=255)                                                                ttaaaaaGATTAATGCGGGAAGATATTTATGGGTAAGGAGAAAATGAAGATTTTTTCGTCACTTTCAAGGTTAAATGCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  1:899780/143‑1 (MQ=255)                                                                tcaaaCAAATTAAGGACGTCAGATTTTTATCGGTAGGTAGAAAACGCAGATCTGTTCGTCAGTTTAAATGGTCAATCCCACCCCCTGG‑CCTCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  2:42263/146‑1 (MQ=255)                                                                   aaCAGATTAATGACGGCAGATATTTATGGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATGCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAc                                                                                    <  1:752596/147‑1 (MQ=255)                                                                   aaCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                   >  1:710857/1‑148 (MQ=255)                                                                    aCAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                 >  1:655700/1‑149 (MQ=255)                                                                          tAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          >  2:1206974/1‑150 (MQ=255)                                                                           aaTGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          >  2:707406/1‑149 (MQ=255)                                                                                 ggaatatatTTATCGGTAAGGAGAAAATGCAGAGCTTTTCTTCACGTTAAATGGTCAATCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                      <  2:302790/142‑1 (MQ=255)                                                                                 ggCAGTTATTTATCGGTAAGGAGAAAAGGCAGATTTTTTCGTCCCGTAAAGTGGTAAATGCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                      <  1:920858/147‑1 (MQ=255)                                                                                 ggCAGATATTTATCGGAAAGGAGAAAATGCAGAGTTATTCGTCACGTTCAAGGTTCAACCCCCCCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                   <  2:373438/150‑1 (MQ=255)                                                                                 ggCAGATATTTAGCGGTAAGGAGGAAATGCAGAGCTATTCGTCACTTTCAATGGTCAATGCCACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                   <  2:255911/150‑1 (MQ=255)                                                                                 ggCAGATATTTAGCGGTAAGGAGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGa                                                                      <  1:148904/147‑1 (MQ=255)                                                                                      atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                             >  2:72338/1‑151 (MQ=255)                                                                                      atatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAc                                                             >  1:1163000/1‑151 (MQ=255)                                                                                       tatTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAAcc                                                            >  2:984805/1‑151 (MQ=255)                                                                                         tttaggggTAAGGAGAAAATGCGGCGTTTTTCGTCCCGTCAAATGTTCAATCCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCaa                                                              <  2:424374/141‑1 (MQ=255)                                                                                         tttATGGGAAAGGAGGAAAGGCAGAGTTTTTCGTCCCTTTCAATGGTAAACCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                           <  2:875038/150‑1 (MQ=255)                                                                                           tATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                         >  1:367639/1‑150 (MQ=255)                                                                                              cGGGAAGAAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCGGAGAAACCAGGATTCGGCAACCGGa                                                         >  2:569013/1‑147 (MQ=255)                                                                                               ggaaaGGAGAAAAGGCAGAGTTTTTCTTCACGTTCAAGGTTCAATCCCACCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgc                                                      <  1:156974/146‑1 (MQ=255)                                                                                                    ggagaaAAGGCAGAGCTTTTCGTCACGTTAAATGGTAAAGCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTa                                                 <  2:832068/149‑1 (MQ=255)                                                                                                      aggaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                             >  2:896843/1‑151 (MQ=255)                                                                                                        gaAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGt                                           >  2:458876/1‑151 (MQ=255)                                                                                                         aaaTGCAGAGTTTTTCGCCCCTTTAAATGGTCAATCCAACCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                          <  1:752031/151‑1 (MQ=255)                                                                                                              cAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  2:1233041/1‑150 (MQ=255)                                                                                                                 agCTTTTCGCCCCTTTCAATGGTAAATCCCCCCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      <  1:815192/147‑1 (MQ=255)                                                                                                                    ttttCGTCACTTTAAAT‑GTCAATTCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGc                                 <  2:784763/147‑1 (MQ=255)                                                                                                                      ttCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                              >  2:1186571/1‑150 (MQ=255)                                                                                                                      ttCGTCACGTAAAATGGTCACTCCCACCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                               <  2:1031433/149‑1 (MQ=255)                                                                                                                       tCGTCACTTTCAAGGGTCAACCCCCCCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                              <  2:771507/149‑1 (MQ=255)                                                                                                                       tCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATcg                            >  1:377743/1‑151 (MQ=255)                                                                                                                              gTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  2:1093322/1‑151 (MQ=255)                                                                                                                               ttaaaTGGTCAACCCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCgcgc                         <  2:97453/143‑1 (MQ=255)                                                                                                                                       tCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt            >  1:792152/1‑151 (MQ=255)                                                                                                                                         attGCCCCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           <  1:749480/148‑1 (MQ=255)                                                                                                                                         aaTGTCACACCCTGG‑CCCCCGCCGATACAAAGCGGTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGAGACCGAGTGCGATCGCGCCACcattgacattg           >  2:835659/1‑150 (MQ=255)                                                                                                                                         aaTGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           >  1:319881/1‑150 (MQ=255)                                                                                                                                         aaTGCCACACCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           <  2:417120/150‑1 (MQ=255)                                                                                                                                              caccCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGa     <  1:47708/151‑1 (MQ=255)                                                                                                                                                 acCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  1:231729/1‑151 (MQ=255)                                                                                                                                                 |                                                                                                                                                        TGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCACACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984297‑2984591

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