Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,463,652	C→G	6.5%	G353A (GGC→GCC)	treC ←	trehalose‑6‑P hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,463,652	0	C	G	6.5%	327.1 / 11.8	92	G353A (GGC→GCC)	treC	trehalose‑6‑P hydrolase
Reads supporting (aligned to +/- strand):  ref base C (45/41);  new base G (2/4);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.45e-01

AATAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGC  >  NC_000913/4463505‑4463796                                                                                                                                                    |                                                                                                                                                 aaTAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGcatgcc                                                                                                                                                 >  4:745353/1‑149 (MQ=255)      cTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTTAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgca                                                                                                                                            >  3:629064/1‑149 (MQ=255)      cTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATACCATgca                                                                                                                                            <  4:695762/149‑1 (MQ=255)       tcgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                           >  3:141896/1‑149 (MQ=255)       tcgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                           >  3:184521/1‑149 (MQ=255)        cgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcagg                                                                                                                                          <  4:903097/149‑2 (MQ=255)        cgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                           >  2:142307/1‑148 (MQ=255)        cgtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgcag                                                                                                                                           >  1:430556/1‑148 (MQ=255)         gtcgGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcaggg                                                                                                                                         <  4:227908/149‑3 (MQ=255)                tCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGc                                                                                                                                  >  4:308319/1‑149 (MQ=255)                      ccATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgca                                                                                                                                            <  4:616755/133‑1 (MQ=255)                      ccATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATgca                                                                                                                                            >  3:616755/1‑133 (MQ=255)                          cGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttttt                                                                                                                        <  4:184521/149‑1 (MQ=255)                          cGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAttttt                                                                                                                        <  3:593381/149‑1 (MQ=255)                          cGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCAtttt                                                                                                                         <  1:142307/148‑1 (MQ=255)                             tGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGc                                                                                                                      >  1:789868/1‑148 (MQ=255)                                   gCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAgg                                                                                                                >  1:330524/1‑148 (MQ=255)                                       ggCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCAcg                                                                                                            <  2:412133/148‑1 (MQ=255)                                          aaaCATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGt                                                                                                        <  4:418627/149‑1 (MQ=255)                                          aaaCATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGt                                                                                                        <  3:343769/149‑1 (MQ=255)                                                atTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCa                                                                                                   >  1:24417/1‑148 (MQ=255)                                                atTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAc                                                                                                  >  4:448222/1‑149 (MQ=255)                                                  tGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGGGGTATTCACCt                                                                                                <  3:514801/149‑1 (MQ=255)                                                   gAGGCTCTCTACGTCGCGATAGTTAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCtt                                                                                               <  4:629064/149‑1 (MQ=255)                                                   gAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGc                                                                                                                      <  4:207556/126‑1 (MQ=255)                                                   gAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGc                                                                                                                      >  3:207556/1‑126 (MQ=255)                                                    aGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTtc                                                                                              >  3:872262/1‑149 (MQ=255)                                                    aGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTtc                                                                                              >  4:269859/1‑149 (MQ=255)                                                           ctACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCAcc                                                                                        >  2:380352/1‑148 (MQ=255)                                                               gTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAg                                                                                    <  2:430556/148‑1 (MQ=255)                                                               gTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGc                                                                                   >  3:156256/1‑149 (MQ=255)                                                                tCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGc                                                                                   >  1:467677/1‑148 (MQ=255)                                                                 cgcgATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAcc                                                                                                 >  1:173682/1‑133 (MQ=255)                                                                 cgcgATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAcc                                                                                                 <  2:173682/133‑1 (MQ=255)                                                                 cgcgATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCga                                                                                 <  4:466562/149‑1 (MQ=255)                                                                     aTAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGaa                                                                              <  2:89648/148‑1 (MQ=255)                                                                         tCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACaa                                                                          <  1:840764/148‑1 (MQ=255)                                                                         tCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACaa                                                                          <  2:24417/148‑1 (MQ=255)                                                                          cAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcagg                                                                                                                                          <  2:348761/83‑2 (MQ=255)                                                                          cAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGGCATgcagg                                                                                                                                          >  1:348761/1‑82 (MQ=255)                                                                            gTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATgc                                                                       >  4:484375/1‑148 (MQ=255)                                                                            gTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATgc                                                                       <  3:484375/148‑1 (MQ=255)                                                                                    cgTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGa                                                               >  2:100273/1‑148 (MQ=255)                                                                                     gTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATc                                                             <  4:152985/149‑1 (MQ=255)                                                                                       gAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCg                                                            <  2:506991/148‑1 (MQ=255)                                                                                        aaaTGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGt                                                           >  1:325553/1‑148 (MQ=255)                                                                                         aaTGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTgg                                                         <  4:583953/149‑1 (MQ=255)                                                                                          aTGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTgg                                                         <  2:490846/148‑1 (MQ=255)                                                                                          aTGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGt                                                        <  3:745353/149‑1 (MQ=255)                                                                                           tGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGt                                                        >  2:267091/1‑148 (MQ=255)                                                                                               ggTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTaca                                                    <  1:390/148‑1 (MQ=255)                                                                                               ggTTGGTCATGCCAATCTCCTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTacac                                                   <  4:616179/149‑1 (MQ=255)                                                                                                  tGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCa                                                 >  2:111602/1‑148 (MQ=255)                                                                                                    gTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGaa                                              >  4:7185/1‑149 (MQ=255)                                                                                                    gTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGaa                                              >  4:11455/1‑149 (MQ=255)                                                                                                     tCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAAc                                             >  3:817144/1‑149 (MQ=255)                                                                                                       aTGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                           <  4:539823/149‑1 (MQ=255)                                                                                                       aTGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACa                                            <  2:467677/148‑1 (MQ=255)                                                                                                        tGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                           <  1:52366/148‑1 (MQ=255)                                                                                                        tGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAg                                          >  4:395914/1‑149 (MQ=255)                                                                                                         gCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAg                                          >  1:910493/1‑148 (MQ=255)                                                                                                          ccAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAgg                                         >  2:491717/1‑148 (MQ=255)                                                                                                          ccAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAgg                                         >  2:609529/1‑148 (MQ=255)                                                                                                           cAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCa                                       >  4:306603/1‑149 (MQ=255)                                                                                                            aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                           <  3:627200/144‑1 (MQ=255)                                                                                                            aaTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACaa                                           >  4:627200/1‑144 (MQ=255)                                                                                                             aTCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCAt                                      >  2:261903/1‑148 (MQ=255)                                                                                                               ctctTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTc                                    >  2:363478/1‑148 (MQ=255)                                                                                                                     gCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTa                             >  3:630200/1‑149 (MQ=255)                                                                                                                            tAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAcc                                                                                                 <  1:684663/74‑1 (MQ=255)                                                                                                                            tAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCAcc                                                                                                 >  2:684663/1‑74 (MQ=255)                                                                                                                            tAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTtgtg                      <  3:306603/149‑1 (MQ=255)                                                                                                                            tAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTtgt                       <  2:325553/148‑1 (MQ=255)                                                                                                                            tAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTtgt                       >  2:151509/1‑148 (MQ=255)                                                                                                                                 ataCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTc                 <  3:7185/149‑1 (MQ=255)                                                                                                                                  taCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTc                 <  2:321238/148‑1 (MQ=255)                                                                                                                                   aCGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTcc                >  1:570064/1‑148 (MQ=255)                                                                                                                                    cGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCt               >  1:341144/1‑148 (MQ=255)                                                                                                                                     ggCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCtt              <  2:570064/148‑1 (MQ=255)                                                                                                                                      gCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCttgt            <  4:156256/149‑1 (MQ=255)                                                                                                                                       cGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCttgt            <  1:380352/148‑1 (MQ=255)                                                                                                                                         ttCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTGATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGc         >  4:148789/1‑149 (MQ=255)                                                                                                                                           tgcTGCATGGCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGa                                               >  1:330422/3‑109 (MQ=255)                                                                                                                                           tgcTGCATGGCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGa                                               <  2:330422/107‑1 (MQ=255)                                                                                                                                           cccTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATgc                                                                       >  4:606963/1‑85 (MQ=255)                                                                                                                                           cccTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATgc                                                                       <  3:606963/85‑1 (MQ=255)                                                                                                                                             cTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGt     >  3:858844/1‑149 (MQ=255)                                                                                                                                             cTGCATGCCATGCAGCACCACCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCa                    >  3:15728/1‑134 (MQ=255)                                                                                                                                             cTGCATGCCATGCAGCACCACCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCa                    <  4:15728/134‑1 (MQ=255)                                                                                                                                              tGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTg    <  4:872262/149‑1 (MQ=255)                                                                                                                                              tGCATGCCATGCAGCACCATCGCAAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGt     >  1:774896/1‑148 (MQ=255)                                                                                                                                               gCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCAt                                      <  2:80543/114‑1 (MQ=255)                                                                                                                                               gCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCAt                                      >  1:80543/1‑114 (MQ=255)                                                                                                                                                catgccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGc  >  3:587714/1‑149 (MQ=255)                                                                                                                                                 atgccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGc  >  2:851803/1‑148 (MQ=255)                                                                                                                                                 atgccatgcAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGc  >  2:767413/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                 AATAACTCGTCGGCATCACGCCCATCGTTGCGCAGCTCGGCAAACATATTGAGGCTCTCTACGTCGCGATAGTCAGTAATGCGCGTGAAATGCGGGTTGGTCATGCCAATCTCTTCGCCCTGGTAGATATACGGCGTTCCCTGCATGCCATGCAGCACCATCGCCAGCATTTTTGCCGCAGGCACGCGGTATTCACCTTCATCACCAAAGCGAGAAACAATGCGCGGCTGATCGTGGTTACACCAGAACAAGGCATTCCATGCTACGTTGTGCATTCCTTGTTGCCAGTGGC  >  NC_000913/4463505‑4463796

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