Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,117,687	A→G	6.2%	F135S (TTC→TCC)	glpF ←	glycerol facilitator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,117,687	0	A	G	6.2%	197.7 / 3.7	65	F135S (TTC→TCC)	glpF	glycerol facilitator
Reads supporting (aligned to +/- strand):  ref base A (25/36);  new base G (2/2);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

CAATCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATAAAAGGAATAA  >  NC_000913/4117541‑4117834                                                                                                                                                   |                                                                                                                                                    cAATCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACGGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGaaa                                                                                                                                                   <  4:384126/149‑1 (MQ=255) cAATCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGaaa                                                                                                                                                   >  2:85947/1‑149 (MQ=255) cAATCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGaaa                                                                                                                                                   >  2:258397/1‑149 (MQ=255) cAATCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGaaa                                                                                                                                                   >  4:314320/1‑149 (MQ=255)  aaTCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAg                                                                                                                                                  <  1:258397/149‑1 (MQ=255)   aTCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGt                                                                                                                                                 >  2:100190/1‑149 (MQ=255)      aGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcc                                                                                                                                              <  3:252741/149‑1 (MQ=255)      aGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcc                                                                                                                                              <  3:433584/149‑1 (MQ=255)       gCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcca                                                                                                                                             <  3:245815/149‑1 (MQ=255)       gCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTgcca                                                                                                                                             >  1:395054/1‑149 (MQ=255)             gAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGa                                                                                                                                       <  2:42283/149‑1 (MQ=255)             gAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGACAGGATCAACCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGa                                                                                                                                       <  3:489323/149‑1 (MQ=255)              aGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGAt                                                                                                                                      >  1:186063/1‑149 (MQ=255)              aGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGAt                                                                                                                                      <  1:110170/149‑1 (MQ=255)                ccAAAGGGCCGTGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGGAAGTgc                                                                                                                                               >  3:330787/1‑138 (MQ=255)                ccAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGGAAGTgc                                                                                                                                               >  1:264096/1‑138 (MQ=255)                ccAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGGAAGTgc                                                                                                                                               <  4:330787/138‑1 (MQ=255)                ccAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGGAAGTgc                                                                                                                                               <  2:264096/138‑1 (MQ=255)                      ggCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACtt                                                                                                                              <  1:295057/149‑1 (MQ=255)                       gCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACttt                                                                                                                             >  4:68514/1‑149 (MQ=255)                        ccGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTc                                                                                                                            <  2:32937/149‑1 (MQ=255)                        acgcgAGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTc                                                                                                                            >  2:258383/2‑149 (MQ=255)                             tGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGc                                                                                                                       <  3:539547/149‑1 (MQ=255)                             tGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGACCAACACTTTCAACGc                                                                                                                       >  1:264685/1‑149 (MQ=255)                                tACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGc                                                                                                                    >  4:178509/1‑149 (MQ=255)                                tACACCGTTGCCATCGTCCGTTAACACCAGGATGAGTCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATACGAGGTTTAGGGTAAGGAGAGAAAGTGACAGCCAGATCAACACTTTCAACGCTGc                                                                                                                    >  2:526630/1‑149 (MQ=255)                                      gTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGaa                                                                                                              <  3:497169/149‑1 (MQ=255)                                      gTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGaa                                                                                                              >  2:492470/1‑149 (MQ=255)                                          ccATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAAt                                                                                                          <  1:109682/149‑1 (MQ=255)                                          ccATCGTCCGTTAACGCCACGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAAt                                                                                                          <  1:569815/149‑1 (MQ=255)                                           cATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAAtg                                                                                                         <  1:85947/149‑1 (MQ=255)                                            aTCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAAtgt                                                                                                        <  1:452607/149‑1 (MQ=255)                                              cGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGtga                                                                                                      <  1:488607/149‑1 (MQ=255)                                              cGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGtga                                                                                                      <  2:437547/149‑1 (MQ=255)                                              cGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGtga                                                                                                      <  3:405397/149‑1 (MQ=255)                                              cGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGtga                                                                                                      <  3:260758/149‑1 (MQ=255)                                                 ccGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGtgatga                                                                                                   <  2:160144/149‑1 (MQ=255)                                                              gATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTc                                                                                                >  1:326670/1‑139 (MQ=255)                                                              gATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTc                                                                                                <  2:326670/139‑1 (MQ=255)                                                                    ccccATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGaaaa                                                                                >  2:260489/1‑149 (MQ=255)                                                                      ccATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGaa                                                                                        <  2:287830/139‑1 (MQ=255)                                                                      ccATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGaa                                                                                        >  1:4093/1‑139 (MQ=255)                                                                      ccATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGaa                                                                                        >  1:287830/1‑139 (MQ=255)                                                                      ccATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGaa                                                                                        <  2:4093/139‑1 (MQ=255)                                                                      ccATCAGAATAGCGGTAATCACCATCTCAACAGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAGAGTGCCAACCAGATCAACACTTTCAACGCTGCCGCGAACAGTGTGATGCGTCTGCTCGAAGTCGAAaaat                                                                              >  3:349419/1‑149 (MQ=255)                                                                              atcgCGGTAATCACTATCTCAACTGCGAAAGCCTACACAAAATTGATATGAGGATTAGGGTTAGTAGAGAAAGTGCCAGCCAGATCAACAGTTTCAACGCTGCCGAGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgta                                                                      <  4:349419/146‑1 (MQ=255)                                                                               tAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTgtag                                                                     <  1:260489/149‑1 (MQ=255)                                                                                       aTCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGccc                                                             <  3:502061/149‑1 (MQ=255)                                                                                          aCCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTa                                                          >  4:453068/1‑149 (MQ=255)                                                                                             aTCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAAc                                                       >  1:549197/1‑149 (MQ=255)                                                                                             aTCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAAc                                                       >  1:180150/1‑149 (MQ=255)                                                                                                    cTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGcc                                                <  3:424262/149‑1 (MQ=255)                                                                                                     tGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCg                                               >  3:544371/1‑149 (MQ=255)                                                                                                        gAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCgcag                                            >  2:554157/1‑149 (MQ=255)                                                                                                                 cacaAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAg                                   >  2:23646/1‑149 (MQ=255)                                                                                                                  acaAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAgc                                  <  2:384921/149‑1 (MQ=255)                                                                                                                      aaTTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCg                              >  2:110674/1‑149 (MQ=255)                                                                                                                        ttGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGc                            >  4:336756/1‑149 (MQ=255)                                                                                                                           atatGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAAc                         <  3:178509/149‑1 (MQ=255)                                                                                                                           atatGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCAAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTTAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAAc                         <  3:452290/149‑1 (MQ=255)                                                                                                                           atatGAGGATTAGGGTAAGAAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAAc                         <  1:526630/149‑1 (MQ=255)                                                                                                                               gAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAAtgtg                                                                                                       <  3:328573/67‑1 (MQ=255)                                                                                                                               gAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAAtgtg                                                                                                       >  4:328573/1‑67 (MQ=255)                                                                                                                                  gATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGaa                  >  1:103051/1‑149 (MQ=255)                                                                                                                                   aTTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGaaa                 >  4:12117/1‑149 (MQ=255)                                                                                                                                    ttAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAAc                <  3:12117/149‑1 (MQ=255)                                                                                                                                           aaGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATaaaa         <  2:553877/149‑1 (MQ=255)                                                                                                                                           aaGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATaaaa         <  2:523665/149‑1 (MQ=255)                                                                                                                                             gTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATAAAAgg       <  3:515684/149‑1 (MQ=255)                                                                                                                                                gagaAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATAAAAGGAAt    >  3:118634/1‑149 (MQ=255)                                                                                                                                                gagaAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATAAAAGGAAt    <  2:438291/149‑1 (MQ=255)                                                                                                                                                  gaAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATAAAAGGACTaa  >  4:517974/1‑149 (MQ=255)                                                                                                                                                   aaaGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAAc                <  2:279078/134‑1 (MQ=255)                                                                                                                                                   aaaGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAAc                >  1:279078/1‑134 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAATCAGCAAGGGAGCCAAAGGGCCGCGTGGTACACCGTTGCCATCGTCCGTTAACGCCAGGATCAGCCCCATCAGAATAGCGGTAATCACCATCTCAACTGCGAAAGCCTGCACAAAATTGATATGAGGATTAGGGTAAGTAGAGAAAGTGCCAGCCAGATCAACACTTTCAACGCTGCCGCGAACAATGTGATGAGTCTGCTCGAAGTCGAAAAATAAATTGTAGTAAAGCCCGTAAACTAAAGCCGCAGCACAGAAAGCGCCGGCAACTTGTGAAACGATAAAAGGAATAA  >  NC_000913/4117541‑4117834

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