Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	937,158	T→A	10.2%	intergenic (+65/‑1952)	yumB → / → USA300HOU_0903	NADH dehydrogenase/NhaC family sodium:proton (Na+:H+) antiporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	937,158	0	T	A	10.2%	199.6 / 8.9	66	intergenic (+65/‑1952)	yumB/USA300HOU_0903	NADH dehydrogenase/NhaC family sodium:proton (Na+:H+) antiporter
Reads supporting (aligned to +/- strand):  ref base T (23/36);  new base A (2/5);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAAAAGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAACATTTAAAT  >  CP000730/937019‑937298                                                                                                                                            |                                                                                                                                             aaaaaaGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                                                                             <  2:551293/141‑1 (MQ=255) aaaaaaGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTc                                                                                                                                             <  2:551310/141‑1 (MQ=255)   aaaaGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCAc                                                                                                                                           >  1:132020/1‑141 (MQ=255)    aaaGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACt                                                                                                                                          <  1:216841/141‑1 (MQ=255)     aaGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTg                                                                                                                                         >  2:494425/1‑141 (MQ=255)      aGTGATTTATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGACTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGa                                                                                                                                        <  2:260679/141‑1 (MQ=255)        tgattgatACACGTGCGGTATTCAAAATCGGTGGTATTGGTTTAGAATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAg                                                                                                                                      <  2:697999/141‑1 (MQ=255)         gattgatACACGTGCGGTATTCAAAATCGGTTGTAGCTGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgg                                                                                                                                     <  2:571070/141‑1 (MQ=255)          attgatACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAggg                                                                                                                                    <  1:186217/141‑1 (MQ=255)          attgatACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGAATTCAAAAACGGTCCATCCACGGCTAAGAGTAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAgga                                                                                                                                    <  1:136837/141‑2 (MQ=255)           ttgatATACGTGATGTATTCAAAATCGGTAGTATATATTTAATATTAAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTTATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                                                                   <  2:356293/141‑1 (MQ=255)           ttgatACCCGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGc                                                                                                                                   <  2:43909/141‑1 (MQ=255)           ttgatACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGCCCTTCAGTGACACTGAAggac                                                                                                                                   <  1:289412/141‑3 (MQ=255)                    gTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                                                                                          <  1:91687/141‑1 (MQ=255)                    gTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTAttgtt                                                                                                                          <  1:217934/141‑1 (MQ=255)                      gCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATtgttgt                                                                                                                        <  1:60459/141‑1 (MQ=255)                      gCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATtgttgt                                                                                                                        <  1:560887/141‑1 (MQ=255)                       cGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTa                                                                                                                       <  2:169983/141‑1 (MQ=255)                       cGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTa                                                                                                                       >  1:426159/1‑141 (MQ=255)                         gTATTCAAAATCGGTGGTATCGGTTTAGCATTCATAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                                                                     <  1:117288/141‑1 (MQ=255)                         gTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACa                                                                                                                     <  2:657531/141‑1 (MQ=255)                           aTTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                                                                   >  1:246532/1‑141 (MQ=255)                           aTTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAAt                                                                                                                   <  2:617484/141‑1 (MQ=255)                               aaaaTCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGa                                                                                                               <  2:426159/141‑1 (MQ=255)                                 aaTCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAg                                                                                                             >  1:366377/1‑141 (MQ=255)                                 aaTCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTAAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAg                                                                                                             >  1:363888/1‑141 (MQ=255)                                  aTCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCa                                                                                                                                            >  2:117254/1‑109 (MQ=255)                                  aTCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCa                                                                                                                                            <  1:117254/109‑1 (MQ=255)                                   tCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTg                                                                                                           >  1:453600/1‑141 (MQ=255)                                     ggtggtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAg                                                                                                         <  2:575093/141‑1 (MQ=255)                                      gtggtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGt                                                                                                        >  2:562778/1‑141 (MQ=255)                                         gtATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAAc                                                                                                     <  1:567775/141‑1 (MQ=255)                                          tATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACa                                                                                                    >  1:310279/1‑141 (MQ=255)                                              ggTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATc                                                                                                <  2:665662/141‑1 (MQ=255)                                                 ttAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCaaa                                                                                             <  1:359045/141‑1 (MQ=255)                                                  tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                                                                            >  1:102755/1‑141 (MQ=255)                                                  tAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAAt                                                                                            >  1:102914/1‑141 (MQ=255)                                                        tCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAagg                                                                                      <  1:540029/141‑1 (MQ=255)                                                                 gTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACAtt                                                                             >  2:540019/1‑141 (MQ=255)                                                                   aaaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTAt                                                                           >  1:659564/1‑141 (MQ=255)                                                                    aaTTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATg                                                                          <  2:125147/141‑1 (MQ=255)                                                                          aGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAAtttt                                                                    >  2:461801/1‑141 (MQ=255)                                                                                taaataaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTttaaatta                                                              <  2:208149/141‑1 (MQ=255)                                                                                    taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa                                                          <  1:152687/141‑1 (MQ=255)                                                                                    taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATa                                                          <  1:152666/141‑1 (MQ=255)                                                                                    taaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTCAAGGGCCTTATTGTTGTACAATTTGAAGCGAGTAACATATCAAATCTAAGCAGGTACATTATGAATTTTAAATTAAATa                                                          <  2:235329/141‑1 (MQ=255)                                                                                     aaaGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATaa                                                         >  2:579989/1‑141 (MQ=255)                                                                                        gTTTCAGCTAAACACATAACATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTAATTGATGTACACTTTGAAGTGAGTACCCTTTCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATGacac                                                      <  2:105516/141‑1 (MQ=255)                                                                                        gTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAacac                                                      <  1:182015/141‑1 (MQ=255)                                                                                              gCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAg                                                >  2:201766/1‑141 (MQ=255)                                                                                                  aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAAc                                            <  1:562778/141‑1 (MQ=255)                                                                                                  aaCACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAAc                                            <  2:468547/141‑1 (MQ=255)                                                                                                          aaattaaattATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTACATTAAATAACACACTAAGCAACGaaataaa                                    <  1:76216/141‑1 (MQ=255)                                                                                                                 attattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTg                             >  2:7104/1‑141 (MQ=255)                                                                                                                  ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                                                                            <  1:205809/93‑1 (MQ=255)                                                                                                                  ttattCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTa                                                                            >  2:205809/1‑93 (MQ=255)                                                                                                                       cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGattatt                       <  2:310279/141‑1 (MQ=255)                                                                                                                       cATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGattatt                       >  1:417289/1‑141 (MQ=255)                                                                                                                             taAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGtata                 <  1:76059/141‑1 (MQ=255)                                                                                                                               aaGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTATATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATAcc               <  1:579989/141‑1 (MQ=255)                                                                                                                                 gcccTTCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAGGTGAGTAACATATCAAATCTATGGAGGTACATTAAGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAg             >  2:323645/3‑141 (MQ=255)                                                                                                                                 gTCGTTCAGTGTCACTGAAGGGCGTTATTGTTGTCCAATTTGAGGTGTGTAACATATCAAATCTAAGGAGGGACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAg             >  2:64503/1‑141 (MQ=255)                                                                                                                                 gTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAg             >  1:279663/1‑141 (MQ=255)                                                                                                                                   ccTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGaa           >  2:408439/1‑141 (MQ=255)                                                                                                                                    cTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAAc          <  2:132020/141‑1 (MQ=255)                                                                                                                                      tCAGTGACACTGAAGGACTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAACAt        >  1:467848/1‑141 (MQ=255)                                                                                                                                         gtgtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAACATTTa     >  2:202477/1‑141 (MQ=255)                                                                                                                                           gtCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAACATTTaaa   >  2:375545/1‑141 (MQ=255)                                                                                                                                            tCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAACATTTAAAt  <  1:202477/141‑1 (MQ=255)                                                                                                                                            |                                                                                                                                             AAAAAAGTGATTGATACACGTGCGGTATTCAAAATCGGTGGTATCGGTTTAGCATTCAAAAAAGGTAAATTCTAGGCTTTAAATAAAGTTTCAGCTAAACACATAAAATTAAATTATTCATTTATAAAGTCCTTCAGTGTCACTGAAGGGCTTTATTGTTGTACAATTTGAAGTGAGTAACATATCAAATCTAAGGAGGTACATTATGAATTTTAAATTAAATAACACACTAAGCAACGAAATAAATACATTGATTATTGGTATACCAGAACATTTAAAT  >  CP000730/937019‑937298

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