Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,491,098	0	A	C	20.3%	133.9 / 11.8	59	Y66S (TAC→TCC)	cybB	cytochrome b561
Reads supporting (aligned to +/- strand):  ref base A (31/16);  new base C (0/12);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 2.72e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCA  >  NC_000913/1490961‑1491229                                                                                                                                          |                                                                                                                                    tATCGCAGCGTTTTGCGCAAGGGATTTTCGTGTTTTCTCCCCACGTAGCGATCGGCCCCTCATCAACATGTTTCATTTTTCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATccc                                                                                                                                   <  2:1118337/140‑1 (MQ=255)    cGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCaa                                                                                                                                >  1:881293/1‑140 (MQ=255)       aGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAAccc                                                                                                                             >  2:1144980/1‑140 (MQ=255)         cGTATTGCGCAATGGAGTTTGGGGGTTTCTTCCCTCGTGGCGATCGGCCACTCATCAAAATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATACCCACCCccg                                                                                                                           <  2:420219/140‑1 (MQ=255)         cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATAccc                                                                                                                                  >  1:428000/1‑133 (MQ=255)         cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATAccc                                                                                                                                  <  2:428000/133‑1 (MQ=255)         cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccg                                                                                                                           <  2:216478/140‑1 (MQ=255)         cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccg                                                                                                                           <  2:397183/140‑1 (MQ=255)         cGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCccg                                                                                                                           >  1:278237/1‑140 (MQ=255)         cGTATTGCGAAATGGGTTTTCGTGTTTTCTTCCCACGTAGCGATCGGCCACTCACCAACATTATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTAATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCACCCccg                                                                                                                           <  2:278237/140‑1 (MQ=255)          gTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCccgc                                                                                                                          <  2:810283/140‑1 (MQ=255)              tGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGAt                                                                                                                      >  2:296429/1‑140 (MQ=255)              tGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGAt                                                                                                                      >  2:1142525/1‑140 (MQ=255)               gcgcAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGAtt                                                                                                                     <  1:669971/140‑1 (MQ=255)                        aGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAg                                                                                                            >  1:1148618/1‑140 (MQ=255)                              gTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCaaaa                                                                                                      <  2:1158766/140‑1 (MQ=255)                                  tttCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCgg                                                                                                  >  1:197467/1‑139 (MQ=255)                                   ttcttcCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCgat                                                                                                 <  2:999455/140‑1 (MQ=255)                                    tcttcCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCgatg                                                                                                <  2:1034702/140‑1 (MQ=255)                                            cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACggg                                                                                        >  1:583084/1‑140 (MQ=255)                                            cGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACggg                                                                                        >  1:429491/1‑140 (MQ=255)                                                gCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTg                                                                                    >  2:693651/1‑140 (MQ=255)                                                 cGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTGCGTCTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTgg                                                                                   <  2:881293/140‑1 (MQ=255)                                                   aTCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGgcg                                                                                 >  1:480062/1‑140 (MQ=255)                                                   aTCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGgcg                                                                                 >  1:363721/1‑140 (MQ=255)                                                      ggCCCCTCATCAACATGATTCATGTTTCCTGGGGCATCTCATTCCTCGTGCTTATGGTCCTTCGTCTTTTTTAAAGGTTGAAGTCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAt                                                                              <  2:197467/140‑1 (MQ=255)                                                       gCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCAtt                                                                             >  1:203655/1‑140 (MQ=255)                                                           cTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTggg                                                                         >  1:963098/1‑140 (MQ=255)                                                              atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACa                                                                      >  1:202569/1‑140 (MQ=255)                                                              atcaACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACa                                                                      >  1:5383/1‑140 (MQ=255)                                                                        ttCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGAtt                                                            >  1:549737/1‑140 (MQ=255)                                                                          cATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTa                                                          >  2:116135/1‑140 (MQ=255)                                                                           ttttttCCTGTGGCATCTAAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                         <  2:963098/137‑1 (MQ=255)                                                                           aTGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                         <  2:1148618/140‑1 (MQ=255)                                                                           aTGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAt                                                         <  2:1190420/140‑1 (MQ=255)                                                                             ttttCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAtct                                                       <  2:111811/139‑1 (MQ=255)                                                                             gTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTAtct                                                       <  1:361342/140‑1 (MQ=255)                                                                                          tttcAAGTCTCGGGCTGTGGGTCGTTGGTTTTTTTTTAAGGCTGAAACCCCTACCCCCGCCGATTATACCTAAGCCAAAGCCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTgcg                                          <  2:413591/138‑1 (MQ=255)                                                                                               aTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                      <  1:19909/139‑1 (MQ=255)                                                                                               aTTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGc                                      >  2:19909/1‑139 (MQ=255)                                                                                                   tcGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTg                                 >  1:882513/1‑140 (MQ=255)                                                                                                     gTGCTGATGGTGGTTCGTCTTTTGTTAGGGCTAAAATACCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                               <  1:8728/140‑1 (MQ=255)                                                                                                     gTGCTGATGGTCGTTCGTTTTTTGTTAAGGCTGAAATCCCCACCCCCGCCGATTATACCTAAACCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                               <  1:1155056/140‑1 (MQ=255)                                                                                                     gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                               >  1:827405/1‑140 (MQ=255)                                                                                                     gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                               >  2:873166/1‑140 (MQ=255)                                                                                                     gTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAt                               >  2:890807/1‑140 (MQ=255)                                                                                                        cTGATGGTCGTTCGTTTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttgg                            <  2:5383/140‑1 (MQ=255)                                                                                                        cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttgg                            <  2:202569/140‑1 (MQ=255)                                                                                                        cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTg                                            >  1:589900/1‑124 (MQ=255)                                                                                                        cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTg                                            <  2:589900/124‑1 (MQ=255)                                                                                                        cTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttgg                            >  1:875139/1‑140 (MQ=255)                                                                                                         tGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGAttggt                           <  2:429491/140‑1 (MQ=255)                                                                                                               tcgttcgtTTTTTTTTAGGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtg                     <  2:955880/140‑1 (MQ=255)                                                                                                                cgttcgtCTTTTGTTAAGGTTGAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                    <  1:296429/140‑1 (MQ=255)                                                                                                                cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                    <  1:1142525/140‑1 (MQ=255)                                                                                                                cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                    <  1:116135/140‑1 (MQ=255)                                                                                                                cgttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                    <  2:549737/140‑1 (MQ=255)                                                                                                                  ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                  >  1:877758/1‑140 (MQ=255)                                                                                                                  ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                    >  2:1047703/1‑138 (MQ=255)                                                                                                                  ttcgtCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtga                    <  1:1047703/138‑1 (MQ=255)                                                                                                                  ttcgtCTTTTGTTAAGGCTAAAATCCCCACCCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGtgatg                  <  2:583128/140‑1 (MQ=255)                                                                                                                        ttttGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGtat            >  1:215549/1‑140 (MQ=255)                                                                                                                           tGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAAc         <  1:890807/140‑1 (MQ=255)                                                                                                                               aaGGCTGAAATCCCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCggg     <  1:693651/140‑1 (MQ=255)                                                                                                                                 ggCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGc   >  2:148882/1‑140 (MQ=255)                                                                                                                                  gCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCa  >  2:841151/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                    TATCGCAGCGTATTGCGCAATGGAGTTTCGTGGTTTCTTCCCACGTAGCGATCGGCCACTCATCAACATGATTCATGTTTCCTGTGGCATCTCAATTCTCGTGCTGATGGTCGTTCGTCTTTTGTTAAGGCTGAAATACCCAACCCCGCCGATTATACCTAAGCCAAAACCGATGATGACGGGACTGGCGCATTTGGGACATTTGGTGATTTATCTCCTGTTTATTGCGCTGCCAGTGATTGGTTTGGTGATGATGTATAACCGGGGCA  >  NC_000913/1490961‑1491229

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