Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,808,173 G→A 30.6% intergenic (‑282/‑273) PP_4234 ← / → dsbD‑II Aldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,808,1730GA30.6% 67.5 / 39.4 49intergenic (‑282/‑273)PP_4234/dsbD‑IIAldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (12/22);  new base A (6/9);  total (18/31)
Fisher's exact test for biased strand distribution p-value = 7.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTG  >  AE015451/4807944‑4808361
                                                                                                                                                                                                                                     |                                                                                                                                                                                            
gCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCG                                                                                                                                                                                        <  1:418151‑M1/236‑1 (MQ=255)
        aaaGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgctc                                                                                                                                                                                <  1:407287‑M1/236‑9 (MQ=255)
                                      aTTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaaca                                                                                                                                                   <  1:34856‑M1/235‑41 (MQ=255)
                                       ttGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttc                                                                                                                                                  <  1:717504‑M1/235‑39 (MQ=255)
                                         gggtgggtAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgg                                                                                                                                                >  1:478475‑M1/1‑195 (MQ=255)
                                                 agagGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctc                                                                                                                                                                                >  1:860415‑M1/1‑182 (MQ=255)
                                                        aaGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaaccc                                                                                                                                <  1:977408‑M1/236‑60 (MQ=255)
                                                         aGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaa                                                                                                                                                                       <  1:831583‑M1/196‑18 (MQ=255)
                                                           tCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggta                                                                                                                                                                        >  2:63167‑M1/1‑177 (MQ=255)
                                                              aGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccc                                                                                                                                         <  1:360923‑M1/221‑51 (MQ=255)
                                                                 tGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaa                                                                                                                            <  1:910879‑M1/231‑64 (MQ=255)
                                                                        gtgtCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaa                                                                                                                                                  >  1:973572‑M1/1‑161 (MQ=255)
                                                                                 tgtgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgac                                                                                                       <  1:626424‑M1/236‑85 (MQ=255)
                                                                                   tgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcg                                                                                                                                                                            <  2:634177‑M1/165‑13 (MQ=255)
                                                                                     tataGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatg                                                                                                   <  1:63169‑M1/236‑86 (MQ=255)
                                                                                     tataGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGGTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactct                                                                                                                    <  2:443213‑M1/219‑72 (MQ=255)
                                                                                              aTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaat                                                                                                                                                         >  1:391576‑M1/1‑139 (MQ=255)
                                                                                              aTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaat                                                                                                                                                         <  2:391573‑M1/173‑35 (MQ=255)
                                                                                              aTGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgt                                                                                                                                                         >  2:967773‑M1/1‑142 (MQ=255)
                                                                                                gggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatc                                                                                         <  1:179496‑M1/235‑96 (MQ=255)
                                                                                                  gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacct                                                                                                                                           >  1:699700‑M1/1‑133 (MQ=255)
                                                                                                  gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacct                                                                                                                                           <  2:699697‑M1/183‑51 (MQ=255)
                                                                                                  gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaa                                                                                                                         >  2:489947‑M1/1‑135 (MQ=255)
                                                                                                  gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGc                                                                                                                                                                                       >  2:498941‑M1/1‑138 (MQ=255)
                                                                                                       aTTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGTGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcggccagtccc                                                                                                                                        <  2:591510‑M1/181‑49 (MQ=255)
                                                                                                                     gatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCCAAAttcaa                                                                                             <  1:62589‑M1/98‑3 (MQ=255)
                                                                                                                     gatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTccacc                                                                     <  1:19803‑M1/122‑3 (MQ=255)
                                                                                                                                         gTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggag                                               >  1:744854‑M1/1‑99 (MQ=255)
                                                                                                                                                  ctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGGAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCCCTCaaaa                                                                                                  <  2:125031‑M1/93‑1 (MQ=255)
                                                                                                                                                     tGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgact                                                         <  2:478472‑M1/214‑128 (MQ=255)
                                                                                                                                                            cgtcgtacaaacaaggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGAccc                                  >  1:125033‑M1/74‑230 (MQ=255)
                                                                                                                                                                 gCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCGCGGGcggcgcttcgcgggtaaacccgc                                                                                                                                                                 <  1:410620‑M1/98‑24 (MQ=255)
                                                                                                                                                                 gCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCGCGGGcggcgcttcgcgggtaaacccgc                                                                                                                                                                 >  2:410617‑M1/1‑75 (MQ=255)
                                                                                                                                                                           gttcAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgactg                                                        <  2:13163‑M1/192‑129 (MQ=255)
                                                                                                                                                                                       ggAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgacccacccccaacacccgcgttaaccttcatttaatcaaccccagccagcatggccccattcccagatctgccaaggcgaagacatgcgcgttctcctgc     <  1:936461‑M1/232‑183 (MQ=255)
                                                                                                                                                                                       gactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTg  >  1:415166‑M1/47‑235 (MQ=255)
                                                                                                                                                                                                   agggactggccgaaacacacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGaa                     >  1:962472‑M1/35‑204 (MQ=255)
                                                                                                                                                                                                    gggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                  <  1:679281‑M1/141‑1 (MQ=255)
                                                                                                                                                                                                    gggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                  >  2:679278‑M1/34‑174 (MQ=255)
                                                                                                                                                                                                    gggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTtact         <  2:415163‑M1/182‑1 (MQ=255)
                                                                                                                                                                                                        gagTTTCAGCGCCTGGGAGATCGAGCGCCGCccgcgcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactct                           <  1:498944‑M1/193‑163 (MQ=255)
                                                                                                                                                                                                          ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAACCACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaa                                                                              <  2:51421‑M1/113‑1 (MQ=255)
                                                                                                                                                                                                          ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCACAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTc     <  2:401835‑M1/186‑1 (MQ=255)
                                                                                                                                                                                                          ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAActc              >  1:377349‑M1/28‑204 (MQ=255)
                                                                                                                                                                                                          ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAActc              <  2:377346‑M1/177‑1 (MQ=255)
                                                                                                                                                                                                             ccaGCGCCAGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgccgcccctcgaaaatccaacaacaaagaccctgcaaacccagaa                                                                                                                            <  2:106679‑M1/90‑64 (MQ=255)
                                                                                                                                                                                                               aaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATcc                                                >  2:570714‑M1/23‑165 (MQ=255)
                                                                                                                                                                                                               aaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATcc                                                <  1:570717‑M1/143‑1 (MQ=255)
                                                                                                                                                                                                                      gggAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgcgggaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaaca                                                       >  2:336935‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                   ccgccCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgacccacccccaacacccgcgttaaccttcatttaatcaaccccagccagcatggccccattcccagatctgccaaggcgaagacatgc                 <  1:39095‑M1/176‑171 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                            
GCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTG  >  AE015451/4807944‑4808361

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