Predicted mutation
evidence seq id position mutation annotation gene description
RA AE015451 3,137,784 (T)9→8 intergenic (+212/‑78) PP_2753 → / → PP_2754 putative Branched‑chain amino acid ABC transporter, ATP‑binding protein/putative OprD family outer membrane porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AE0154513,137,7760T.98.0% 224.7 / ‑1.5 50intergenic (+204/‑86)PP_2753/PP_2754putative Branched‑chain amino acid ABC transporter, ATP‑binding protein/putative OprD family outer membrane porin
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base . (24/25);  total (24/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CATCCTGGAAAACGGCCGGGTGGTAGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATGGG  >  AE015451/3137476‑3138072
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                        
cATCCTGGAAAACGGCCGGGTGGTAGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGt                                                                                                                                                                                                                                                                                                          >  4:650651/1‑301 (MQ=255)
         aaaCGGCCGGGTGGTAGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttttttttgt                                                                                                                                                                                                                                                                                                 >  3:221957/1‑299 (MQ=255)
          aaCGGCGGGGTGGGAGGGGAGGGCGGTGCTGCGGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttttttttgtt                                                                                                                                                                                                                                                                                                <  2:183086/301‑4 (MQ=255)
           aCGGCCGGGTGGTAGGTGAGGGAAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGttttttttgttg                                                                                                                                                                                                                                                                                               <  4:303943/301‑5 (MQ=255)
               gcGGGGGGTAGGGGAGGGCAGGGCTGCCGAGCTGGCCGGGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCAt                                                                                                                                                                                                                                                                                          <  2:650440/300‑1 (MQ=255)
                   gtggtAGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGgtt                                                                                                                                                                                                                                                                                      <  1:188689/301‑1 (MQ=255)
                                 gtgtgctcttccgatctccGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACggg                                                                                                                                                                                                                                                                        <  3:56140/284‑1 (MQ=255)
                                          ggAGCTGGCCGGGCGGGAAGATATTCAGCATGTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAgc                                                                                                                                                                                                                                                               <  2:726649/300‑1 (MQ=255)
                                               tGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTg                                                                                                                                                                                                                                                          <  4:618029/301‑1 (MQ=255)
                                                  ccGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGagatcggaagagcgtcg                                                                                                                                                                                                                                                       >  1:56138/1‑285 (MQ=255)
                                                       cgGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGgcg                                                                                                                                                                                                                                                  >  2:783830/1‑301 (MQ=255)
                                                              atatTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCt                                                                                                                                                                                                                                           <  3:543587/301‑1 (MQ=255)
                                                                atTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCg                                                                                                                                                                                                                                         <  4:784110/301‑1 (MQ=255)
                                                                        atcttACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCa                                                                                                                                                                                                                                 <  4:42199/297‑1 (MQ=255)
                                                                            tACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCTTGCAACGGCGCCTAGCTCGTTATTCCAAgag                                                                                                                                                                                                                             >  2:42188/1‑299 (MQ=255)
                                                                                 gggCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACt                                                                                                                                                                                                                        <  4:344283/301‑1 (MQ=255)
                                                                                       cAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCAt                                                                                                                                                                                                                   <  4:489366/301‑1 (MQ=255)
                                                                                        aaGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATgc                                                                                                                                                                                                                 >  3:353548/1‑301 (MQ=255)
                                                                                                     gCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTg                                                                                                                                                                                                    <  1:221959/301‑1 (MQ=255)
                                                                                                               cGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGcaca                                                                                                                                                                                          >  2:318588/1‑301 (MQ=255)
                                                                                                                  cTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGcc                                                                                                                                                                                       >  1:375323/1‑300 (MQ=255)
                                                                                                                                cagacgtgtgctcttccgatcttAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCAt                                                                                                                                                                         <  4:259862/279‑1 (MQ=255)
                                                                                                                                aGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCAt                                                                                                                                                                         >  3:156243/1‑301 (MQ=255)
                                                                                                                                 gTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATc                                                                                                                                                                        <  1:353547/301‑1 (MQ=255)
                                                                                                                                   cATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATTCGCTCTACCGTTTGCTACAGCACCGCTCTTTCTCTGGCCATCct                                                                                                                                                                      >  2:149399/1‑301 (MQ=255)
                                                                                                                                                    tgTAGGTGATGTTCCATTCAACGCGGCCCAGGCTGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCtt                                                                                                                                                     <  1:609480/300‑1 (MQ=255)
                                                                                                                                                     gTAGGTGATGTTCCATTCAACGCGGCCCAGGCTGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCtta                                                                                                                                                    >  3:609482/1‑300 (MQ=255)
                                                                                                                                                      tAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCAtagatcggcagagcgtcgtgtag                                                                                                                                                   >  2:259804/1‑279 (MQ=255)
                                                                                                                                                       gggTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTg                                                                                                                                                  <  1:156241/300‑1 (MQ=255)
                                                                                                                                                             tGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACa                                                                                                                                            <  4:318661/301‑1 (MQ=255)
                                                                                                                                                                        aCGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGaa                                                                                                                                 <  3:375323/301‑1 (MQ=255)
                                                                                                                                                                             gCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAATGCtt                                                                                                                            >  2:685930/1‑301 (MQ=255)
                                                                                                                                                                                    ccGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGAc                                                                                                                     <  4:686160/301‑1 (MQ=255)
                                                                                                                                                                                     ccacggcgctcttccgatcggAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTGGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACg                                                                                                                    <  2:141649/281‑1 (MQ=255)
                                                                                                                                                                                        gTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGcc                                                                                                                 >  3:114082/1‑300 (MQ=255)
                                                                                                                                                                                        gTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGcc                                                                                                                 >  2:715078/1‑300 (MQ=255)
                                                                                                                                                                                         tgagacgtgtgcacttccgatctcggtgGCGCAGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGAGCTGGAGCGCCTGGAAAGGCGCCTACCTCGTCATGCCAAGGATCACTGCCAGGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCAc                                                                                                                <  4:62447/276‑1 (MQ=255)
                                                                                                                                                                                                         gAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACgagatcggaagagcacacgtc                                                                                                >  4:141678/1‑281 (MQ=255)
                                                                                                                                                                                                              tgcTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCg                                                                                           <  4:360242/300‑1 (MQ=255)
                                                                                                                                                                                                               gcTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCACCCgg                                                                                          >  2:360158/1‑300 (MQ=255)
                                                                                                                                                                                                                cTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTTCTGGCCAGCACCGCTTTTGCCGCCAGCGATGATGAAGGCTTCATCGACGGCAcagatcggaatagcgtcgggtagg                                                                                         >  2:62430/1‑278 (MQ=255)
                                                                                                                                                                                                                  cctttgccctacacgacgctcttccgatctGCTTCGCAGACGTAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAAt                                                                                       <  2:718647/272‑1 (MQ=255)
                                                                                                                                                                                                                    ggCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAAttt                                                                                     <  1:114079/301‑1 (MQ=255)
                                                                                                                                                                                                                                     ttGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACt                                                                    >  3:39514/1‑301 (MQ=255)
                                                                                                                                                                                                                                                gCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAAtagatcggaagagcacacgtctgaactccag                                                         >  4:718892/1‑271 (MQ=255)
                                                                                                                                                                                                                                                gCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGc                                                         >  4:225460/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  catAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACCATCAATccc                                       >  2:671504/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                       tgtTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCTCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAATGCTTCATCGACGGCACGCGGGTCCCTGTTCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCTCAACGTCTACTTTACCCTCAATCCCaatac                                  >  2:267464/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                         ttGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATg                                <  4:715323/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                 gttgttCAATTGCTTGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCgg                        >  1:351536/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                               tGAATAACAATGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCg          >  1:250058/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                       aaTGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATggg  >  2:412073/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                       aaTGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATggg  <  4:149437/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                       aaTGGTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATggg  <  4:671726/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                        
CATCCTGGAAAACGGCCGGGTGGTAGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATGGG  >  AE015451/3137476‑3138072

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