Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,808,173 G→A 20.7% 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,1730GA20.7% 117.3 / 31.4 58intergenic (‑282/‑273)PP_4234/dsbD‑IIAldehyde oxidase and xanthine dehydrogenase/thiol/disulfide interchange protein 2
Reads supporting (aligned to +/- strand):  ref base G (13/33);  new base A (3/9);  total (16/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.64e-01

CCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCG  >  AE015451/4807940‑4808381
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                
ccTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCC                                                                                                                                                                                                                >  1:1279588‑M1/1‑236 (MQ=255)
       tGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcg                                                                                                                                                                                                          <  1:858548‑M1/235‑8 (MQ=255)
           aaaaGCTGACCTTGGCGTCAACTTGACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgct                                                                                                                                                                                                     <  1:79064‑M1/236‑8 (MQ=255)
            aaaGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcg                                                                                                                                                                                                       <  2:1232788‑M1/233‑6 (MQ=255)
                                 tttACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCC                                                                                                                                                                                                               >  2:653547‑M1/1‑204 (MQ=255)
                                              ggtgggtAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgc                                                                                                                                                                                                               >  2:60456‑M1/1‑189 (MQ=255)
                                              ggtgggtAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCGCGGGcggcgcttcgcgggtaaacccgctcctacgtttgtttcggcc                                                                                                                                                                  <  1:669062‑M1/236‑43 (MQ=255)
                                                    tAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcggccagtccc                                                                                                                                                            <  1:782333‑M1/236‑49 (MQ=255)
                                                     agcggcaaAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggcagtccc                                                                                                                                                            <  1:1100835‑M1/232‑49 (MQ=255)
                                                     agagGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctg                                                                                                                                                           >  1:971819‑M1/1‑184 (MQ=255)
                                                        ggCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacct                                                                                                                                                               <  1:935948‑M1/229‑51 (MQ=255)
                                                            aaGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaa                                                                                                                                                                   <  1:1073279‑M1/221‑45 (MQ=255)
                                                                      ggcaaTGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccg                                                                                                                                               <  2:1279585‑M1/229‑62 (MQ=255)
                                                                       aCAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaa                                                                                                                                                        <  1:467693‑M1/221‑56 (MQ=255)
                                                                              gtCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatct                                                                                                                                                                                               >  1:350708‑M1/1‑157 (MQ=255)
                                                                                    gtgtgtATAGGGCGATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacc                                                                                                                                                                                        <  1:210850‑M1/176‑21 (MQ=255)
                                                                                    gtgtgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGGCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctga                                                                                                                            <  1:689280‑M1/236‑84 (MQ=255)
                                                                                    gtgtgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgttt                                                                                                                                                                       <  1:820665‑M1/193‑38 (MQ=255)
                                                                                    gtgtgtATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgc                                                                                                                                                 <  1:1070434‑M1/215‑60 (MQ=255)
                                                                                                   tGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaaca                                                                                                                                              <  1:166636‑M1/203‑66 (MQ=255)
                                                                                                    gggCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTTGGAGGTGTTCGCCGGGGAATTTTCAGGGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctggatctcacaggg                                                                                                                                                                                        >  1:1203400‑M1/1‑135 (MQ=255)
                                                                                                      gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacac                                                                                                                                                                 >  2:1216351‑M1/1‑133 (MQ=255)
                                                                                                      gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgg                                                                                                                                                                    >  2:1251875‑M1/1‑138 (MQ=255)
                                                                                                      gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgg                                                                                                                                                                    <  1:1251878‑M1/178‑41 (MQ=255)
                                                                                                      gCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacg                                                                                                                                                                              >  2:492435‑M1/1‑138 (MQ=255)
                                                                                                         gAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactct                                                                                                                                        <  1:63870‑M1/203‑72 (MQ=255)
                                                                                                          aaTTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaa                                                                                                                                    <  2:1149817‑M1/206‑73 (MQ=255)
                                                                                                               gcgcTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgc                                                                                                                                                                                                      <  2:493478‑M1/135‑7 (MQ=255)
                                                                                                                                       cGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttc                                                                                                                                                                      <  2:1029209‑M1/143‑39 (MQ=255)
                                                                                                                                                   aGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgact                                                                             <  2:398376‑M1/220‑128 (MQ=255)
                                                                                                                                                       ccTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggc                                                                                           <  2:953868‑M1/202‑114 (MQ=255)
                                                                                                                                                          gggAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcct                                                                                                                              <  1:348625‑M1/164‑82 (MQ=255)
                                                                                                                                                           ggAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaagg                                                                                            <  1:703288‑M1/197‑113 (MQ=255)
                                                                                                                                                                     gCCGGGGGATTTGCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgttt                                                                                                                                                                       <  2:362572‑M1/112‑38 (MQ=255)
                                                                                                                                                                             aggcggtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGACGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCt                                                                                                                                                 >  2:250479‑M1/61‑126 (MQ=255)
                                                                                                                                                                             aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaagattccaacgcagcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaa                                                                                                                                             <  2:1221874‑M1/130‑67 (MQ=255)
                                                                                                                                                                             aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgacccacccccaacacccgcgttaaccttcattt                                                                                            <  1:1233938‑M1/179‑116 (MQ=255)
                                                                                                                                                                             aTTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgacccacccccaacacccgcgttaaccttcattt                                                                                            >  2:1233935‑M1/1‑64 (MQ=255)
                                                                                                                                                                                gtcAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccaattattcctgcgggatttgcgcgcgaccgccttggggcatgtctcgatatcgcgtcgtacaaacaaggcggtcgcgcgc                                                                                 <  1:21217‑M1/186‑124 (MQ=255)
                                                                                                                                                                                 agtcgcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                                      <  2:1213167‑M1/141‑1 (MQ=255)
                                                                                                                                                                                 agtcgcgcgcgactgctgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAAt                                                                      >  1:1213170‑M1/57‑197 (MQ=255)
                                                                                                                                                                                     gcgcgcgactggtgcaggagggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACaaa                                                           <  1:176703‑M1/152‑1 (MQ=255)
                                                                                                                                                                                     cgcCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgacc                                 <  1:60456‑M1/230‑177 (MQ=255)
                                                                                                                                                                                        cTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaaga                                                                       >  2:454781‑M1/1‑51 (MQ=255)
                                                                                                                                                                                                    cGCCGGGAGCGTGTCAGCGCCGGGGAGATCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaa                                                                                                                                                <  2:38017‑M1/104‑64 (MQ=255)
                                                                                                                                                                                                        gggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAaca        <  1:1204623‑M1/203‑1 (MQ=255)
                                                                                                                                                                                                        gggactggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACACCTGACCCACCCCCaa          <  1:659578‑M1/201‑1 (MQ=255)
                                                                                                                                                                                                              ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCACAATTcc                                                                                                                  <  2:937869‑M1/97‑1 (MQ=255)
                                                                                                                                                                                                              ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAAcgc                                                                                                             >  2:420109‑M1/28‑129 (MQ=255)
                                                                                                                                                                                                              ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAAcgc                                                                                                             <  1:420112‑M1/102‑1 (MQ=255)
                                                                                                                                                                                                              ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCaa                                                                                                  <  2:409927‑M1/113‑1 (MQ=255)
                                                                                                                                                                                                              ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCt                                                                             <  1:164596‑M1/134‑1 (MQ=255)
                                                                                                                                                                                                              ggccgaaacaaacgtaggagcgggtttACCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCcgcg  <  1:1246008‑M1/209‑1 (MQ=255)
                                                                                                                                                                                                                            gagaTCGAGCGCCGCCCcctcaaaattccaacgccgcctcccctcgaaaagccaacaacaaagaccctgcaaacccagaacaaactcttactactcctgacccacccccaacacccgcgttaaccttcatttaatcaaccccagccagcatggccccattcccagatctgccaaggcgaagacatgcgcgttctcctgctcttcctcacactcctactggc   >  2:881075‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                             agaTCAAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaactccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaat                                                                                    <  2:640819‑M1/139‑126 (MQ=255)
                                                                                                                                                                                                                               aTCGAGCGCCGCCCGCGcggcgctcgatctcacagg                                                                                                                                                                                         >  2:503559‑M1/1‑17 (MQ=255)
                                                                                                                                                                                                                                  gAGCGCCGCgcgggcggcgctcgatctcacaggcgctgaaaatcctccggcaaacacctcaaaattccaacgccgtccactcaaaattccaacgccgccccctcaaaattccaacgccgcctcccctcgaaaatccaacaacaaagaccctgcaaacccagaacaaactct                                               <  1:112056‑M1/171‑163 (MQ=255)
                                                                                                                                                                                                                                  gAGCGCCGCCCGCGcggcgcttcgcgggtaa                                                                                                                                                                                           <  2:269208‑M1/31‑18 (MQ=255)
                                                                                                                                                                                                                                      gccgccCGCGcggcgcttcgcgggtaaacccgctcctacgtttgtttcgggccagttattcctgtgggacgcgcacgcgagccccctggcgcacgcctcgatatcgcgtcgtacaaacaaggcggtcgcgcgcgcctggcacaggagggattggcca                                                         <  1:506724‑M1/157‑148 (MQ=255)
                                                                                                                                                                                                                                      gccgccCGCGcgccgcttcgcgggtaaacccgctcctacgtttgtttcggcc                                                                                                                                                                  <  2:137619‑M1/52‑43 (MQ=255)
                                                                                                                                                                                                                                         |                                                                                                                                                                                                                
CCTGGCTTGAGAAAAGCTGACCTTGGCGTCAACTTTACAGCAATTGGGTGGGTAGAGGCAAAGTCAAGCTGACAATGTGTCGCGGTGTGTATAGGGCAATGGGCGGAATTCGCGCTTTGGCGGGTAGGTGTTTGGCGGGAGGTTTTCAGCGCCTGGGAGGTGTTCGCCGGGGGATTTTCAGCGCCTGGGAGGTGTTCGCCGGGAGAGTTTCAGCGCCTGGGAGATCGAGCGCCGCCCGCGAAGCGCCGCGCGGGCGGCGCTCGATCTCACAGGCGCTGAAAATCCTCCGGCAAACACCTCAAAATTCCAACGCCGTCCACTCAAAATTCCAACGCCGCCCCCTCAAAATTCCAACGCCGCCTCCCCTCGAAAATCCAACAACAAAGACCCTGCAAACCCAGAACAAACTCTTACTACTCCTGACCCACCCCCAACACCCGCG  >  AE015451/4807940‑4808381

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