Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,679,304	0	A	C	43.1%	99.5 / 20.8	72	Y293S (TAC→TCC)	ydgH	DUF1471 family periplasmic protein
Reads supporting (aligned to +/- strand):  ref base A (30/11);  new base C (0/31);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 2.13e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGC  >  NC_000913/1679157‑1679451                                                                                                                                                    |                                                                                                                                                    aTGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTAcc                                                                                                                                                   >  2:182550/1‑150 (MQ=255)     cACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGCTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATc                                                                                                                                              <  2:170129/151‑1 (MQ=255)         aaaGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCAccc                                                                                                                                          >  2:1050658/1‑151 (MQ=255)         aaaGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCAccc                                                                                                                                          <  1:510037/151‑1 (MQ=255)                   aaCTAAAAAAAGCGACCGCGGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCAAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCa                                                                                                                                <  1:585448/151‑1 (MQ=255)                      tGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAgg                                                                                                                              >  2:80143/1‑150 (MQ=255)                        aaCAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAAc                                                                                                                           >  2:52328/1‑151 (MQ=255)                          cAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACg                                                                                                                          >  1:519602/1‑150 (MQ=255)                           aaaGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTg                                                                                                                        >  1:777723/1‑151 (MQ=255)                                  ccGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCGGGCAACTCCGGCAACATGACCGAATTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGtaa                                                                                                                    <  2:1068547/148‑1 (MQ=255)                                   cGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAAc                                                                                                                >  2:231936/1‑151 (MQ=255)                                       gCGATGATGGTCCCGTTCGACATCACAAAATTCTCGGGCAACTAGGGCAACATGCCCGAATTCTCCTATCAGGTTGCGAACCGTGCCGCGAAGAAAGGTGCTAAGTTTTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGa                                                                                                            <  2:140590/151‑1 (MQ=255)                                       gCGATGATGGTCCCGTTCGACAGCATCAAATTCTCGGGCAACTACGGCAACATGACCGAAGTCCCCTATCAGGTTGCGAAACGTGCCGCAAAGAAAGGTGCTAAGTATTACCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCt                                                                                                              <  2:484845/149‑1 (MQ=255)                                        cGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAc                                                                                                           >  1:1069354/1‑151 (MQ=255)                                         gatgatGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGAcc                                                                                                          >  2:1056933/1‑151 (MQ=255)                                               ggTCCCGTTCGACAGCATCAATTCCTCTGGCATTTACGGCAATATGGCCGAAGTCTCCTACCGGTTTGCGAAACTGCCCGCGAAAAAAGGGGTTAACTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgc                                                                                                    <  2:160077/151‑1 (MQ=255)                                                gTCCCGTCCGACAGCATCAAATTCTCTGGCAACTACGGAAACAGGACCGAAGTCCCCTATCAGTTTGCAAAACGGGCCGCAAAGAAAGGTGCAAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAg                                                                                                     <  1:1083052/149‑1 (MQ=255)                                                  cccTTTCGACAGCACAAAATTCTCTGGCAACTCCGGAACAAGGCCCGAAGTCTCCTATCAGGTTGCAAAACGTCCCGCGAAAAAAGGTGCTAAGTTTTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCa                                                                                                 <  2:202961/151‑1 (MQ=255)                                                  cccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                                                                                  <  2:441120/150‑1 (MQ=255)                                                   ccGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAg                                                                                                <  2:243486/151‑1 (MQ=255)                                                   ccGTTCGACAGCATAAAATTCTGTGGCAATTGCGGCAACATGCCGGAAGTCTCCTATCAGTTGGCAAAACGTGCCGCAAAAAAAGGTGGTAAGTATTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAg                                                                                                <  2:914097/151‑1 (MQ=255)                                                    cGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAgcgc                                                                                                  >  1:160601/1‑148 (MQ=255)                                                         gACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCCCCTATCAGGTTGCGAAACGTGCCGCAAAGAAAGGTGCTAATTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATc                                                                                             <  2:983384/148‑1 (MQ=255)                                                            accATCAATTTCTCTGGCAATTCCGGAAAAAGGACCGAATTCTCCTATCAGGTTGCAAAACGTGCGGGAAAAAAGGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGtata                                                                                       <  1:978161/149‑1 (MQ=255)                                                              cATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataa                                                                                      >  1:414211/1‑150 (MQ=255)                                                               ataaaaTTCTCGGGCAACTCCGGAAACATGACGGAAGTCTCTTATCAGGTTGCGAACCGTGCCGCGAAGAAAGGTCCTAAGTATTCCCACAACACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataa                                                                                      <  2:497398/146‑1 (MQ=255)                                                                  aaaTTCTCTGGCACCGCCTGCAAAATGCCCGAATTCTCCTCCCAGTTTGGCAACGGTGCCGCGAAAAAGGGTGCTAAGTATTCCCAATCCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAAc                                                                                 <  2:1052853/151‑1 (MQ=255)                                                                  aaaTTCTCTGGCAACTCCGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGta                                                                                         <  2:825090/143‑1 (MQ=255)                                                                  aaaTTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGTGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTataaataa                                                                                  >  1:170528/1‑150 (MQ=255)                                                                     ttCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTAt                                                                              >  1:173141/1‑151 (MQ=255)                                                                     ttCTCTGGCAACTACGGCAACAGGCCCGAAGTCTCTTATCAGTTTGCAAAACGGCCCGCGAAGAAAGGTGCTAAGTATTACCACTCCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTa                                                                               <  1:988165/150‑1 (MQ=255)                                                                          tGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGt                                                                         >  1:545091/1‑151 (MQ=255)                                                                           ggCAACTACGGAAACATGCCCGAAGTCTCCTATCAGGTTGCAAAACGTGCCGGGAAGAAGGGTGCTAAGTATTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt                                                                        <  1:762401/151‑1 (MQ=255)                                                                           gccAACTGCGGCAACAGGACCGAAGTCCCCTTGCAGGTTGCAAAACGGACCGCGAAGAAGGGTGCTAATTATTCCCACATCGTCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCa                                                                           <  1:715850/146‑1 (MQ=255)                                                                             cAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGtt                                                                        >  1:963952/1‑149 (MQ=255)                                                                              aactccgaaaaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGGGCCGCGAAAAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTAc                                                                      <  1:347257/141‑1 (MQ=255)                                                                                   cGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCtat                                                                >  1:466110/1‑151 (MQ=255)                                                                                   cGGAAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAAAAAGGTGTTAAGTATTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTcc                                                                   <  2:937438/148‑1 (MQ=255)                                                                                     gCAACATGACCAAAGTCTCCTAACAGGTTGGAAAACGTGCCGAGAAGAAAGGTGCAAATTTTTCCCACACCAGCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTcc                                                                   <  2:1184888/146‑1 (MQ=255)                                                                                       aaCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa                                                              >  1:547851/1‑149 (MQ=255)                                                                                        aCATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATaa                                                              >  1:144266/1‑148 (MQ=255)                                                                                                 aattctcTTATCAGTTGGCGAAACGGGCCGCAAAAAAAGGTGCTAATTATTCCCCCACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGAt                                                   <  1:976114/147‑1 (MQ=255)                                                                                                 aaGTCTCCTATCGGGTTGCAAAACGTGCGGCAAAGAAGGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCtt                                                      <  2:316987/147‑1 (MQ=255)                                                                                                 aaGTCTCCTATCAGGTTGCAAAACGTGCCGCAAAGAAAGGTGCTAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATg                                                  <  1:815062/151‑1 (MQ=255)                                                                                                     ctcCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgc                                              >  2:952575/1‑151 (MQ=255)                                                                                                       ccTATCAGGTTGCGAACCGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgc                                              <  2:841546/149‑1 (MQ=255)                                                                                                       ccTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCt                                            >  1:824853/1‑151 (MQ=255)                                                                                                       ccTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCt                                            >  2:314031/1‑151 (MQ=255)                                                                                                       ccTATCAGGTGGCGAAACGTGCCGCGAAAAAGGGTGCAAAGTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATgccgcc                                             <  1:576782/150‑1 (MQ=255)                                                                                                             gggTTGCGAAGCGTGCCGCGAAGAAAGGTGCTAAGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCa                                      <  1:665858/150‑1 (MQ=255)                                                                                                                ttGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCAttt                                   >  2:1001069/1‑151 (MQ=255)                                                                                                                ttGCAAAACGTGCCGCGAAGAAAGGTACTAATTTTTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACACCTATAAGGCGGCTTGATGCCGCCTTTCGCAtt                                    <  1:10420/150‑1 (MQ=255)                                                                                                                     aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGt                              >  1:379805/1‑151 (MQ=255)                                                                                                                     aaaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGt                              >  2:432521/1‑151 (MQ=255)                                                                                                                      aaCGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTc                             >  1:1278236/1‑151 (MQ=255)                                                                                                                        cGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTcac                           >  1:10427/1‑151 (MQ=255)                                                                                                                         gggCCGAGAAGAAAGTTGCTAATTTTTCCCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTcaca                          <  1:1031803/149‑1 (MQ=255)                                                                                                                            ccGAGAAGAAGGGTGATAAGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaa                         <  1:699512/149‑1 (MQ=255)                                                                                                                               caaagaaAGGTGTTAATTATTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACaaa                        <  2:1002857/145‑1 (MQ=255)                                                                                                                                 aagaagGGTGCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAtt                      <  1:184633/147‑1 (MQ=255)                                                                                                                                 aagaaAGGTGCTAAGTTTTCCCACACCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCTCCTTTCGCATTTTTAGTCACAAATTTc                    <  1:619057/149‑1 (MQ=255)                                                                                                                                 aagaaAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTccc                  >  1:1269657/1‑151 (MQ=255)                                                                                                                                 aaaaaaGGTGCAAAGTATTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAtt                      <  2:486831/144‑1 (MQ=255)                                                                                                                                 aaaaaaGGTCCTAAGTATTCCCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAAttt                     <  1:395277/145‑1 (MQ=255)                                                                                                                                    aaaGGTGCTAATTTTTACCACACCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCcgc               <  1:1278120/151‑1 (MQ=255)                                                                                                                                      aGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCcgcg              >  2:1191701/1‑150 (MQ=255)                                                                                                                                          gCTAAGTATTCCCCCATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCAt         <  2:515115/151‑1 (MQ=255)                                                                                                                                             aaGTTTTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc      <  2:1041878/151‑1 (MQ=255)                                                                                                                                             aaGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgc      <  2:1175914/151‑1 (MQ=255)                                                                                                                                             aaGTATTCCCACATCCCCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtg       <  2:855005/150‑1 (MQ=255)                                                                                                                                                tATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatg   >  1:372938/1‑151 (MQ=255)                                                                                                                                                 aTTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATtgcatgc  >  1:300129/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATGGCACTAAAGTCGAAGAACTGAACAAAGCGACCGCAGCGATGATGGTCCCGTTCGACAGCATCAAATTCTCTGGCAACTACGGCAACATGACCGAAGTCTCCTATCAGGTTGCGAAACGTGCCGCGAAGAAAGGTGCTAAGTATTACCACATCACCCGCCAGTGGCAGGAACGTGGTAATAACCTGACCGTCAGCGCAGATCTGTATAAATAACTATTCAGTTACTCCTATAAGGCGGCTTGATGCCGCCTTTCGCATTTTTAGTCACAAATTTCCCCGCGACCATTGCATGC  >  NC_000913/1679157‑1679451

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