Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 4,586,030 2 bp→TC 100% intergenic (+113/+101) 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,0300CT98.0% 165.9 / ‑4.1 49intergenic (+113/+102)PP_4061/PP_4063conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (23/25);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
*AE0154514,586,0310TC98.0% 166.2 / ‑3.3 49intergenic (+114/+101)PP_4061/PP_4063conserved protein of unknown function/putative Long‑chain‑fatty‑acid‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (23/25);  minor base G (1/0);  total (24/25)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GACCCGAGCTTGCCGCCACCCGTACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTT  >  AE015451/4585801‑4586251
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                               
gacCCGAGCTTGCCGCCACCCGTACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgg                                                                                                                                                                                                                            >  1:774397/1‑235 (MQ=255)
           gccgccACCCGTACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGc                                                                                                                                                                                                                   >  1:676041/1‑234 (MQ=255)
                       aCCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGAcc                                                                                                                                                                                                    <  1:1109245/236‑2 (MQ=255)
                                                ccGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGc                                                                                                                                                                           <  1:59786/236‑1 (MQ=255)
                                                              cGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgggtaaa                                                                                                                                                                                                                       <  2:68285/178‑6 (MQ=255)
                                                                  ccgccgACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATgcg                                                                                                                                                                                  <  2:897433/211‑1 (MQ=255)
                                                                        aCAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCt                                                                                                                                                                       <  2:818289/216‑1 (MQ=255)
                                                                            tGGGCGCCCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTAc                                                                                                                                                                                                             <  1:350181/175‑1 (MQ=255)
                                                                                  ccccGCCGGGCCAGGAGCCACCGGTGAAGGCTTACGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgggta                                                                                                                                                                                                                         <  2:176517/155‑4 (MQ=255)
                                                                                  aCCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAAcc                                                                                                                                                                                                                      <  1:298361/160‑1 (MQ=255)
                                                                                  aCCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                                                                                    <  2:861455/162‑1 (MQ=255)
                                                                                  aCCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                                                                                    >  1:861456/1‑162 (MQ=255)
                                                                                    ccgccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAg                                                                                                                                                                                                          <  2:216886/170‑1 (MQ=255)
                                                                                     cgccgAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCAGGGACCGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAggcc                                                                                                                                         <  2:1050340/233‑1 (MQ=255)
                                                                                      gccgGGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                                                                                    <  1:472650/158‑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)
                                                                                                                      ggTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCgggtaa                                                                                                                                                                                                                        >  2:445907/1‑117 (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)
                                                                                                                                                                                gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                                                                                    >  2:331460/1‑68 (MQ=25)
                                                                                                                                                                                gATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCg                                                                                                                                                                                                                    <  1:331461/68‑1 (MQ=25)
                                                                                                                                                                                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)
                                                                                                                                                                                                       gAAGATCTTGCATTGATCGGGCCGGCCTCTTCGCGGGTAAACCCGCt                                                                                                                                                                                                                  >  2:5395/1‑47 (MQ=14)
                                                                                                                                                                                                          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)
                                                                                                                                                                                                                                      ||                                                                                                                                                                                                                               
GACCCGAGCTTGCCGCCACCCGTACCGGAAGAAGACCCGGACGAGCTGCCGATCAAGCCGACCGTGCCGCCGACAGTGGGCGACCCGCCGAGCCAGGAGCCACCGGTGAAGGCTTAAGGGTAGACCACGGCACCGGCGTGAACAGCTGATCTGCGCTGCTCCATTGAAAGGTTGAGGATGACCGGGTGCTTGGCAAGGTGAAGATCTTGCATTGATCGGGCCGGCCTCTCTGCGGTAAACCCGCTCCTACCAGGGACGTGCCGGCCTTGAATGCGGCGCTGCCCCTGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCGCGCCAGCACATCAACCCGCAGATGCAGCCGCAATCTCCGCCACACTGATCTCACGCATGCGGAACTTCTGCACCTTGCCGGTCACGGTCATCGGGTACTCGTCGACAAAACGTATATACCGCGGCACCTT  >  AE015451/4585801‑4586251

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: