Predicted mutation
evidence seq id position mutation freq annotation gene description
RA AE015451 278,742:1 +A 100% intergenic (‑237/+77) PP_0224 ← / ← sctC Monooxygenase, DszC family/sulfur compound ABC transporter ‑ ATP‑binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE015451278,7421.A95.4% 159.1 / 0.4 43intergenic (‑237/+77)PP_0224/sctCMonooxygenase, DszC family/sulfur compound ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (36/5);  total (37/6)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-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.

CTCAATGCCCTGGCGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCGG  >  AE015451/278527‑278963
                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
cTCAATGCCCTGGCGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                                                                                                             <  2:313884/218‑1 (MQ=255)
   aaTGCCCTGGCGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                                                                                                           <  2:786983/217‑3 (MQ=255)
       cccTGGCGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCaa                                                                                                                                                                                                     <  1:788870/236‑1 (MQ=255)
                                                    atcTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                                                                                                           <  2:127147/166‑3 (MQ=255)
                                                            aaaTAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCCAGATCTGTGCACTGTAGGAGCg                                                                                                                                                                                >  1:25712/1‑204 (MQ=255)
                                                                   tATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCTAAAGGGCCGGAACAgg                                                                                                                                                    >  1:380220/1‑225 (MQ=255)
                                                                     tGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACCGGTTTACTCGGAGGc                                                                                                                                       >  1:528223/1‑236 (MQ=255)
                                                                          cTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGa                                                                                                                                                         >  2:767430/1‑213 (MQ=255)
                                                                              aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCAcg                                                                                                                                                                                         >  1:314369/1‑176 (MQ=255)
                                                                              aGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTccc                                                                                                                              >  1:725613/1‑236 (MQ=255)
                                                                                                       gtccggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCGGCTCCTACACACATGCCTCAAGATCTGTGCACTGTa                                                                                                                                                                                      >  2:368852/1‑155 (MQ=255)
                                                                                                          cggtggtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCcgcg                                                                                                    >  2:424196/1‑234 (MQ=255)
                                                                                                                    gtgtatataCCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGc                                                                                                                                       >  2:782436/4‑189 (MQ=255)
                                                                                                                    gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTa                                                                                                  >  2:557148/1‑226 (MQ=255)
                                                                                                                                  aTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcctt                                                                                       >  2:613564/1‑223 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGGAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGa                                                                                                                         >  2:142264/1‑187 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGa                                                                                                                         >  2:280903/1‑187 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGa                                                                                                                         <  1:280905/187‑1 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCgg                                                                                    >  2:777567/1‑224 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAg                                                                        >  1:742232/1‑236 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAg                                                                        >  1:428480/1‑236 (MQ=255)
                                                                                                                                    gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTACCcg                                                                                                                                                                      >  1:46305/1‑141 (MQ=255)
                                                                                                                                           aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTg                                                                 >  1:702964/1‑236 (MQ=255)
                                                                                                                                           aaaaTTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGAGGGTAAAACCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAgga                                                                                                                                                   >  2:267419/1‑153 (MQ=255)
                                                                                                                                                                     aaaatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGcc                                                              >  2:433488/4‑213 (MQ=255)
                                                                                                                                                                      atatTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACg                                         <  2:702962/233‑1 (MQ=255)
                                                                                                                                                                         tttATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTAAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTg                                   >  1:117446/1‑236 (MQ=255)
                                                                                                                                                                                     ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCACACAGGGCCt                                                                   >  1:44672/1‑192 (MQ=255)
                                                                                                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCg                                                                                                                                            >  1:279219/1‑118 (MQ=255)
                                                                                                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGCGGa                                                                                                                         >  2:173706/1‑137 (MQ=255)
                                                                                                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGAt                                             >  1:705747/1‑213 (MQ=255)
                                                                                                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGc                                        >  1:521606/1‑218 (MQ=255)
                                                                                                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCGCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCtt                                                       >  1:659286/1‑203 (MQ=255)
                                                                                                                                                                                      gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAGACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCa                      >  1:481078/1‑236 (MQ=255)
                                                                                                                                                                                                   ggTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTTgac                              >  1:442025/1‑214 (MQ=255)
                                                                                                                                                                                                       ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAgg                                                                                                                          >  1:7507/1‑119 (MQ=255)
                                                                                                                                                                                                       ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCa          >  2:577127/1‑231 (MQ=255)
                                                                                                                                                                                                       ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGt       >  2:672514/1‑234 (MQ=255)
                                                                                                                                                                                                       ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTcc     >  1:560303/1‑236 (MQ=255)
                                                                                                                                                                                                       ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGAAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCaa                                                                            >  1:678487/1‑165 (MQ=255)
                                                                                                                                                                                                       ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGAAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCaa                                                                            <  2:678485/165‑1 (MQ=255)
                                                                                                                                                                                                         gggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTc                                                   >  2:600980/1‑188 (MQ=255)
                                                                                                                                                                                                         gggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCg   >  1:776821/1‑236 (MQ=255)
                                                                                                                                                                                                          ggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCgg  <  1:424198/236‑1 (MQ=255)
                                                                                                                                                                                                          gccctttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGGTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGGACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCCCGCCTTTGt                                  >  1:178213/6‑204 (MQ=255)
                                                                                                                                                                                                             tctTCGCGGGTAGACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCgg  <  2:481076/233‑1 (MQ=255)
                                                                                                                                                                                                                        |                                                                                                                                                                                                                              
CTCAATGCCCTGGCGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCGG  >  AE015451/278527‑278963

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