Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,449,050	T→A	15.1%	intergenic (+32/+26)	feaB → / ← tynA	phenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,449,050	0	T	A	15.1%	253.6 / 34.1	99	intergenic (+32/+26)	feaB/tynA	phenylacetaldehyde dehydrogenase/tyramine oxidase, copper‑requiring
Reads supporting (aligned to +/- strand):  ref base T (51/33);  new base A (8/7);  total (59/40)
Fisher's exact test for biased strand distribution p-value = 7.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTA  >  NC_000913/1448952‑1449145                                                                                                   |                                                                                                ccGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTg                                                                                                >  1:691534/1‑100 (MQ=255) ccGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTg                                                                                                >  2:1017303/1‑100 (MQ=255)  cGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTgg                                                                                               >  1:742567/1‑100 (MQ=255)   gTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACTGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTg                                                                                                <  1:1097675/98‑1 (MQ=255)           ggCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttt                                                                                      <  1:3213413/100‑1 (MQ=255)           ggCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttt                                                                                      <  2:2637920/100‑1 (MQ=255)            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttt                                                                                     >  1:1522766/1‑100 (MQ=255)            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttt                                                                                     >  2:410951/1‑100 (MQ=255)            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttt                                                                                     <  2:717201/100‑1 (MQ=255)            gCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttt                                                                                     <  2:557802/100‑1 (MQ=255)            gCCCCGACTGGCGGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAGTCTGGTTCGCTCATAAGGAAAAACCGGCACCTGGTGCCGttttt                                                                                     <  1:377528/100‑1 (MQ=255)             ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttttt                                                                                    <  1:3067190/100‑1 (MQ=255)             ccccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGtttttt                                                                                    <  2:1357933/100‑1 (MQ=255)              cccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTATGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttttt                                                                                   >  2:1071851/1‑100 (MQ=255)              cccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttttt                                                                                   >  1:1319992/1‑100 (MQ=255)               ccGTCTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGCTCGCTTATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                                                                                  <  1:2271596/100‑1 (MQ=255)               ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttttttt                                                                                   >  2:1551619/1‑99 (MQ=255)               ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                                                                                  <  1:900624/100‑1 (MQ=255)               ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                                                                                  <  2:2161941/100‑1 (MQ=255)               ccGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttttac                                                                                  <  2:482677/100‑3 (MQ=255)                cGACTGGCTGGACGGTTGGTGTGAATCTAAGTCTGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCAGGTGCCGttttttact                                                                                 <  1:1117621/100‑4 (MQ=38)                cGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGt                                                                                 >  2:3013547/1‑100 (MQ=255)                   ctggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCt                                                                               <  1:2290894/99‑1 (MQ=255)                   ctggctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                                                                              >  1:2444161/1‑100 (MQ=255)                       ctggACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAAc                                                                          >  2:127620/1‑100 (MQ=255)                           aCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATc                                                                      >  1:2793575/1‑100 (MQ=255)                            cGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCa                                                                     <  2:1752024/100‑1 (MQ=255)                             gattGGTGTGAAAATAAGTCGGTGTGTGTACAGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCAc                                                                    >  2:2288098/3‑100 (MQ=255)                              gTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACt                                                                   <  2:19975/100‑1 (MQ=255)                                tGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTa                                                                 <  2:1383217/100‑1 (MQ=255)                                  gtgtGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc                                                               <  1:3013521/100‑1 (MQ=255)                                    gtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctt                                                             >  2:1915977/1‑100 (MQ=255)                                    gtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctt                                                             >  1:1138533/1‑100 (MQ=255)                                    gtgAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctt                                                             >  2:875613/1‑100 (MQ=255)                                         aCTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTttcttc                                                        <  2:3282820/100‑1 (MQ=255)                                          cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCa                                                       >  2:1854594/1‑100 (MQ=255)                                          cTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCa                                                       >  1:2236515/1‑100 (MQ=255)                                            aaGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAgc                                                     >  2:1188475/1‑100 (MQ=255)                                            aaGGCGGTGTGTGTACGGTATTAATCTGGTTTGCTTATAAGTAAAAAACGGCACCTGGTGCCGTTTCTTTGTCTGAAACAATCACTTATCTTTCTTCAgc                                                     <  1:2239996/100‑1 (MQ=255)                                                      gtgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCg                                           >  2:2352337/1‑100 (MQ=255)                                                      gtgtACGGTATTAATCTGGTTCGCGCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCg                                           <  1:2994854/100‑1 (MQ=255)                                                       tgtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGt                                          >  2:1698184/1‑100 (MQ=255)                                                        gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                         <  2:3293397/100‑1 (MQ=255)                                                        gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                         >  2:1602814/1‑100 (MQ=255)                                                        gtACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGtt                                         >  2:920282/1‑100 (MQ=255)                                                          aCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTgg                                       >  1:2510686/1‑100 (MQ=255)                                                            ggTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCg                                     >  2:2720992/1‑100 (MQ=255)                                                            ggTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCg                                     <  1:1055453/100‑1 (MQ=255)                                                            ggTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCg                                     <  2:3103730/100‑1 (MQ=255)                                                                 tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                <  1:887165/100‑1 (MQ=255)                                                                 tAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCg                                >  2:2745899/1‑100 (MQ=255)                                                                  aaTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGt                               >  1:2650931/1‑100 (MQ=255)                                                                   aTCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTGTGTCTGAAACAATCGCTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTc                              <  1:88360/100‑1 (MQ=255)                                                                     cTGGTTCGCTCATAATTAAAAAACGGCATCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCCGCGGCCCCAGCGGTGGCGTTCCGGCaa                            >  1:882915/1‑100 (MQ=255)                                                                     cTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCaa                            >  2:303626/1‑100 (MQ=255)                                                                       ggTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAaag                          <  1:3061809/100‑1 (MQ=255)                                                                        gTTTGCTCAAAAGTAAAAAACGGCACCTGGTACCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCCCCCCTAGCGTTGGCGTTTCGTCAaag                          <  1:2154206/99‑1 (MQ=255)                                                                         ttCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCTCCCCTAGCGTTGGCGTTTCGTCAaaga                         >  1:2194602/1‑99 (MQ=255)                                                                           cGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGt                      >  2:685146/1‑100 (MQ=255)                                                                           cGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCATTGGCGTTTCGTCAaagacg                       >  1:1386106/1‑97 (MQ=255)                                                                            gCTCCTAAGTAAAAAACGGCACCTGGTGCCGGTGTTTTCTCTGAAACAATCACTAATCTTTCGTCAGCGCCCCTACCGTTGGCGTCTCGTCAAAGCAGt                      >  1:566382/1‑99 (MQ=255)                                                                             cTCATAAGTAAAAAACGGCACATGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTc                    <  2:85496/100‑1 (MQ=255)                                                                              tCATAAGTAAAAAACTGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTcc                   <  2:3339574/100‑1 (MQ=255)                                                                              tCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTcc                   <  2:1450351/100‑1 (MQ=255)                                                                               cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCa                  >  1:1200942/1‑100 (MQ=255)                                                                               cATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCa                  <  1:171448/100‑1 (MQ=255)                                                                                aTAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCAt                 <  2:943593/100‑1 (MQ=255)                                                                                 agacaaaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGTAGTTCCATg                >  1:3172885/6‑100 (MQ=38)                                                                                 agacaaaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATg                >  2:2624881/6‑100 (MQ=255)                                                                                  gacaaaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGCAACAATCACTTATCTTTCTTCAGCGCTCCTAGCGTTGGCGTTTCGTCAAAGACGTTCCATgg               >  1:1805163/5‑100 (MQ=38)                                                                                  aaGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATgg               >  1:903497/1‑100 (MQ=255)                                                                                  aaGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATgg               >  1:2679617/1‑100 (MQ=255)                                                                                    gTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt             >  1:513982/1‑100 (MQ=255)                                                                                    gTAAAAAACGGCACCTGGTGCCCTTTTTTTGGCTGAACCAATCACTTGTCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt             <  1:2288607/100‑1 (MQ=255)                                                                                    caaaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGtt             >  1:648196/3‑100 (MQ=255)                                                                                     tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGttt            >  1:2602668/1‑100 (MQ=255)                                                                                     tAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGttt            <  2:2929570/100‑1 (MQ=255)                                                                                     aaaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGttt            >  2:37466/2‑100 (MQ=255)                                                                                      aaaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTc           >  1:1575333/1‑100 (MQ=255)                                                                                      aaaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTc           >  2:87057/1‑100 (MQ=255)                                                                                      aaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTc           >  2:2035682/1‑100 (MQ=255)                                                                                      aaaaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTc           >  1:1853762/1‑100 (MQ=255)                                                                                       aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTca          <  2:691537/100‑1 (MQ=255)                                                                                       aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTca          >  2:2714555/1‑100 (MQ=255)                                                                                       aaaaaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTca          >  2:1673496/1‑100 (MQ=255)                                                                                       aaaaaCGACACCTGGTGCCGTATTTTTGTCTGAAACAATCACCTATCTTCCTTCAGCACCCCTAGTATTTGCGTTTCGTCAACGAAATACCAGGGTCTca          >  1:1377292/1‑100 (MQ=255)                                                                                       aaaaaCGACACCTAGTGCCGTATTTTTGTCTGAAACAAGCACTTATCATTCTTCAGCGCACCTAGAGTTGGCGTATCGTCAAAGACGTCCCATgatctc           >  1:2096572/1‑94 (MQ=255)                                                                                        aaaaCGGCACCAGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcag         >  2:1898713/1‑100 (MQ=255)                                                                                          aaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagca       >  1:2544409/1‑100 (MQ=255)                                                                                          aaCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTcagca       <  1:1071848/100‑1 (MQ=255)                                                                                           aCGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCGGCTCCACTAGCGTTGCCGTTTCGTCAAAGAAGTTCCATGGTTTcagcag      >  1:3360916/1‑100 (MQ=255)                                                                                            cGGCACCTGGTGCCGTTTTTTTTTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAATTTCCATGGTTTCAGCaga     >  1:2734936/1‑100 (MQ=255)                                                                                            cGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCaga     >  2:279923/1‑100 (MQ=255)                                                                                             ggCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCagag    <  2:2089005/100‑1 (MQ=255)                                                                                              gCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGGGt   >  1:2544516/1‑100 (MQ=255)                                                                                               cACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCagag    <  1:3038145/98‑1 (MQ=255)                                                                                               cACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTa  >  1:777769/1‑100 (MQ=255)                                                                                               cACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTa  <  2:184415/100‑1 (MQ=255)                                                                                                 ccTGGTGCCGTTTTTTTTTTTTTTTTTATCACTCATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTGGTCAAAGGAGTTCCATTGTTTcagc        >  1:2027621/1‑92 (MQ=255)                                                                                                   |                                                                                                CCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTCTTCAGCGCCCCTAGCGTTGGCGTTTCGTCAAAGAAGTTCCATGGTTTCAGCAGAGTA  >  NC_000913/1448952‑1449145

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