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.A97.9% 178.9 / ‑3.2 47intergenic (‑237/+77)PP_0224/sctCMonooxygenase, DszC family/sulfur compound ABC transporter ‑ ATP‑binding subunit
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base A (33/13);  total (33/14)
Fisher's exact test for biased strand distribution p-value = 2.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-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.

CGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCGGGCA  >  AE015451/278540‑278966
                                                                                                                                                                                                           |                                                                                                                                                                                                                                 
cGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCt                                                                                                                                                                                                           <  2:473846/227‑1 (MQ=255)
                                    cTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGc                                                                                                                                                                                               >  2:528095/1‑203 (MQ=255)
                                          tCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCTGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCTCGAAAGGGCCGGAACAg                                                                                                                                                        >  1:353842/1‑236 (MQ=255)
                                            acaAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGt                                                                                                                                                      >  1:358655/1‑236 (MQ=255)
                                              aaaaTAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                                                                                                                 <  2:579746/159‑1 (MQ=255)
                                              aaaaTAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                                                                                                                 >  1:579747/1‑159 (MQ=255)
                                                 aTAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTAc                                                                                                                                                                            >  1:3523/1‑209 (MQ=255)
                                                       aTGTGGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAGAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                                                                                                                 <  2:749607/150‑1 (MQ=255)
                                                        tGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCt                                                                                                                                                                                                           <  2:95880/171‑1 (MQ=255)
                                                         gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                                                                                                                 >  2:559750/1‑148 (MQ=255)
                                                         gTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaa                                                                                                                                                                                                                                 <  1:559751/148‑1 (MQ=255)
                                                                         gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                                                                                                               <  2:269214/134‑3 (MQ=255)
                                                                         gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTaaac                                                                                                                                                                                                                               >  1:269215/1‑132 (MQ=255)
                                                                         gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCATAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                                                                                                                   <  2:547431/162‑1 (MQ=255)
                                                                         gcgACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCATAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCt                                                                                                                                                                                                   >  1:547432/1‑162 (MQ=255)
                                                                                        tggtccggtggtccggTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATATGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGt                                                                                                                                                      >  2:302507/1‑192 (MQ=255)
                                                                                                  gtccggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGctc                                                                                                >  1:720466/1‑236 (MQ=255)
                                                                                                     cggtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                                               <  2:590216/170‑1 (MQ=255)
                                                                                                       gtATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGt                                                                                                      >  1:497455/1‑225 (MQ=255)
                                                                                                                cccTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCgg                                                                                                                                              >  2:458947/1‑176 (MQ=255)
                                                                                                                    aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCTGGGCTCTTCGCGTGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGTAGCGGGTTTACCCGCGACAGGGCCGGaa                                                                                                                                                           >  1:184615/1‑159 (MQ=255)
                                                                                                                    aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTcc                                                                                                                                  >  1:500743/1‑184 (MQ=255)
                                                                                                                    aaTGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGCACTGCCGCGTACGCTCTTCCTTCGGCGCTGCaaa                                                                              >  1:740428/1‑236 (MQ=255)
                                                                                                                       gATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGCTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAAGCCGCAGTCCCGAg                                                                                                                              >  2:43691/1‑185 (MQ=255)
                                                                                                                        aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacaca                                                                                                                                                                                                                  <  1:439477/100‑1 (MQ=255)
                                                                                                                        aTTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTacaca                                                                                                                                                                                                                  >  2:439476/1‑100 (MQ=255)
                                                                                                                                     ttATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCt                                                                >  2:609071/1‑233 (MQ=255)
                                                                                                                                                 gttttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGTCTCTTCGCGGGTAACCCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGc                                                        <  1:452449/228‑1 (MQ=255)
                                                                                                                                                    ttGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcct                                                                                           <  1:538856/191‑1 (MQ=255)
                                                                                                                                                                       gggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGGTTACTCGGCGGCCGCAGTCCCGAGGAACTTgcgcca                                                                                                                  >  2:147844/1‑147 (MQ=255)
                                                                                                                                                                        ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                                                                                                                     >  2:74361/1‑65 (MQ=255)
                                                                                                                                                                        ggATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGAt                                                                                                                                                                                                     <  1:74361/65‑1 (MQ=255)
                                                                                                                                                                         gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTcttcct                                                                                           >  1:230554/1‑170 (MQ=255)
                                                                                                                                                                         gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCt                                                                >  2:350291/1‑197 (MQ=255)
                                                                                                                                                                         gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCtt                                                          >  1:685172/1‑203 (MQ=255)
                                                                                                                                                                         gATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGaaa                             >  2:449786/1‑232 (MQ=255)
                                                                                                                                                                            ttGTTACAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCAccc                      <  1:528096/236‑1 (MQ=255)
                                                                                                                                                                                  cAGAGGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCg                >  1:748945/1‑236 (MQ=255)
                                                                                                                                                                                   agagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCAc                                                   >  1:632383/1‑200 (MQ=255)
                                                                                                                                                                                   agagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCAc                                                   <  2:632382/200‑1 (MQ=255)
                                                                                                                                                                                    gagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCg                                                                                        >  2:30483/1‑162 (MQ=255)
                                                                                                                                                                                    gagGTGCCGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCg                                                                                        <  1:30483/162‑1 (MQ=255)
                                                                                                                                                                                          ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCCAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGAGCCGGAACAGGCTTACTCGTCGGcc                                                                                                                                         >  2:26720/1‑107 (MQ=37)
                                                                                                                                                                                          ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGc                                                                  >  2:625374/1‑178 (MQ=255)
                                                                                                                                                                                          ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTcc        >  1:165699/1‑236 (MQ=255)
                                                                                                                                                                                          ccGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCGCTTCCTTCGGCGCTGCAAACAGGGCCGTgg                                                                   >  1:214978/1‑177 (MQ=255)
                                                                                                                                                                                           cGGGCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTccc       <  1:690223/236‑1 (MQ=255)
                                                                                                                                                                                            gggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTg                                                                                  >  2:88648/1‑160 (MQ=255)
                                                                                                                                                                                             ggCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTc                                                      >  2:690222/1‑187 (MQ=255)
                                                                                                                                                                                              gCTCTTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCggg    <  1:88648/236‑1 (MQ=255)
                                                                                                                                                                                                 ctTCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCggg    <  2:748944/233‑1 (MQ=255)
                                                                                                                                                                                                  ttCGCGGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGAGAAAGTGCCGGAACAGGTTTACTCg                                                                                                                                               >  2:692859/1‑93 (MQ=38)
                                                                                                                                                                                                    cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGccct                                        >  2:149780/1‑192 (MQ=255)
                                                                                                                                                                                                    cgcgGGTAAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCGGGCa  >  2:276072/1‑232 (MQ=255)
                                                                                                                                                                                                           |                                                                                                                                                                                                                                 
CGTTACGCACTGGGGTGATTGTGATCCTGACCATTCCTGACCTCACAAAATAGTTATGTCGCTACAGCTGATCGCGACAGACGTGTGCTGGTCCGGTGGTCCGGTATATATACCCTAATGATTTATAAAATTTTTATGAACTAACTTTTTGGAATATTTATAGAAGGGGGATTTGTTACAGAGGTGCCGGGCTCTTCGCGGGTAACCCGCTCCTACACACATGCCTCAAGATCTGTGCACTGTAGGAGCGGGTTTACCCGCGAAAGGGCCGGAACAGGTTTACTCGGAGGCCGCAGTCCCGAGGAACTTGCGCAAGAACTGCCGCGTACGCTCTTCCTTCGGCGCTGCAAACAGGGCCTTGGCCTCGCCTTGCTCCACGATCACGCCTTTGTCGAAGAAAATCACCCGGTTCGCCACGTCCCGGGCA  >  AE015451/278540‑278966

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