Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,257,774	G→A	C198C (TGC→TGT)	psuG ←	pseudouridine 5'‑phosphate glycosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,257,774	0	G	A	98.9%	305.1 / ‑4.5	91	C198C (TGC→TGT)	psuG	pseudouridine 5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (40/50);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 4.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

GGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGAT  >  NC_000913/2257677‑2257870                                                                                                  |                                                                                                 ggCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACa                                                                                                <  2:1095900/100‑1 (MQ=255) ggCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACa                                                                                                <  1:2295821/100‑1 (MQ=255)  gCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAg                                                                                               <  2:1950831/100‑1 (MQ=255)  gCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAg                                                                                               >  1:2784119/1‑100 (MQ=255)       cGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTAAAACGGGCTGGTGCGACAGaaaaa                                                                                          <  2:1126077/100‑1 (MQ=255)        gTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAAc                                                                                         <  1:1167966/100‑1 (MQ=255)        gTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAAc                                                                                         >  1:2672811/1‑100 (MQ=255)        gTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAAc                                                                                         <  2:538612/100‑1 (MQ=255)         ttCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACg                                                                                        <  1:2064281/100‑1 (MQ=255)         ttCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACg                                                                                        <  2:2973557/100‑1 (MQ=255)           cAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAAAgca                                                                                      >  2:1312662/1‑100 (MQ=255)              aCCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACgcaggc                                                                                   <  1:892092/100‑1 (MQ=255)              aCCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACgcaggc                                                                                   >  1:450606/1‑100 (MQ=255)               ccGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACgcaggca                                                                                  >  1:1449914/1‑100 (MQ=255)               ccGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACgcaggca                                                                                  >  1:1545073/1‑100 (MQ=255)                     tGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCtt                                                                            <  2:2671764/100‑1 (MQ=255)                      gCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCttt                                                                           <  2:1147969/100‑1 (MQ=255)                      gCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCttt                                                                           <  1:2266464/100‑1 (MQ=255)                        cATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAg                                                                         <  2:588376/100‑1 (MQ=255)                        cATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAg                                                                         >  2:80177/1‑100 (MQ=255)                        cATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACACAGGCAGCGCTTTAg                                                                         <  2:2406167/100‑1 (MQ=255)                         aTTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGt                                                                        >  1:2549150/1‑100 (MQ=255)                            tCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTg                                                                     >  1:641140/1‑100 (MQ=255)                             caccacCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGa                                                                    >  2:949176/1‑100 (MQ=255)                              accaccATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGAt                                                                   <  1:2935535/100‑1 (MQ=255)                              accaccATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGAt                                                                   <  1:2087307/100‑1 (MQ=255)                               ccaccaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTCGTCTGATa                                                                  >  2:294880/1‑100 (MQ=255)                                caccaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAg                                                                 >  1:1629937/1‑100 (MQ=255)                                caccaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAg                                                                 <  2:1953152/100‑1 (MQ=255)                                caccaTTGCACGGGCAATTTCACTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAg                                                                 >  2:1998768/1‑100 (MQ=255)                                 accaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGc                                                                >  2:282948/1‑100 (MQ=255)                                  ccaTTGCACAGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGc                                                                >  1:403278/1‑99 (MQ=255)                                   caTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGcc                                                               <  2:38726/99‑1 (MQ=255)                                    aTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCaa                                                             <  1:286213/100‑1 (MQ=255)                                    aTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCaa                                                             <  2:1290141/100‑1 (MQ=255)                                     ttGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAAt                                                            <  2:530561/100‑1 (MQ=255)                                     ttGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAAt                                                            >  1:1243177/1‑100 (MQ=255)                                      tGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAAtt                                                           <  2:2825006/100‑1 (MQ=255)                                        naCGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTaa                                                         >  2:1906885/2‑100 (MQ=255)                                        cACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTaa                                                         >  1:2922717/1‑100 (MQ=255)                                          cGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACg                                                       >  1:2608245/1‑100 (MQ=255)                                          cGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACg                                                       <  1:171009/100‑1 (MQ=255)                                           gggCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACgg                                                      >  1:1166870/1‑100 (MQ=255)                                           gggCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACgg                                                      <  1:1382767/100‑1 (MQ=255)                                             gCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAACAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGc                                                     <  1:608257/99‑1 (MQ=255)                                                 tttCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACAcc                                                >  1:535305/1‑100 (MQ=255)                                                 tttCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACAcc                                                <  1:982726/100‑1 (MQ=255)                                                    cGCTGGCGCTGTCGAGACGAATGCTGACGACAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTACCGGCACACCGaa                                             >  2:77922/1‑100 (MQ=255)                                                    cGCTGGCGCTGTCGAGAAGAATGCTGACGTCAAACGGGATGCTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGaa                                             <  2:1008109/100‑1 (MQ=255)                                                     gCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAg                                            <  1:2268669/100‑1 (MQ=255)                                                     gCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAg                                            <  2:2797715/100‑1 (MQ=255)                                                      cTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAgg                                           >  1:2993612/1‑100 (MQ=255)                                                         gcgcTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGttt                                        >  2:533910/1‑100 (MQ=255)                                                          cgcTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTc                                       >  2:2703809/1‑100 (MQ=255)                                                              gTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTaaa                                   >  1:1350304/1‑100 (MQ=255)                                                                 gagaCGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATAc                                >  2:703688/1‑100 (MQ=255)                                                                 gagaCGAATGCTGACGTCAAACGGGCCGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATAc                                >  1:1841616/1‑100 (MQ=255)                                                                  agaCGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGAGCTTTAGTCGGATAGCCAATTAACGGCCCACCGAAGGTTTCTAAATACt                               >  2:1262637/1‑100 (MQ=255)                                                                   gaCGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAGATACTc                              <  2:207821/100‑1 (MQ=255)                                                                   gaCGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTc                              <  2:1545060/100‑1 (MQ=255)                                                                     cGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAg                            >  1:3011214/1‑100 (MQ=255)                                                                     cGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAg                            >  1:34316/1‑100 (MQ=255)                                                                     cGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAg                            <  2:2308019/100‑1 (MQ=255)                                                                       aaTGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTg                          >  1:2797395/1‑100 (MQ=255)                                                                        aTGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTgg                         <  2:467784/100‑1 (MQ=255)                                                                        aTGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTgg                         <  1:622874/100‑1 (MQ=255)                                                                           cTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTa                      <  1:2922263/100‑1 (MQ=255)                                                                           cTGACGTCAAACGGGCTGGTGAGAAAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTa                      <  2:1811028/100‑1 (MQ=255)                                                                            tGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTaa                     <  1:2695889/100‑1 (MQ=255)                                                                             gACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAAt                    >  2:1584736/1‑100 (MQ=255)                                                                              aCGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATc                   <  2:2608222/100‑1 (MQ=255)                                                                               cGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATcc                  <  2:2483503/100‑1 (MQ=255)                                                                               cGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATcc                  <  2:1237510/100‑1 (MQ=255)                                                                               cGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATcc                  <  1:1934788/100‑1 (MQ=255)                                                                                gTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCt                 >  1:2917880/1‑100 (MQ=255)                                                                                 tCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTa                <  1:1210595/100‑1 (MQ=255)                                                                                  cAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTaa               >  1:2514774/1‑100 (MQ=255)                                                                                   aaaCGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGGTTCTAAATACTCAGTGGTTAATCCTa                <  2:2500164/98‑1 (MQ=255)                                                                                    aaCGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAAt             >  1:1330901/1‑100 (MQ=255)                                                                                      cGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCg           >  1:361787/1‑100 (MQ=255)                                                                                       gggCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCACTGGTTAATCCTAAATCg           >  2:2785539/1‑99 (MQ=255)                                                                                        ggCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag         <  1:949185/100‑1 (MQ=255)                                                                                        ggCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag         <  2:2549127/100‑1 (MQ=255)                                                                                         gCTGGTGCGACAGAATAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgaga        <  1:2983905/100‑1 (MQ=255)                                                                                         gCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAAt             <  2:1705313/95‑1 (MQ=255)                                                                                         gCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgaga        <  2:1679420/100‑1 (MQ=255)                                                                                         gCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgaga        <  1:1146386/100‑1 (MQ=255)                                                                                          cTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGaa       >  2:1461198/1‑100 (MQ=255)                                                                                            ggTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATa     >  1:2708268/1‑100 (MQ=255)                                                                                              tGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGa   <  2:3011187/100‑1 (MQ=255)                                                                                               gCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGAt  >  2:2221878/1‑100 (MQ=255)                                                                                                  |                                                                                                 GGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGAT  >  NC_000913/2257677‑2257870

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