Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,247,144	T→A	5.2%	D78V (GAT→GTT)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,247,144	0	T	A	5.2%	294.7 / 2.4	88	D78V (GAT→GTT)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base A (3/2);  total (44/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTTTCTTC  >  NC_000913/1246997‑1247286                                                                                                                                                    |                                                                                                                                               gTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGCGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAt                                                                                                                                                <  4:166628/148‑1 (MQ=255)   tAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAt                                                                                                                                                >  1:137090/1‑146 (MQ=255)       aCCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGc                                                                                                                                            >  1:169590/1‑146 (MQ=255)            ccAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                                     <  4:319781/148‑1 (MQ=255)              aaaGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatcat                                                                                                                                   >  4:103126/1‑148 (MQ=255)                aGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGatcat                                                                                                                                   >  4:386803/1‑146 (MQ=255)                 gTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAg                                                                                                                                <  8:6277/148‑1 (MQ=255)                  tCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGc                                                                                                                               <  6:169413/148‑1 (MQ=255)                       cAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCGAGGATCATCAGCGGATc                                                                                                                          >  7:403347/1‑148 (MQ=255)                                 gcgcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggc                                                                                                                <  5:276318/148‑1 (MQ=255)                                  cgcAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggca                                                                                                               <  5:271405/148‑1 (MQ=255)                                      gTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacggc                                                                                                           <  2:151230/148‑1 (MQ=255)                                       tGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTcggcacggca                                                                                                          >  8:274659/1‑148 (MQ=255)                                           tGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                                     >  3:392639/1‑117 (MQ=255)                                           tGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGatc                                                                                                                                     <  4:392639/117‑1 (MQ=255)                                            gCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCgg                                                                                                      >  8:284093/1‑147 (MQ=255)                                               ccTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCaaa                                                                                                  <  7:415087/148‑1 (MQ=255)                                                  ctGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               >  1:422407/1‑148 (MQ=255)                                                  ctGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGATAATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               >  6:49155/1‑148 (MQ=255)                                                     gCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               <  1:113262/145‑1 (MQ=255)                                                      cGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATc                                                                                                        >  2:420935/1‑135 (MQ=255)                                                      cGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATc                                                                                                        <  1:420935/135‑1 (MQ=255)                                                       ggAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCGAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGtttttt                                                                                          <  8:403347/148‑1 (MQ=255)                                                           cATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGCTCGGCACGGCATCGGCAAAGGTTTTTTGGTc                                                                                      >  7:415501/1‑148 (MQ=255)                                                              atTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGCCAAAGGTTATTTGGGCCg                                                                                    >  8:329685/1‑147 (MQ=255)                                                              atTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGATAACCAAATGGCGCAGATCAAATCCGCTCTGGTTATGCTGAATCCGATAATCAGAAAGGATCATCAGCGTATCGCTGATCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                                    >  2:354/1‑147 (MQ=255)                                                               tttCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCgg                                                                                   <  7:274659/147‑1 (MQ=255)                                                                ttCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGa                                                                                  >  4:300318/1‑147 (MQ=255)                                                                    cGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAg                                                                             >  5:949/1‑148 (MQ=255)                                                                     gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATTGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  2:169590/145‑1 (MQ=255)                                                                     gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  2:422407/145‑1 (MQ=255)                                                                     gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAg                                                                             <  6:949/147‑1 (MQ=255)                                                                     gCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  5:49155/145‑1 (MQ=255)                                                                       cttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                                    <  3:48340/138‑1 (MQ=255)                                                                       cttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCg                                                                                    >  4:48340/1‑138 (MQ=255)                                                                        ttcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTg                                                                         >  3:72731/1‑148 (MQ=255)                                                                         tcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGttttt                                                                                           <  4:290355/129‑1 (MQ=255)                                                                         tcttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGttttt                                                                                           >  3:290355/1‑129 (MQ=255)                                                                         acttACGGCAGGGTGAAATAGACGTAAACTAAATGGCCCAGATCAAATCCGCTCTCGTTTTGCTGCATCCGATAAACAGCAAGCATCATCAGCGGATCGCTCTACGGCACGGAATCGGCAAAGGTTTTTTGCTCCCGAAAAAGATTgg                                                                        >  2:362301/2‑148 (MQ=255)                                                                          cttTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTg                                                                         <  2:137090/146‑1 (MQ=255)                                                                             tCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAAACAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGt                                                                     <  3:386803/147‑1 (MQ=255)                                                                                  aGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGca                                                               <  4:72731/148‑1 (MQ=255)                                                                                     gTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTg                                                                 >  6:402831/1‑143 (MQ=255)                                                                                     gTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTg                                                                 <  5:402831/143‑1 (MQ=255)                                                                                             gACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTATTGCACATCATTaaa                                                     >  3:306615/1‑147 (MQ=255)                                                                                                 ttAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCg                                                >  7:440362/1‑148 (MQ=255)                                                                                                  tAACGAAATGGCGCAGATCAAATCCGCACTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGTCATCGGCAAAGGTTTTTAGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCgg                                               >  5:11634/1‑148 (MQ=255)                                                                                                   aaCGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCgg                                               <  4:409733/147‑1 (MQ=255)                                                                                                    aCGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGc                                              <  7:284093/147‑1 (MQ=255)                                                                                                      gAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     <  7:258602/138‑1 (MQ=255)                                                                                                      gAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     >  8:258602/1‑138 (MQ=255)                                                                                                      gAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGccc                                            >  8:92952/1‑147 (MQ=255)                                                                                                       aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATGAGCGGCTCGCTGTTCGGCAAGGCATCGGCAAAGGTTTTTTGGTCCGGAAATAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCta                                          <  1:354/148‑1 (MQ=255)                                                                                                       aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTc                                                                                                                   <  7:376166/75‑1 (MQ=255)                                                                                                       aaaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTc                                                                                                                   >  8:376166/1‑75 (MQ=255)                                                                                                        aaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               >  7:88987/1‑94 (MQ=255)                                                                                                        aaTGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGt                                                                                               <  8:88987/94‑1 (MQ=255)                                                                                                             cgcAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATaaaa                                    >  8:398825/1‑148 (MQ=255)                                                                                                               cAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAata                                  >  6:82882/1‑148 (MQ=255)                                                                                                                 gATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATc                                <  5:82882/148‑1 (MQ=255)                                                                                                                  aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGc                                                                >  5:207569/1‑115 (MQ=255)                                                                                                                  aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGc                                                                <  6:207569/115‑1 (MQ=255)                                                                                                                  aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAatat                                 >  5:261443/1‑146 (MQ=255)                                                                                                                  aTCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAatat                                 >  8:242102/1‑146 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaa                                                      <  2:64513/124‑1 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaa                                                      >  1:64513/1‑124 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  8:58930/1‑148 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  4:408765/1‑148 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  8:13158/1‑148 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGACATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAg                              >  5:163189/1‑148 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGAACGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               >  1:30820/1‑99 (MQ=255)                                                                                                                   tCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGAACGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGaaaa                                                                               <  2:30820/99‑1 (MQ=255)                                                                                                                       aTCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     >  3:276981/1‑121 (MQ=255)                                                                                                                       aTCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTaaa                                                     <  4:276981/121‑1 (MQ=255)                                                                                                                        tCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAggcggc                         <  7:13158/148‑1 (MQ=255)                                                                                                                            ctctGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCCCGGCATCGGCAAAGGTTTTTTGGTCCGGCATAAGGTTGGCGTTTTGCACAGCATTACCCAGCGGCCCTAATAAAATATCAGGCGGCtgtg                     <  1:362301/148‑1 (MQ=255)                                                                                                                            ctctGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCtgtg                     <  3:300318/148‑1 (MQ=255)                                                                                                                               tGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGCTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTgtgg                    <  8:415501/146‑1 (MQ=255)                                                                                                                                ggTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGgtgt                 >  7:161649/1‑148 (MQ=255)                                                                                                                                 gTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGgtgt                 <  8:161649/147‑1 (MQ=255)                                                                                                                                    ttGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTAc              <  7:242102/147‑1 (MQ=255)                                                                                                                                     tgctgcATCCGATAATCAGCAAGGATCTTCAGCGGATCGCTGTTCGGCACGCCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCGCTAATAAAATATCAGGCGGCTGTGGTGTTACCg            >  6:331063/1‑148 (MQ=255)                                                                                                                                      gctgcATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTcaa                                                     <  2:307990/106‑1 (MQ=255)                                                                                                                                      gctgcATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTcaa                                                     >  1:307990/1‑106 (MQ=255)                                                                                                                                          cATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGgtg         <  8:440362/146‑1 (MQ=255)                                                                                                                                          cATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGtt       <  7:126624/148‑1 (MQ=255)                                                                                                                                             ccGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAat                                   <  3:399884/117‑1 (MQ=255)                                                                                                                                             ccGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAat                                   >  4:399884/1‑117 (MQ=255)                                                                                                                                              cGATAATCAGCAAGGATCATCAGCGGATCGGTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTttctt   <  8:97636/148‑1 (MQ=255)                                                                                                                                               gATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTttcttc  >  3:112238/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                               GTTAATACCGGCCAAAGTCCGTCAATATGTTCGCGCAGTGACTGCCCCTCTGGCGGAACATATTTCTCGCCTTCTTTCGGCAGGGTGAAATTGACGTTAACGAAATGGCGCAGATCAAATCCGCTCTGGTTTTGCTGCATCCGATAATCAGCAAGGATCATCAGCGGATCGCTGTTCGGCACGGCATCGGCAAAGGTTTTTTGGTCCGGAAAAAGTTTGGCGTTTTGCACATCATTAAACAGCGGCCCTAATAAAATATCAGGCGGCTGTGGTGTTACCGGTGTTTCTTC  >  NC_000913/1246997‑1247286

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