Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,202,012	A→T	16.7%	Q41L (CAG→CTG)	cca →	fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,202,012	0	A	T	16.7%	221.9 / 34.2	90	Q41L (CAG→CTG)	cca	fused tRNA nucleotidyl transferase/2'3'‑cyclic phosphodiesterase/2'nucleotidase and phosphatase
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base T (13/2);  total (54/36)
Fisher's exact test for biased strand distribution p-value = 2.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

GTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTT  >  NC_000913/3201913‑3202111                                                                                                    |                                                                                                    gtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                     <  2:263307/100‑1 (MQ=255) gtgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagca                                                                                                     <  1:244098/100‑1 (MQ=255)  tgCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACcagcag                                                                                                    >  1:2152249/1‑100 (MQ=255)   gCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAAACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAgg                                                                                                   >  1:1807775/1‑100 (MQ=255)     tGTTCGGGATGCATTGTTAGAGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGCGATGCTCGACGCGGGCTACCAGAAGGTa                                                                                                 >  2:981572/1‑100 (MQ=255)      gTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAg                                                                                                >  2:334834/1‑100 (MQ=255)      gTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAg                                                                                                >  1:2488289/1‑100 (MQ=255)        tCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGc                                                                                              <  2:568825/100‑1 (MQ=255)         cGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGcc                                                                                             <  2:338773/100‑1 (MQ=255)           ggATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCcgc                                                                                           <  2:984835/100‑1 (MQ=255)            gATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCcgcg                                                                                          >  1:1461403/1‑100 (MQ=255)              tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAt                                                                                        <  1:2491928/100‑1 (MQ=255)              tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAt                                                                                        <  1:296686/100‑1 (MQ=255)              tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGCCTACCTGCTGGTAGCCcgcgtc                                                                                        <  2:2557758/100‑3 (MQ=38)              tGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGGGGTCGGCAGTACGCCACAGGAGATGTTCGACGCGGGCTACCAGCAGGTAGGCCGCGAt                                                                                        <  2:538141/100‑1 (MQ=255)                ttccgatctggCTACCGGTCAATGACAGAGACTGGGTCGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAttt                                                                                      <  2:1727315/91‑1 (MQ=255)                 aTTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGAtttt                                                                                     >  1:2402652/1‑100 (MQ=255)                   tGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGTCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTcc                                                                                   <  2:1804369/100‑1 (MQ=255)                   tGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTcc                                                                                   >  1:1804383/1‑100 (MQ=255)                   tGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTcc                                                                                   <  2:1767734/100‑1 (MQ=255)                   tGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTcc                                                                                   <  2:618441/100‑1 (MQ=255)                    gTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCt                                                                                  >  2:2384975/1‑100 (MQ=255)                      tAGGGCTACCGGTCAAAGACAGAGATTGGGTGTTGGTCGGCAGTACGCCACACGAGATGCACGACGCGGCCTACCAGCAGCTCGCCCGCGAt                                                                                        >  2:1976215/1‑92 (MQ=255)                         ggCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGtt                                                                             >  1:1888540/1‑100 (MQ=255)                          gCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGttt                                                                            >  2:1065744/1‑100 (MQ=255)                             aCCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                         <  1:835308/100‑1 (MQ=255)                             aCCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTg                                                                         <  1:59037/100‑1 (MQ=255)                                ggTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCAt                                                                      >  1:1071788/1‑100 (MQ=255)                                 gTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATc                                                                     <  2:67138/100‑1 (MQ=255)                                  tCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATcc                                                                    <  2:262404/100‑1 (MQ=255)                                   cAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCg                                                                   >  2:2110798/1‑100 (MQ=255)                                       gACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCaaa                                                               <  1:1065752/100‑1 (MQ=255)                                       gACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCaaa                                                               >  2:2501331/1‑100 (MQ=255)                                       gACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCaaa                                                               >  2:1097426/1‑100 (MQ=255)                                         cAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACg                                                             >  1:2312845/1‑100 (MQ=255)                                           gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                           >  2:361131/1‑100 (MQ=255)                                           gagaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCa                                                           <  1:918031/100‑1 (MQ=255)                                             gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                         >  1:42796/1‑100 (MQ=255)                                             gaTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATg                                                         <  2:1468461/100‑1 (MQ=255)                                              aTTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGa                                                        <  1:334838/100‑1 (MQ=255)                                                tGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAag                                                      >  2:582680/1‑100 (MQ=255)                                                 gggTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAaga                                                     >  1:2117578/1‑100 (MQ=255)                                                   gtggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGt                                                   >  2:1378647/1‑100 (MQ=255)                                                    tggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTa                                                  <  2:776917/100‑1 (MQ=255)                                                    tggtggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTa                                                  <  1:767398/100‑1 (MQ=255)                                                       tggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgc                                               <  1:2183148/100‑1 (MQ=255)                                                        ggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgcg                                              >  2:200492/1‑100 (MQ=255)                                                        ggtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgc                                               >  2:1475785/1‑99 (MQ=255)                                                         gtCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATgcgc                                             <  2:2038532/100‑1 (MQ=255)                                                          tCGGCAGTACGCCACATGAGATGCTCGACGCGCGCTACCAGCAGGTAGGCCGCGATTTTCCTCTCTTTCTTCATCCTCAAACGCATGAAGAGTATGCGCt                                            <  2:2389746/100‑1 (MQ=255)                                                           cGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTACAGTGTTTCTGCATCCGCAAACGAATGAAGAGTAttagct                                            >  2:744157/1‑94 (MQ=255)                                                            ggCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTgg                                          >  1:306833/1‑100 (MQ=255)                                                            ggCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTgg                                          <  2:111763/100‑1 (MQ=255)                                                                gTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACg                                      >  2:1461721/1‑100 (MQ=255)                                                                   cGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAc                                   >  1:682664/1‑100 (MQ=255)                                                                    gCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTAcc                                  >  1:657859/1‑100 (MQ=255)                                                                     ccACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCg                                 >  1:1554959/1‑100 (MQ=255)                                                                     ccACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCg                                 >  1:680250/1‑100 (MQ=255)                                                                          ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                            >  1:1349686/1‑100 (MQ=255)                                                                          ggAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACgg                            >  2:1047495/1‑100 (MQ=255)                                                                           gagaTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGa                           <  1:35459/100‑1 (MQ=255)                                                                             caggaaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACTGATGAAGAGTATGCGCTGGCACGTACCGAACGGaaa                         >  2:1829955/10‑100 (MQ=21)                                                                              aTGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAt                        <  1:200494/100‑1 (MQ=255)                                                                               tGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtc                       >  2:444108/1‑100 (MQ=255)                                                                               tGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtc                       >  2:1967776/1‑100 (MQ=255)                                                                               ggaaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtc                       >  2:2098419/8‑100 (MQ=25)                                                                                gcgggACGCGGGCTACCAGCAGGTAGGCCGCGATGTTCCTGGGTGTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAgcc                      <  2:1701667/96‑1 (MQ=255)                                                                                 cTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                     >  1:1451101/1‑100 (MQ=255)                                                                                 aaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                     >  1:1428768/6‑100 (MQ=37)                                                                                 aaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                     >  2:549786/6‑100 (MQ=37)                                                                                 aaaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccg                     >  1:826403/6‑100 (MQ=37)                                                                                  tCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTTTTTTTTTTTTTCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                    >  2:2484657/1‑100 (MQ=255)                                                                                  aaatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccgg                    >  2:710129/5‑100 (MQ=37)                                                                                   cGACGTGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                   <  2:2038222/100‑1 (MQ=255)                                                                                   aatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                   >  1:2002893/4‑100 (MQ=38)                                                                                   aatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                   >  2:523151/4‑100 (MQ=38)                                                                                   aatcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggt                   <  1:549793/97‑1 (MQ=38)                                                                                    atcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTATCGAAAGGAAAtccgg                    >  2:858532/3‑98 (MQ=37)                                                                                    atcgcgGCCTACCTGCTGGTAGCCCGCGACTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggtt                  >  2:1030840/3‑100 (MQ=38)                                                                                     tcgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttc                 >  1:62632/2‑100 (MQ=38)                                                                                      cgcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCCTGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttcc                >  1:995253/1‑100 (MQ=255)                                                                                      cgcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttcc                >  1:2183517/1‑100 (MQ=38)                                                                                       gcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg               >  1:1359134/1‑100 (MQ=255)                                                                                       gcgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg               <  2:51748/100‑1 (MQ=255)                                                                                       gcgGCCTACCTGCTGGTAGCCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccg               >  1:1204819/1‑100 (MQ=255)                                                                                        cgGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg              >  2:1845567/1‑100 (MQ=255)                                                                                         gggCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccggt             <  1:1311117/100‑1 (MQ=255)                                                                                          ggCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATccggttccggtt            <  1:1633852/100‑1 (MQ=255)                                                                                           gCTCCCAGCAGTTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAAtccggttccgg              <  1:1616036/97‑1 (MQ=255)                                                                                           gCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTa           >  2:128747/1‑100 (MQ=255)                                                                                            cTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTac          >  1:2161510/1‑100 (MQ=255)                                                                                              aCCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTacac        <  1:1147396/100‑1 (MQ=255)                                                                                               ccAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACAcc       <  1:66307/100‑1 (MQ=255)                                                                                                cagcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCg      >  2:2399981/1‑100 (MQ=255)                                                                                                 agcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCgg     >  1:2320071/1‑100 (MQ=255)                                                                                                 agcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCgg     <  2:363701/100‑1 (MQ=255)                                                                                                  gcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGt    <  1:1378656/100‑1 (MQ=255)                                                                                                  gcagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGATCCGGTTACACCGGt    >  1:2525813/1‑100 (MQ=255)                                                                                                   cagGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGtt   <  2:642294/100‑1 (MQ=255)                                                                                                    agGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGttt  >  1:1079067/1‑100 (MQ=255)                                                                                                    |                                                                                                    GTGCTGTTCGGGATGCATTGTTAGGGCTACCGGTCAAAGACAGAGATTGGGTGGTGGTCGGCAGTACGCCACAGGAGATGCTCGACGCGGGCTACCAGCAGGTAGGCCGCGATTTTCCTGTGTTTCTGCATCCGCAAACGCATGAAGAGTATGCGCTGGCACGTACCGAACGGAAATCCGGTTCCGGTTACACCGGTTT  >  NC_000913/3201913‑3202111

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