Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 6,013,885 T→C intergenic (+90/+20) xpt → / ← PP_5266 Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154516,013,8850TC97.0% 107.2 / ‑6.4 34intergenic (+90/+20)xpt/PP_5266Xanthine phosphoribosyltransferase/putative acetyl‑CoA hydrolase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (14/18);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CGATGACTTCCTGGCCAACGGCAAAGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCCGCGCATGCTGTGCCCTGATCGCCT  >  AE015451/6013591‑6014179
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      
ggatcaCTTCCCGGCCAAAGGCAAAGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTATGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAAt                                                                                                                                                                                                                                                                                                  <  1:728334/296‑1 (MQ=255)
            ggCCAACGGCAAAGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGCGCAGCACAGGTCTCGACCGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGCATGGCCGCAACCCt                                                                                                                                                                                                                                                                                      >  1:692013/1‑301 (MQ=255)
                              cccaGCGCTGATCTCGGTCATCAAGCAGCCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGGTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCa                                                                                                                                                                                                                                                                    <  2:655808/298‑1 (MQ=255)
                                   cgcTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGCATGGCCGCAACCCTTTTACCTGGTTGCCGCCAGcccc                                                                                                                                                                                                                                                               >  2:25837/1‑301 (MQ=255)
                                     cTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGc                                                                                                                                                                                                                                                             >  2:545575/1‑301 (MQ=255)
                                       tatctcgcacatcaAGCAGGCGGGTGCCCACGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCca                                                                                                                                                                                                                                                           <  2:197337/292‑1 (MQ=255)
                                                 atcaAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCTCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCCACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAgcc                                                                                                                                                                                                                                                 >  1:703095/1‑300 (MQ=255)
                                                                   aCCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCttt                                                                                                                                                                                                                               <  1:485527/301‑1 (MQ=255)
                                                                     gggggCGGGCCTGGGTATCGTGATCGAGAAGGCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTaa                                                                                                                                                                                                                             <  1:397455/298‑1 (MQ=255)
                                                                                   gtaccgtgctcggcaaGTCGTTCCAGCGCGTCCCTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGCGGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTCCGCCGGCACGCCCGCTCCCACAGGTGCAGCCCAGGTCTCGAAAGCTGTGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGc                                                                                                                                                                                                               <  1:14617/287‑1 (MQ=255)
                                                                                         tGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGCCACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTCCCTGTGGGAGCGGGCGTGCCCGCGAACGGCCGCAACCCTTTTATCTGGTTGCCGCCTGCCCCGCCAATAACAAGCGCTGGTAACGTTCTTCTCTAAGCCCCTGCGGCTGCCCCAAc                                                                                                                                                                                                         >  1:745670/1‑301 (MQ=255)
                                                                                              gagaAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGGCGCTGGAAGGTGGGGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCAt                                                                                                                                                                                                    <  2:877361/301‑1 (MQ=255)
                                                                                                gaAGTCGTTCCAGGAGGGCCTGTCCGAACTGGACAGCCCGGGTTACCGCGTTGAATCGCTGGCACGCGTGAAGTCCCTGGAAGGTGGTGTGGTCAGCGTCATCGGATAACCGCGACAGCTTCGTGGCGGGCACGCCCTCGCCCACAGGGGCATCACAGGTCTCGAACTCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATcc                                                                                                                                                                                                  <  2:44048/301‑1 (MQ=255)
                                                                                                         ccAGGGCCGTCGTGCCCAACTGGACAGCCAGGGTTACCGCGTTGAATCGGAGGCCCGGGTGAAGTCGCTGGAAGGTGCTGTGGTCAGCTTCATCGAATCAACGCGACACCTGGTTCGCGTGCACGCCCGCTCCCACAAGTGCAGCACAGGTCGTGATCGCCGTGCCGGACCTGTGGGAGCGGGCGTTCCTGCGATTGGCCGCAACCCTGTTATCTGGTTGCCGCCAGCCCCGCCAACAGCAAGCGCTGGTAAAGTTCTTCTTTAAACCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGt                                                                                                                                                                                         <  1:85682/301‑1 (MQ=255)
                                                                                                                                    gCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGt                                                                                                                                                               <  1:545761/301‑1 (MQ=255)
                                                                                                                                      cAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGTTGc                                                                                                                                                             >  1:660896/1‑301 (MQ=255)
                                                                                                                                         ggTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTCCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGTCCCCGGTGCATc                                                                                                                                                          >  1:460673/1‑301 (MQ=255)
                                                                                                                                          gTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTCCATcc                                                                                                                                                         >  1:663057/1‑301 (MQ=255)
                                                                                                                                                             ggCCCGAGTGAAGTCGCTGGAAGGTGGGGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCt                                                                                                                                      <  1:114642/301‑1 (MQ=255)
                                                                                                                                                                        ggTCGCTGGAAGTTGGTGTGGTCAGCTTCAGCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCg                                                                                                                           <  2:50053/300‑1 (MQ=255)
                                                                                                                                                                        aGTCGCTCGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTCCCCGTCCGATCGGGCGTGCCCGCGTATTGCCGCAACCCTTTCATCTGGTTCCCGCGAGCGCCGCCAACAACAAGCGCTGGTAACTTTCTTCTTTACGCCTCTACGCCTGCACCAACTGCATCCGCCCCAGGTGGCTTTCGTACCACTCCGGCCCCGGTTCATCCAACGCTACCTCGCCCACCTGCGGCACCTcc                                                                                                                           >  2:259392/1‑300 (MQ=255)
                                                                                                                                                                           cGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGACTGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGATGCATCCAACGCTGCCTTGCCCAGCTCCAGCCCCTCGccc                                                                                                                        >  1:526665/1‑299 (MQ=255)
                                                                                                                                                                             cTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAg                                                                                                                      <  1:25843/301‑1 (MQ=255)
                                                                                                                                                                                 agggtggtGGGGTCAGCTTCATCGAATGAACGCGACCCCTGCTGCGCGGGCACGCCCGCTCCCACAGGTGCAGGACAGGGCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTg                                                                                                                  <  2:662827/299‑1 (MQ=255)
                                                                                                                                                                                 aagggggggtggtCAGCTTCATCGAATGAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTg                                                                                                                  <  2:789397/293‑1 (MQ=255)
                                                                                                                                                                                             tCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGTAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCATCTCGCTCAGCTTGAACTTGCTCTTc                                                                                                      <  1:597044/301‑1 (MQ=255)
                                                                                                                                                                                               cgCTTCATCGAATAAACGGGCGCCCTGCTTCGCGTGCACGCCCGGTTGCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGCGCGGGCGTGCACGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCaa                                                                                                    <  2:156589/300‑1 (MQ=255)
                                                                                                                                                                                                              cccggACACCTGCTTCGCGTGCCCGCCCGCTCCCACAGGTGCAGCACAGGTCTCCAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCg                                                                                     <  2:460519/297‑1 (MQ=255)
                                                                                                                                                                                                                  gACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGATGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAGGTTCTTCTTTAAGCCCCTGCGGCTGCCCCACCTGCATCCGCCCCAGGTGCCTTTCATCCCCCTCCGGCCCCGGTGCATCCACCGCTCCCTTGCCCAGCCCCAGCCCCTCGCTCAGCTTGAACTTGCCCTTTAACCAGCTCACCGCCCCCAAc                                                                                 >  2:426428/1‑301 (MQ=255)
                                                                                                                                                                                                                                gggCAAGCCAGCTCGCACCGGGGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCAGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGctcctc                                                                   <  1:124429/301‑1 (MQ=255)
                                                                                                                                                                                                                                     cGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCg                                                              <  1:64678/301‑1 (MQ=255)
                                                                                                                                                                                                                                            tCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGCTCCCGCTCCTCGGCCGTaaattc                                                       >  1:375520/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                               cGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCAGGGAACAGCCACGCac                    >  2:230840/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                 cGTGCCCGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGCTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCCGCGCCCGCAACAGATCCCGCTCCTCGCCCGTAAAGTCCCTGCCCACCGGCTACTCCGGGAACACCCgcgcctcttgtcccctgctccccc  >  2:237595/1‑281 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                      
CGATGACTTCCTGGCCAACGGCAAAGCCTCCCAAGCGCTGATCTCGATCATCAAGCAGGCAGGTGCCACCGTGGCGGGCCTGGGTATCGTGATCGAGAAGTCGTTCCAGGGCGGCCGTGCCGAACTGGACAGCCAGGGTTACCGCGTTGAATCGCTGGCCCGGGTGAAGTCGCTGGAAGGTGGTGTGGTCAGCTTCATCGAATAAACGCGACACCTGCTTCGCGGGCACGCCCGCTCCCACAGGTGCAGCACAGGTCTCGAACGCTGTGCCGTACCTGTGGGAGCGGGCGTGCCTGCGAATGGCCGCAACCCTTTTATCTGGTTGCCGCCAGCCCCGCCAACAACAAGCGCTGGTAAAGTTCTTCTTTAAGCCCCTGCGGCTGCCCCAACTGCATCCGCCCCAGGTGGCTTTCATACCCCTCCGGCCCCGGTGCATCCAACGCTGCCTTGCCCAGCTCCAGCACCTCGCTCAGCTTGAACTTGCTCTTCAACCAGTTCAGCGCCCGCAACAGATCCCGCTCCTCGGCCGTAAAGTCCGTGCCCATCGGATACTCCGGGAACAGCCGCGCATGCTGTGCCCTGATCGCCT  >  AE015451/6013591‑6014179

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