Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,708,652	G→C	9.1%	A578A (GCG→GCC)	rsxC →	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,708,652	0	G	C	9.1%	157.5 / 2.3	69	A578A (GCG→GCC)	rsxC	SoxR iron‑sulfur cluster reduction factor component; putative membrane‑associated NADH oxidoreductase of electron transport complex
Reads supporting (aligned to +/- strand):  ref base G (19/43);  new base C (3/4);  total (22/47)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTG  >  NC_000913/1708503‑1708800                                                                                                                                                      |                                                                                                                                                     gAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                                                                                                                      <  2:1182784/150‑1 (MQ=255)     cccGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGCCGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCCCGCAAGCTGGAACAGCAACAGGCTAATGCGGacc                                                                                                                                                  >  1:1603249/1‑148 (MQ=255)        gTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCGAGCTGGAACAGCAACAGGCTAATGCGGAACCAg                                                                                                                                               <  2:1582371/150‑1 (MQ=255)         tAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGCACCAGa                                                                                                                                              >  1:1492609/1‑150 (MQ=255)          aaaGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATAGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGAGAATGCGGAGCCAGaa                                                                                                                                             >  2:490375/1‑150 (MQ=255)                    ccAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAAAGCGGAACCAGAACAACAGGTCg                                                                                                                                   >  2:877339/1‑150 (MQ=255)                    ccAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGCTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTCGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGtgg                                                                                                                                   >  2:224414/1‑148 (MQ=255)                      agagCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCCTAACCAGCACAACCGGtacct                                                                                                                                 >  1:1004540/1‑146 (MQ=37)                           gAAACAGGCAGAACTGCAGCAAACTAACGACTCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcgcg                                                                                                                            <  2:341027/150‑1 (MQ=255)                           gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcgcg                                                                                                                            <  1:851245/150‑1 (MQ=255)                           gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcgcg                                                                                                                            <  2:1750081/150‑1 (MQ=255)                           gAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCcgcg                                                                                                                            <  2:341026/150‑1 (MQ=255)                            aaaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCgcgc                                                                                                                           >  1:627507/1‑150 (MQ=255)                            aaaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCgcgc                                                                                                                           <  1:1205297/150‑1 (MQ=255)                            aaaCAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCgcgc                                                                                                                           <  2:705041/150‑1 (MQ=255)                                 ggCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAgc                                                                                                                      <  2:542956/150‑1 (MQ=255)                                 ggCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTGGAAGCAGCTATCGCCCGCGCCAAAGAGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAgc                                                                                                                      >  1:1014277/1‑150 (MQ=255)                                   cAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATACACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTAGATCCACGAAAAcccc                                                                                                                    >  2:1729258/1‑149 (MQ=25)                                      aaCTTCAGCAAACTCACGACCCAGCAACAGAGGCTGATCCTCGAAAACCAGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGccgccg                                                                                                                 <  2:654070/150‑1 (MQ=17)                                             gcaAACTTACGAAGCAGCAACAGAGGCTGATCCTCGAAAACCAGCCGTAGAAGCAGCTATTGCCCGCGCCAAAGCACGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc                                                                                                          <  1:1186641/150‑1 (MQ=11)                                             gcaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTCGAACAGAAACAGGCTAATTCCGAGCCAGAACAAGAGGTCGATCCGCGAAAAGCCTACGTAGAAga                                                                                                          >  2:570116/1‑149 (MQ=255)                                             ccaAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATAGCCCGCGCCAAACCGCGCAAGCGTGAACAGAAACAGGCTAATGCCGAAACATAACAAAAGGTCGATCCGCGCAAAGCAGCCGTCGAAGc                                                                                                          >  1:1039726/2‑150 (MQ=255)                                              caAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                         >  2:1058046/1‑150 (MQ=255)                                              caAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                         <  1:1086657/150‑1 (MQ=255)                                                    cacgaccCAGCACCAGATGCCGATCCCCGTAAACCTGCCGTAGAAGCAGCTATTGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTAtt                                                                                                   <  1:1344829/149‑1 (MQ=14)                                                      cgacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAGCAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGc                                                                                                 >  1:1677135/1‑150 (MQ=255)                                                       gacgCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                                                >  2:925711/1‑150 (MQ=255)                                                       gaatcagcaCACGTTGTTGATCCTCGAAAACCAGCCGTCGAAGCCGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGAGGAACCAGAACAACAGGGCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGcc                                                                                                <  2:1677135/146‑1 (MQ=11)                                                           cagcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTg                                                                                            <  1:925711/150‑1 (MQ=255)                                                             gcaACCGTTGCTGATCCACGTAAAACTTCCGTTGAAGCAGCTATCGCCCGCTCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                                          <  1:1727466/150‑1 (MQ=255)                                                             gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGCAGCCGCTATTGCCCGAGcc                                                                                          >  1:856559/1‑150 (MQ=255)                                                             gcaACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGcc                                                                                          <  2:1339432/150‑1 (MQ=255)                                                                    ttgttgATCCACGCAAAGCAGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCTCGCAAGCTGGAACAGCAACAGGCTAATGCCGAGCCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAgcgc                                                                                   <  1:570116/150‑1 (MQ=11)                                                                     tGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcg                                                                                  <  1:431707/150‑1 (MQ=255)                                                                      gCTGATCCACGTATAACTGCAGTTGAAGCAGCTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                 <  2:283230/150‑1 (MQ=255)                                                                      gCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAgcgcgc                                                                                 <  2:38767/150‑1 (MQ=255)                                                                           tCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACCGGCTAGTGCGGAACCGGACCCGGCGGACGAACCGCGCCAGGCCGCCGTCGAAGCCGCTCTTGCCCGTGCCAAAGCACGCAGGcg                                                                            >  1:17076/1‑149 (MQ=17)                                                                            ccACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTg                                                                           <  2:748071/150‑1 (MQ=255)                                                                             cACTTAAAACTGCCGTTGAAGCAGACATCGCCCGCGCCAAAGCGCGCAAGCGGGAACAGCAACAGGCTAATGCGGACCCAGAACAACAGGTCGATCCGCGCAAAGCGGCCGGCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTgg                                                                          <  2:1712250/150‑1 (MQ=17)                                                                              cctttaaaCTGCCCTTCTACCAGCTAACGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                         <  1:611573/145‑1 (MQ=14)                                                                              ccGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                         <  1:771066/149‑1 (MQ=255)                                                                              aCGTAAAACTGCAGTTGAAGCAGCTATCTACCGCGCAAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAGCAGAACAAGAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGa                                                                         <  2:685049/150‑1 (MQ=17)                                                                                gTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                       <  1:590321/150‑1 (MQ=255)                                                                                ctAAAACTGCCGTTGAAGCAGCTATCGACCGCTCCAAAGCGCGCAAGCTGGAACAGCAACAGGATAATGCGGAACCAGAACAACATGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAAc                                                                       <  2:822698/149‑1 (MQ=255)                                                                                   aaaCTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGc                                                                    <  2:461868/150‑1 (MQ=255)                                                                                    taCTGCCGGAGTAGCAGCTATGTCCCGCGCCAACTCTCGCAATCTGGAACATCAACAGGAAACTCCGGAACCAGCACACCAGGCCGATCCGCGCAAACCCGCCGTCGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCa                                                                   <  2:1151136/149‑1 (MQ=19)                                                                                    aaCTGCCGTAGAAGCAGCTATTTCCCGCGCCTAAGCGCGCAAGCTGGAACACCAACAGGCTAATCAGGAACCAGAACCCCAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCACGCAAGCTGGAACAGCa                                                                   <  2:40444/150‑1 (MQ=11)                                                                                                 gcagcTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCg                                                      <  1:639889/150‑1 (MQ=255)                                                                                                     cTATCGCCCGCGCCAACGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                  <  1:1500528/150‑1 (MQ=255)                                                                                                     cTATCGCCCGCGCCAAAGCGCGGAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAAc                                                  <  1:1814140/150‑1 (MQ=255)                                                                                                       aTCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCGCAGGCGCGCAAGCTGGAACAGCACCCGGCTAATGCGGAACCa                                                >  1:632185/1‑150 (MQ=14)                                                                                                               cgcCAAATCGCGCAATCTGGAACAGCATCAGGCTAATGCGGAACCAGAACACCTGGTCGATCCGCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAAcg                                        <  1:1799872/150‑2 (MQ=14)                                                                                                               cgcCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACa                                        <  1:1058046/150‑1 (MQ=255)                                                                                                                   aaatggcgcAATCTGGAAGATCTACAGGATAATGACGAACCAGAAAAACAGGTCGATCCTCGCAAACCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                    <  1:1037723/145‑1 (MQ=14)                                                                                                                   aaaGCTCGCAAGCTGGAACATCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAATCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTc                                    <  1:1802796/150‑1 (MQ=255)                                                                                                                   aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGATCCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                                                      <  2:1190516/84‑1 (MQ=18)                                                                                                                   aaaGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCt                                                                                                      >  1:1190516/1‑84 (MQ=33)                                                                                                                      tgccgcAAGCTGGATCATCAACAGGATAATTAGTAACCATTACAAATTTTAGATCCTCGTAAATCCTCCTTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACACCAACTGGCTAATGCGGAACCAGAACAACAGGTCGAt                                 <  2:1161740/147‑1 (MQ=2)                                                                                                                      gctcgcAATCTGGAACATCAACAGGATAATTCGGAACCAGCACAACAGGTCGATCCGCGCAAATCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGAt                                 <  1:911692/147‑1 (MQ=14)                                                                                                                       ggcggAATCTGGAACATCAACAGGAAAATTAGGAACCAGAACAACAGGTCGATCCTCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                <  1:388762/149‑1 (MQ=14)                                                                                                                       ggcgcAATCTGGAACATCATCAGGCAAATTCGGAACCAGAAACACAGGTAGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATc                                <  1:395639/149‑1 (MQ=14)                                                                                                                          gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                         <  2:86523/74‑1 (MQ=34)                                                                                                                          gcAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGcc                                                                                                         >  1:86523/1‑74 (MQ=34)                                                                                                                                 ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                        <  1:562783/84‑1 (MQ=18)                                                                                                                                 ggAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCaa                                                                                        >  2:562783/1‑84 (MQ=18)                                                                                                                                             ggCTAATTCGGAACCAGAACCACAGGCAGATCCTCGAAAATCTGCCGTAGAAGCAGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc          <  2:526279/150‑1 (MQ=14)                                                                                                                                             ggCTAATTCGGAACCAGAACAAGAGGTCGATCCGCGAAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGc          <  2:1460507/150‑1 (MQ=255)                                                                                                                                             ggCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGAGCCAAAGCGCGCaaac                                                                             >  2:1719390/1‑81 (MQ=9)                                                                                                                                                     cGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCgag                                                                                            <  2:1083356/60‑3 (MQ=25)                                                                                                                                                     cGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTAGAAGCCGCTATTGCCCGAGCCAAAGCGCGCAAACTCGAACAGCAACAGGCTAATGCCGAACCAGAAGAACAGGTCGATCCGCGCAAAGCAGCCGTTGAAGCGGCTATTg  >  2:578448/1‑150 (MQ=11)                                                                                                                                                      |                                                                                                                                                     GAAGCCCGTAAAGCGCAAGCCAGAGCGAAACAGGCAGAACTGCAGCAAACTAACGACGCAGCAACCGTTGCTGATCCACGTAAAACTGCCGTTGAAGCAGCTATCGCCCGCGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAACAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTGCCCGTGCCAAAGCGCGCAAGCTGGAACAGCAACAGGCTAATGCGGAACCAGAAGAACAGGTCGATCCGCGCAAAGCCGCCGTCGAAGCCGCTATTG  >  NC_000913/1708503‑1708800

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