Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,604,3500CT34.3% 45.2 / 14.7 35S64F (TCT→TTT) bglJDNA‑binding transcriptional regulator BglJ
Reads supporting (aligned to +/- strand):  ref base C (2/21);  new base T (12/0);  total (14/21)
Fisher's exact test for biased strand distribution p-value = 1.09e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATC  >  NC_000913/4604111‑4604589
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
aGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGaa                                                                                                                                                                                                                                          <  2:310743/251‑1 (MQ=255)
             acttaacttGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCagaga                                                                                                                                                                                                                             <  1:190770/251‑1 (MQ=255)
                      gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATGCCTATACCAGGCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTGTGCGGCGGTCATTTTTTCTTTTTCTCC‑CATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                                                                                                                     <  2:407895/250‑1 (MQ=255)
                      gCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCAAAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTAt                                                                                                                                                                                                                     <  2:129842/250‑1 (MQ=255)
                         ttCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGC‑CCTGAGAAGTGAGCGCAGAGAGGGGTTATCttt                                                                                                                                                                                                                 >  1:333306/1‑250 (MQ=255)
                                                  tGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTTTTC‑CATTAGGAGTGAGGGCAGGGGGGGGTTTTCTAGCCTGAATGAACTTGCGGTTacattt                                                                                                                                                                                        >  2:294324/1‑246 (MQ=255)
                                                                  tAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGC‑CATGAGAAGTGGGCGCAGAGAGGGTTTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑tt                                                                                                                                                                       >  1:146579/1‑251 (MQ=255)
                                                                            aatttCGCACTCATAGAAAAATGCGTCATGAGGAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCAGTTCAGTATGCCATGTCGCGGATCTCGTTTGCGTCGGTCATTGTTTCTTGGTCTGC‑CATGAGAAGTGAGAGCAGAGATGGATTATCTTTCCTTACTTAGCTCGCTATTAACTTTGCGCGTACCCGGCG‑TTTATGTATTGc                                                                                                                                                             <  2:267538/247‑1 (MQ=255)
                                                                                       cATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTTATTTTTT‑TTTTTTTGCTCTTTAGAAGTGAGCGCAGAGGGGGGTTATTTTGTATTAGCGGACTTGCGATTAAAGTTCCCAGTGCCCGGCG‑TTTAGTTATTTCTggtgatgat                                                                                                                                                   >  2:302014/1‑250 (MQ=255)
                                                                                           gAAAAATGCGTCATGAGTAGTATCGGTATTGAGGGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGa                                                                                                                                              <  1:41439/251‑1 (MQ=255)
                                                                                            aaaaaTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTTC‑CATGAGGAGGGAGCGCCGGGAGGGGTTTTCTTTCCTGACTGAACTGGCGGTTAAGTTTTCGCGTACCCGGCG‑TTTTGTTATTGCGGATGATGATAttta                                                                                                                                              >  1:128219/1‑248 (MQ=255)
                                                                                                aTGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTC‑CCTGAGAAGTGGGCGCAGAGAGGGAGTATATTTCCTGACTGAACTGGGGCTTAAGGTTCCTCGGACCCGGCG‑TTTAGTTTTTGCgggtggtggtttttaagacc                                                                                                                                         >  1:278678/1‑233 (MQ=255)
                                                                                                             agtaTCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGTTCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTTCCCGGCG‑TTTAGTTATTGCGGATGATGTTATTGAAGCCCGGCTGATT‑GGTTc                                                                                                                           <  2:241547/251‑1 (MQ=255)
                                                                                                               taTCGGCATTGAGCGTTTATTCAAAAAGTTTGCGGGTAACCCTAATACGCTACATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGCAATTTTTTTTTTTTTTGC‑CCTGCGAAGTGCGCGCAGCAAGGGGTTATATTGCCAGACTGGACTGCCGATTTAGTTTCCGCGCACCCGGCGGTTTA‑TTCTTTCGGGTGTTGCTATTGTAGCGCGGCTGATT‑TGtttat                                                                                                                         >  2:384965/1‑248 (MQ=255)
                                                                                                                        tGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTGTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGTCTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGAGGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCa                                                                                                                 <  1:28794/250‑1 (MQ=255)
                                                                                                                         gagaggttaggcAGAGAGTGTGCGGGTAACTCTTCTCAGCTCCGGACCTATACCGGTCAGTTGTCAGTTCAGGATGCCTGGTTGCGGAGCTGGGTTGGGGCGGTCATTTTTTCTTTTTCTGC‑CTTGAGAATTGAGCGCAGAGAGTGGTTATGTTGCGTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGTG‑TTGAGTTATTGCGTATGATGAAATTGAATGTCGGCTGATT‑CGTTCATTGTCGCCAt                                                                                                                <  2:313764/239‑1 (MQ=255)
                                                                                                                          agagTTTATTCAGAAAGATTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGGCATTTTTTTTTTTTTTGT‑CATGAGAAGGGAGCGCCGAGAGGGATTATTTTTCCTTACTGGACTTGAGGTTAAGTTTTCACGGAGCCGGCG‑TTTAGTTTTTTCGGGTGGAGATATTGGAGGGCGGCTGGTT‑TGTTCATTTTCGCCATcc                                                                                                              >  2:135573/1‑250 (MQ=255)
                                                                                                                               ttATTCAGAAAGTGTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGc                                                                                                         <  1:302155/250‑1 (MQ=255)
                                                                                                                                   tCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTGC‑CCTGAGAAGGGGGCGCAGAGGGGGGTTATCTTTCCTGACTGGACTGGGGATTAAGTTTCCGCGTAACCGGCG‑CTTTGTTTTTTCGGGTGGTGATTTTTAAAGTCGGGTGATT‑GGTTCCTTTTCGCCCTC‑ACCGCTgg                                                                                                      >  2:178078/1‑249 (MQ=255)
                                                                                                                                              gCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAg                                                                                         <  2:374211/251‑1 (MQ=255)
                                                                                                                                                gggTAACCCTTATAAACTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTCTTC‑CATGAGAAGGGAGCGCAGAGAGGGAATATCTTTACTGACTGAACTGGCGATTAAGGTTCCGCGGACCCGGCG‑GTTAGGTATTGGGGGTGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAgg                                                                                        >  1:319659/1‑249 (MQ=255)
                                                                                                                                                          tataAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCa                                                                               <  2:338600/249‑1 (MQ=255)
                                                                                                                                                                    aTACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCGCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGAtt                                                                   <  2:11984/251‑1 (MQ=255)
                                                                                                                                                                       ccTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGAttttt                                                                <  2:305641/251‑1 (MQ=255)
                                                                                                                                                                            aCCAGTTAGGAGTCATTTCAGGCTGCCATGTCGCGGATCTCGTTTGCGGCGGTCCTTTTTTTTTTTTTTGC‑CCTGAGAAGGGAGCGCAGCGGGGGGTTATCTTTCCTTACTGAACTTGCGGTTAAGTTTCCGCGTACCCGGCG‑TTTTGTTATTGCGGGTTGTGATATTGAAGCGCGGCTGATTTGGTT‑ATTTTCGCCATCCACC‑CTGGGCGGGGTATTAAGTAAAGCGGCAACGCTTGGGATTTTtt‑atcc                                                          >  1:114710/1‑250 (MQ=255)
                                                                                                                                                                                gTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTTTTTTTTTTC‑CCTGGGAAGGGAGCGGAGGGAGGGATTTTTTTTCCTGGCTGAAATGGGGCTTTAGTTTTCTCGTACCCCGCG‑TTTTGTTTTTATGGGTGCGGATATTGTAACTCGGCTGATT‑GGTTTTTTGTTGCCATC‑ACACCTGGGCGGGGTTTTAAGTAAAGCGTGAACGCTGGGGATTTTTT‑ATCAGGa                                                       >  2:195170/1‑250 (MQ=255)
                                                                                                                                                                                        tCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGACATTTTTTTTTTTTTTGC‑CATGAGGAGAGTGCGCCGCGAGGGGTTATATTTCCTCACTGAACTTGGGGGTATAATTCCGCGTCCACGGCG‑TTTATTTATTTCGGATACTGATATTGGAGCACAGCTGTTT‑GGGTTATTTTCGGCCTT‑ACCGCCGGGCGGGGGATTAAGGAAAACGTCAACCCTGGGGATTTTTCTAT‑AGGAACTCttttt                                              >  1:303493/1‑251 (MQ=255)
                                                                                                                                                                                          aTTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGt                                            <  2:235118/251‑1 (MQ=255)
                                                                                                                                                                                                     ccATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGgt                                 <  1:423801/251‑1 (MQ=255)
                                                                                                                                                                                                      cATGTCGCGGATCGCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGgtg                                <  2:184522/251‑1 (MQ=255)
                                                                                                                                                                                                       aTGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGgtgt                               <  2:66844/251‑1 (MQ=255)
                                                                                                                                                                                                         ttCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGTCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATTATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑TCCGCTGGACGGTGTATTAAGTAAATCGTCAATGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTa                              <  1:386007/249‑1 (MQ=255)
                                                                                                                                                                                                         gTCGCGGATCGCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCGCGGCTGATT‑GGTTCATTGGCGCCATC‑ACCGCTGGACGGTTTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTAc                             <  1:271517/251‑1 (MQ=255)
                                                                                                                                                                                                            gcgGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTc                          <  2:333306/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTg         >  1:344685/1‑250 (MQ=255)
                                                                                                                                                                                                                                     ttttctttttctGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATc  >  1:95621/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    
AGGGTGCAAGAGAACTTAACTTGCCTTCTGTTTATGGAGATGCCGCAGAATGGAACACAGCCGAATTAAGAAGAGAAATGTCGCACTCATAGAAAAATGCGTCATGAGTAGTATCGGTATTGAGAGTTTATTCAGAAAGTTTGCGGGTAACCCTTATAAGCTCCATACCTATACCAGTCAGGAGTCATTTCAGGATGCCATGTCGCGGATCTCGTTTGCGGCGGTCATTTTTTCTTTTTCTGC‑CATGAGAAGTGAGCGCAGAGAGGGATTATCTTGCCTGACTGAACTGGCGATTAAGTTTCCGCGTACCCGGCG‑TTTAGTTATTGCGGATGATGATATTGAAGCTCGGCTGATT‑GGTTCATTGTCGCCATC‑ACCGCTGGACGGTGTATTAAGTAAAGCGTCAACGCTGGAGATTTTTC‑ATCAGGAACTCTTTTTGTCATTAAATGGTGTACGTCAGGCGACCGACCGACTGAACAATC  >  NC_000913/4604111‑4604589

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