Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,381,596	A→G	G241G (GGT→GGC)	frdA ←	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,381,596	0	A	G	98.0%	365.0 / ‑1.9	98	G241G (GGT→GGC)	frdA	anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (46/50);  total (48/50)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

TTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCA  >  NC_000913/4381450‑4381740                                                                                                                                                   |                                                                                                                                                 ttGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGcc                                                                                                                                                <  1:202409/149‑1 (MQ=255) ttGTCGCG‑GGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGcca                                                                                                                                               <  1:702146/149‑2 (MQ=255)    tCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCGGCGGAGTTTCCGGCCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGaa                                                                                                                                            <  1:339608/150‑1 (MQ=255)    tCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAAc                                                                                                                                           >  2:1134912/1‑151 (MQ=255)      gcgTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCt                                                                                                                                         >  1:843857/1‑151 (MQ=255)        gTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCt                                                                                                                                         <  1:499967/149‑1 (MQ=255)        gTGGACCCAGTTCCATATATTTGTTTTGCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTgg                                                                                                                                       <  2:258949/151‑1 (MQ=255)         tGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGc                                                                                                                                      >  2:469924/1‑151 (MQ=255)         tGGACCCAGTTCCATATATTTGTTTTACGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTTACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTgg                                                                                                                                       <  2:1484126/150‑1 (MQ=255)          ggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGc                                                                                                                                      <  2:1569877/150‑1 (MQ=255)          ggACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGc                                                                                                                                      >  2:809167/1‑150 (MQ=255)           gACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGc                                                                                                                                      >  1:666729/1‑149 (MQ=255)             cccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGAc                                                                                                                                  <  1:906936/151‑1 (MQ=255)              ccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGAcc                                                                                                                                 <  1:1180964/151‑1 (MQ=255)              ccAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGAcc                                                                                                                                 >  2:1103235/1‑151 (MQ=255)               cAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCg                                                                                                                                <  1:1682044/151‑1 (MQ=255)                   tCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGtt                                                                                                                             <  2:579072/150‑1 (MQ=255)                   tCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTg                                                                                                                            >  1:691545/1‑151 (MQ=255)                   tCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTg                                                                                                                            >  2:642157/1‑151 (MQ=255)                    ccATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTgg                                                                                                                           >  1:1452739/1‑151 (MQ=255)                    ccATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTgg                                                                                                                           <  1:365568/151‑1 (MQ=255)                         tatTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGAt                                                                                                                      >  1:587539/1‑151 (MQ=255)                          atTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGAt                                                                                                                      <  2:370589/150‑1 (MQ=255)                           tttGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATa                                                                                                                     <  1:835266/150‑1 (MQ=255)                            ttGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATAc                                                                                                                    <  2:1067011/150‑1 (MQ=255)                                 tttCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAAc                                                                                                               <  1:1027339/150‑1 (MQ=255)                                 tttCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACg                                                                                                              >  2:341479/1‑151 (MQ=255)                                    cGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAt                                                                                                           <  1:388800/151‑1 (MQ=255)                                     ggCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAAtt                                                                                                          <  1:218782/151‑1 (MQ=255)                                          gCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg                                                                                                     >  1:703061/1‑151 (MQ=255)                                          gCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg                                                                                                     >  1:1008103/1‑151 (MQ=255)                                          gCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg                                                                                                     >  2:1020054/1‑151 (MQ=255)                                          gCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATg                                                                                                     >  2:601335/1‑151 (MQ=255)                                            ccAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTc                                                                                                   >  2:925789/1‑151 (MQ=255)                                            ccAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTc                                                                                                   <  2:1261711/151‑1 (MQ=255)                                            ccAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGCAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATAAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTc                                                                                                   >  1:767526/1‑151 (MQ=255)                                             cAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTc                                                                                                   >  1:629358/1‑150 (MQ=255)                                                 ggAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGc                                                                                               <  1:1496436/150‑1 (MQ=255)                                                  gAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGATGGGTGATACTGAACGAATTCCATGTCACGc                                                                                               >  1:1084405/1‑149 (MQ=255)                                                    gTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGc                                                                                            <  1:963971/150‑1 (MQ=255)                                                    gTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCg                                                                                           >  2:1593933/1‑151 (MQ=255)                                                       tCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGa                                                                                         >  1:1405163/1‑150 (MQ=255)                                                        ccGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGaa                                                                                        <  1:873592/150‑1 (MQ=255)                                                        ccGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGa                                                                                         <  2:234753/149‑1 (MQ=255)                                                        ccGGCCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGa                                                                                         <  2:890754/149‑1 (MQ=255)                                                         cGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACg                                                                                      <  2:976263/151‑1 (MQ=255)                                                           ggCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGcc                                                                                    <  2:407718/151‑1 (MQ=255)                                                            gCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCg                                                                                   >  2:888859/1‑151 (MQ=255)                                                            gCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCg                                                                                   >  2:842364/1‑151 (MQ=255)                                                               cATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTgg                                                                                >  1:970871/1‑151 (MQ=255)                                                                aTGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGc                                                                               >  1:778503/1‑151 (MQ=255)                                                                aTGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGc                                                                               >  2:159703/1‑151 (MQ=255)                                                                 tGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCt                                                                              >  2:739068/1‑151 (MQ=255)                                                                  gCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCt                                                                              <  1:711331/150‑1 (MQ=255)                                                                   ccGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCtt                                                                             >  1:176693/1‑150 (MQ=255)                                                                            tGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGccat                                                                     <  1:1117910/149‑1 (MQ=255)                                                                            tGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGcc                                                                       <  2:1066619/147‑1 (MQ=255)                                                                             gCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGccatacc                                                                  <  2:581138/151‑1 (MQ=255)                                                                              cAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGccataccc                                                                 <  2:850564/151‑1 (MQ=255)                                                                                  tAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGccatacctatac                                                             >  1:771293/1‑151 (MQ=255)                                                                                   aaCGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCcatacccatacc                                                            >  1:713789/1‑151 (MQ=255)                                                                                     cGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGt                                                          <  1:1124174/151‑1 (MQ=255)                                                                                       gTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCa                                                        <  2:136294/151‑1 (MQ=255)                                                                                        tAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCAc                                                       >  1:1377556/1‑151 (MQ=255)                                                                                         aGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCAcc                                                      <  1:339998/151‑1 (MQ=255)                                                                                         aGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCAcc                                                      >  2:1070036/1‑151 (MQ=255)                                                                                         aGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCAcc                                                      >  2:1611541/1‑151 (MQ=255)                                                                                         aGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCAcc                                                      >  2:375339/1‑151 (MQ=255)                                                                                           ccATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCgg                                                    <  1:1028648/151‑1 (MQ=255)                                                                                            cATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGt                                                   <  2:1591286/151‑1 (MQ=255)                                                                                            cATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGt                                                   <  1:1512088/151‑1 (MQ=255)                                                                                            cATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGt                                                   >  1:1307848/1‑151 (MQ=255)                                                                                                tttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAc                                                <  2:176704/150‑1 (MQ=255)                                                                                                tttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAAc                                                <  1:220251/150‑1 (MQ=255)                                                                                                tttGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACg                                               >  2:1278434/1‑151 (MQ=255)                                                                                                 ttgttgACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACg                                               <  2:756149/150‑1 (MQ=255)                                                                                                   gttgACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGAt                                             <  2:741699/150‑1 (MQ=255)                                                                                                   gttgACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGAt                                             <  1:1567262/150‑1 (MQ=255)                                                                                                       aCCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATgccgc                                        >  2:860063/1‑151 (MQ=255)                                                                                                        ccAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATgccgcc                                       >  2:1259619/1‑151 (MQ=255)                                                                                                          aGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGt                                     <  2:1452837/151‑1 (MQ=255)                                                                                                           gAATCCCGCTTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGAAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATCCCCATACCGTCACCGGTAACGATGCCGCCGtt                                    <  2:666788/151‑1 (MQ=255)                                                                                                               aCCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGgtg                                >  1:1336599/1‑151 (MQ=255)                                                                                                                ccGCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGgtgt                               >  1:1706543/1‑151 (MQ=255)                                                                                                                  gCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGtgttg                             >  2:850933/1‑151 (MQ=255)                                                                                                                  gCCTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGtgttg                             >  1:790294/1‑151 (MQ=255)                                                                                                                   ccTTCACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGtgttgt                            >  1:45187/1‑151 (MQ=255)                                                                                                                       cACCGCGGCAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAAc                         <  1:127340/150‑1 (MQ=255)                                                                                                                          cgcgGCAACCTTCGGTCATCATGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTCGGTTAGCGCCATACCCATACCGTCCCCGGTAACGATGCCGCCGTTGGTGTTGTAACGATa                     <  2:278956/151‑1 (MQ=255)                                                                                                                             ggCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATaaa                   >  2:1151742/1‑150 (MQ=255)                                                                                                                               cAACCTTCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAAcgc                >  2:1445948/1‑151 (MQ=255)                                                                                                                                    ttCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGAccc           >  2:58868/1‑151 (MQ=255)                                                                                                                                    ttCGGTCATCAGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGAccc           >  2:1610877/1‑151 (MQ=255)                                                                                                                                          catcaGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCAcc      <  1:159695/150‑1 (MQ=255)                                                                                                                                          catcaGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCAcc      <  1:375311/150‑1 (MQ=255)                                                                                                                                           atcaGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGc    <  2:1147781/151‑1 (MQ=255)                                                                                                                                             caGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCa  <  1:1134824/151‑1 (MQ=255)                                                                                                                                              aGGATGCCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGcc   >  2:1325821/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                 TTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTAGCCATTTTTGTTGACCAGAATACCGCCTTCACCGCGGCAACCTTCGGTCATCAGGATACCGGAACCTGGCAGACCGGTTGGGTGATACTGAACGAATTCCATGTCACGCAGCGGAACGCCGTGGCTTAGCGCCATACCCATACCGTCACCGGTAACGATGCCGCCGTTGGTGTTGTAACGATAAACGCGACCCGCACCGCCA  >  NC_000913/4381450‑4381740

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