Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,586,056:1 +C 100% intergenic (+139/+76) PP_4061 → / ← PP_4063 conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154514,586,0561.C91.3% 148.5 / 8.0 46intergenic (+139/+76)PP_4061/PP_4063conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base . (4/0);  new base C (19/23);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGAT  >  AE015451/4585824‑4586263
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                
aCCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAcc                                                                                                                                                                                                                <  1:1109245/236‑2 (MQ=255)
                         ccGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGc                                                                                                                                                                                       <  1:59786/236‑1 (MQ=255)
                                           ccgccgACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                                                                                                                                                                                              <  2:897433/211‑1 (MQ=255)
                                                 aCAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                                                                                                   <  2:818289/216‑1 (MQ=255)
                                                              cgccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                                                                                                     <  2:1050340/233‑1 (MQ=255)
                                                               gccgGGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCgg                                                                                                                                                                        <  2:1124787/213‑1 (MQ=255)
                                                                     ccAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggccgg                                                                                                                                                   <  1:873815/228‑1 (MQ=255)
                                                                                            aaGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGACGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAgg                                                                                                                                                                             <  1:946782/179‑1 (MQ=255)
                                                                                              gggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGttt                                                                                                                                                                    <  2:988747/186‑1 (MQ=255)
                                                                                                gTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTACACCCGCTCCTACCAGGGACCGTGCCGGCCtt                                                                                                                                                                                                     >  2:425035/1‑151 (MQ=255)
                                                                                                       aCGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATGTTGCACTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgc                                                                                                                                                                                               >  2:1070393/1‑150 (MQ=38)
                                                                                                          gCACCGGCGGGAACAGCTGAGCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAcc                                                                                                                                                                                                                <  2:269196/130‑2 (MQ=38)
                                                                                                            aCCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCTTGCCGGCCTTGCATgcggcg                                                                                                                                                                                           >  2:306325/1‑149 (MQ=38)
                                                                                                                                  cgcTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGc                                                                                                                 <  2:948472/201‑1 (MQ=255)
                                                                                                                                     tgctCCATTGAAAGGTTGAGGATGACCGGGGGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGac                                                                                                                                                           >  2:322638/1‑155 (MQ=38)
                                                                                                                                          cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGAGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTTGGAGCGGGTTTAcc                                                                                                                                                                 >  1:172976/1‑145 (MQ=255)
                                                                                                                                          cATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACGCGCTCCTACCAGCGACTGTGCCGGcc                                                                                                                                                                                                       >  2:899678/1‑107 (MQ=37)
                                                                                                                                           aTTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCaa                                                                                                               <  1:528895/194‑1 (MQ=255)
                                                                                                                                               aaaGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACGCTCTCGCAGCCGCGACTGATCCGGCCTTGAATGcgccgcaccc                                                                                                                                                                                      >  1:830496/1‑115 (MQ=255)
                                                                                                                                                cgggTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGTTCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACAt                                                                                                                                     <  1:681833/165‑1 (MQ=38)
                                                                                                                                                    ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACTCTCTCCTACCAGGGACTGTACCTGCCTtgcgtgcgc                                                                                                                                                                                             >  1:753087/1‑106 (MQ=255)
                                                                                                                                                    ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGaa                                                                                >  1:328342/1‑216 (MQ=255)
                                                                                                                                                    ttGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGCAGCGGGTTTCCCCGCGCAGAGGCCGGCCGCGCCac                                                                                                                                          >  2:59240/1‑157 (MQ=38)
                                                                                                                                                       aGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGAtc                                                                                              <  2:396618/199‑1 (MQ=255)
                                                                                                                                                       aGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGAtc                                                                                              >  1:396619/1‑199 (MQ=39)
                                                                                                                                                         gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcggtcacggt                                                        >  1:595662/1‑235 (MQ=255)
                                                                                                                                                            gACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAACAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                                                                                                   >  1:65218/1‑109 (MQ=37)
                                                                                                                                                                gggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTg                                                                                                                                                                                                    >  2:175958/1‑88 (MQ=21)
                                                                                                                                                                gggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGcccg                                                                                                                                          >  2:373355/1‑144 (MQ=37)
                                                                                                                                                                 ggTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCGGGACCGTGCCGGCCTTGAATGCGCCGCTGCCCCTGGTAGGAGc                                                                                                                                                                          >  1:47071/1‑113 (MQ=37)
                                                                                                                                                                 ggTGCTTGGCAAAGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGaa                                                                                >  2:632470/1‑203 (MQ=39)
                                                                                                                                                                                   gATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAATCCGCTCCTAGCAGGGACTGTGCATGCCTTGGATgc                                                                                                                                                                                               >  1:634920/1‑74 (MQ=9)
                                                                                                                                                                                       ttGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTAc                                              <  1:960940/215‑1 (MQ=255)
                                                                                                                                                                                                   ggccggccTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcgg                                                               <  1:1070394/186‑1 (MQ=255)
                                                                                                                                                                                                     ccggccTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCtt              <  2:595661/233‑1 (MQ=255)
                                                                                                                                                                                                             cttcgcgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGCATATACcg                       >  1:297461/8‑217 (MQ=255)
                                                                                                                                                                                                              ttcgcgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg                                                                                     <  1:897187/148‑1 (MQ=255)
                                                                                                                                                                                                              ttcgcgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATg                                                                                     >  2:897186/7‑154 (MQ=255)
                                                                                                                                                                                                                 gcgggTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATAcc                        <  1:851196/209‑1 (MQ=255)
                                                                                                                                                                                                                       aaaCCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCcggtc                                                             <  1:632471/169‑1 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGcc                                                                                                       >  1:942302/1‑116 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGcc                                                                                                       <  2:942301/116‑1 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACtcgtc                                         >  2:448834/1‑178 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACtcgtc                                         <  1:448835/178‑1 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGtatata                          <  1:479241/193‑1 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTTTGCGAt  >  2:1079809/1‑217 (MQ=255)
                                                                                                                                                                                                                                  tACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGAt  <  1:1079810/217‑1 (MQ=255)
                                                                                                                                                                                                                                     cAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAAc                                 <  2:635040/183‑1 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                
ACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTTGAAGTGTGCGAT  >  AE015451/4585824‑4586263

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: