Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 499,203 A→G intergenic (+47/‑249) uhpA → / → PP_0411 DNA‑binding response regulator UhpA/Polyamine ABC transporter, ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451499,2030AG97.0% 106.3 / ‑5.0 33intergenic (+47/‑249)uhpA/PP_0411DNA‑binding response regulator UhpA/Polyamine ABC transporter, ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (15/17);  minor base . (0/1);  total (15/18)
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

TGGCCGACCGCATGGTGTTCGAAGTAGGCCTCACCGAATCCCGGCCGTTGCACACACTGCTCTCGGAGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCAGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACAAGAAGGTGCAGGCATGAGCGAGGCGAAGTCGAACGCTGCAACCGGCGAAACGCTGGTCA  >  AE015451/498912‑499497
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                       
gggCCGACCGCATGGTGTTCGAAGTAGGCCTCACCGAATCCCGGCCGTTGCACACACTGCTCTCGGAGCGTGGATTTTCGGTGTTCGAGCGCCTGGCACAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGccc                                                                                                                                                                                                                                                                                                <  2:409003/300‑1 (MQ=255)
 ggCCGACCGCAGTGTGTTCGAAGTAGGCCTAACGGAATCCCGGCCGTTGCACACACTGGTCTCGGAGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTAGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCcg                                                                                                                                                                                                                                                                                               <  2:82052/301‑1 (MQ=255)
              gtgtTCGAAGTAGGCCTCACCGAATCCCGGCCGTTGCACACACTGCTCTCGGAGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGggcctgg                                                                                                                                                                                                                                                                                  <  1:81459/301‑1 (MQ=255)
                      aGTAGGCCTCACCGAATCCCGGCCGTTGCACACACTGCTCTCGGAGCGGGAATTTTCGGTTTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTAATACAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCACGAACTGGCCTGGCctgctgt                                                                                                                                                                                                                                                                          >  1:82825/1‑300 (MQ=255)
                                            ccGTTGCACACACTGCTCTCGGAGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGACTGCTGATCTCTATGGCTGCGCCGGCctc                                                                                                                                                                                                                                                    >  2:475413/1‑301 (MQ=255)
                                                    cacacTGCTCTCGGGGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCggg                                                                                                                                                                                                                                            <  1:75573/301‑1 (MQ=255)
                                                       acTGCTCTCGGAGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCAt                                                                                                                                                                                                                                         >  1:680425/1‑301 (MQ=255)
                                                                              cGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCACGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGAATGCCTGCTCCCACAGGTTTCgggg                                                                                                                                                                                                                  >  1:308839/1‑298 (MQ=255)
                                                                                     ggAGCGCCTGGCCCAGGGCGGCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTc                                                                                                                                                                                                          <  2:680135/300‑1 (MQ=255)
                                                                                            cTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAggg                                                                                                                                                                                                   <  1:56874/301‑1 (MQ=255)
                                                                                             tGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCCGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAgggg                                                                                                                                                                                                  >  2:219039/1‑300 (MQ=255)
                                                                                                           aaCGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAATGGGCCTGGCCTGCTGCTCTCTATGGCTGAGCCGGCCTCTTCGGGGGCATTCCCGCTCCCACCGGTTTCGGGCTTGGCCTCGCACGGAC‑TTCCAATCCCggc                                                                                                                                                                                    >  2:35339/1‑299 (MQ=255)
                                                                                                              gTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATc                                                                                                                                                                                 >  1:450513/1‑301 (MQ=255)
                                                                                                                 aaCGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCtt                                                                                                                                                                              >  1:110541/1‑301 (MQ=255)
                                                                                                                    gACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCCTCGCGTGCATGCCCGCTCCCACAGGTTTCGTTGTTGG‑CTCGCAATGGCATTCTAATCCCGCCATCCTTGtt                                                                                                                                                                           >  2:458452/1‑300 (MQ=255)
                                                                                                                     acgtcgccccgcggtcgacccTGTCGAGCAAGACCATCAGCCCCCACCAGGCACGCCTGATGCGGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGTAGCCCAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGG‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAg                                                                                                                                                                          <  2:479917/281‑1 (MQ=255)
                                                                                                                       cgctgccagagggtGGCGCGGTCGAGCAAGACCAGCAGCACCACCAAGGCAATCCGGATGCAGAAG‑TTGAAAGTGAACTCGCTGGAGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGAGGCCCGAAACTTGCCGAAGCCTATGGGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGATCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAggg                                                                                                                                                                        <  2:556908/288‑1 (MQ=255)
                                                                                                                                 agTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTac                                                                                                                                                              >  1:748298/1‑301 (MQ=255)
                                                                                                                                           gTCGTGCAAGCCCATCACCACCCACAAGGCACGCCTGCTTCAGAAGTTTGAGA‑TGAACTCGCCGGCGGCGCTGGTGATGTACGCCATGGAGCACAAGCTGATCTGAGTGCGACATACGGGTGCCCGACACTTGCCGAAGCCGTTGTGGGG‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCtt                                                                                                                                                    <  2:497801/301‑1 (MQ=255)
                                                                                                                                                     aCCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCCCCCTCCTTGTAGGGCGATCCCTACAACAAATCTTCCATACggtt                                                                                                                                          >  2:310014/1‑299 (MQ=255)
                                                                                                                                                                  acaAGGTACGCATGATGCCGAAG‑TTGAAAGTGAACTCGATGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAACCTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCACGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATtctct                                                                                                                             <  2:308720/301‑1 (MQ=255)
                                                                                                                                                                     gggCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGATGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAg                                                                                                                          <  1:247006/300‑1 (MQ=255)
                                                                                                                                                                           gCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATTGCATTCTCGCAGTCcccc                                                                                                                    >  1:251759/1‑301 (MQ=255)
                                                                                                                                                                           gCCTGATGCAGAAG‑GTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTAcccc                                                                                                                    <  1:475615/301‑1 (MQ=255)
                                                                                                                                                                             cTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCCATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTGTGGTGGTTGG‑CTCGCATGAGCATTCCAATCCCGCCATCCTTTTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTTAGTACCCCCg                                                                                                                  >  1:183372/1‑301 (MQ=255)
                                                                                                                                                                                                                   tGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGACATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCt                                                                             <  1:280807/301‑1 (MQ=255)
                                                                                                                                                                                                                              gCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATc                                                                  >  1:200328/1‑301 (MQ=255)
                                                                                                                                                                                                                              gCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCACCGATTTTCGGGCGCTTGCGTCTCTTCTCATCTTTTGCCTACGCAGTCATc                                                                  >  2:417192/1‑301 (MQ=255)
                                                                                                                                                                                                                                        acaAGCTAATCTGATTGCGAACTACGGGTGCCGCAAACTTGCCGAAGACTATGTGGGA‑GCGGGCATGGCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACaagaag                                                        <  1:387194/301‑1 (MQ=255)
                                                                                                                                                                                                                                          aaGCGGATCTCACTGCGACATACGGGTGCACGAAACTTGCCGAAGCAGATGTGTGA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACAAGAAGGt                                                      <  1:70374/301‑1 (MQ=255)
                                                                                                                                                                                                                                               tatctgttTGCGACATAAGGGTGCCCTAAGATTGGCGAAGCCGTTGTGGGACGC‑GGCATGCCCGGGAACGGGCCTGGCCTGCTGATCGCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACAAGAAGGTGCAgg                                                 <  2:450323/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       ccTATGTGAGA‑GCGCGCATGCCCGCGAACGGGCCTGGCCTGCTGAGATCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACAAGAAGGTGCAGGCATGAGCGAGGCGAAGTCGAACGCTGCAACCGGCGAAACg         <  2:251659/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              gggA‑GCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCAGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACAAGAAGGTGCAGGCATGAGCGAGGCGAAGTCGAACGCTGCAACCGGCGAAACGCTGGTCa  <  1:219124/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                                                                       
TGGCCGACCGCATGGTGTTCGAAGTAGGCCTCACCGAATCCCGGCCGTTGCACACACTGCTCTCGGAGCGGGAATTTTCGGTGTTCGAGCGCCTGGCCCAGGGCGCCAACGTCAACGACATCGCCCAGCAGTTGGCGCTGTCGAGCAAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAG‑TTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGA‑GCAGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGG‑CTCGCAAGGGCATTCCAATCCCGCCATCCTTGTAGGGCGATCCCTACAACAAATCTTCCATCCGGATGATGGCATTCTCTCAGTACCCCCGATTTTCGGGCGCTTGCGTCTCTTCTACTCTTTTGCCTACGCAGTCATCCAACAAGAAGGTGCAGGCATGAGCGAGGCGAAGTCGAACGCTGCAACCGGCGAAACGCTGGTCA  >  AE015451/498912‑499497

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