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.96.8% 133.5 / ‑1.6 31intergenic (+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 (1/0);  new base . (16/14);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

AGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATGGG  >  AE015451/3137500‑3138072
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                        
gggTGAGGGCAGTGCTGCCGGGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTgg                                                                                                                                                                                                                                                                                 <  2:703378/300‑1 (MQ=255)
 ggTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGa                                                                                                                                                                                                                                                                                >  1:725826/1‑301 (MQ=255)
       ggCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACg                                                                                                                                                                                                                                                                          <  1:835937/301‑1 (MQ=255)
              ctgcCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTg                                                                                                                                                                                                                                                                   >  1:361124/1‑301 (MQ=255)
                                gggAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGgcgc                                                                                                                                                                                                                                                 >  1:207683/1‑301 (MQ=255)
                                         tgtcgcgtttatactggggCGGGAAGGGCAGACGGGGCGGTGGCGACGGGTGTGGAGCGGCGCAGTAATACCGCCTTTTACTGAGTAGGTGATGTTCCATTGAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGGCGAAAGCGTTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGt                                                                                                                                                                                                                                        <  2:419543/285‑1 (MQ=255)
                                                                            ggCGGTGGCGACGGCTTTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGTTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGttt                                                                                                                                                                                                     <  2:725573/300‑1 (MQ=255)
                                                                                                      gcAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTTCTACAGAACAGCTCTTTCTCTGGcc                                                                                                                                                                           >  2:84927/1‑301 (MQ=255)
                                                                                                                       tACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGTGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGccc                                                                                                                                                          >  1:765455/1‑301 (MQ=255)
                                                                                                                        acttaGTATGTGCTGTTCCATTCAACGCGGGCCATGCCGCGTTTTTTTGGCCCTGACGAAAGCGCTGTGGCGCCTGCAACAGAAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGcccc                                                                                                                                                         <  2:493981/297‑1 (MQ=255)
                                                                                                                           gAGTAGGTGATGTTCCATTCAACGCGGCCCATGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCt                                                                                                                                                      <  2:735091/301‑1 (MQ=255)
                                                                                                                           gAGTAGGGGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCt                                                                                                                                                      <  2:560667/301‑1 (MQ=255)
                                                                                                                               aGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCttttt                                                                                                                                                  >  2:721714/1‑300 (MQ=255)
                                                                                                                                                        cAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCCCTGGCCCTCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGCTGAAGGCtttat                                                                                                                         >  2:535760/1‑298 (MQ=255)
                                                                                                                                                                     tttCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCgg                                                                                                            >  1:872892/1‑301 (MQ=255)
                                                                                                                                                                         ccccTTACGAAAGGGCTGCGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTTGCTGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TGTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTACCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTcc                                                                                                        <  2:361011/300‑1 (MQ=255)
                                                                                                                                                                           ccTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCAt                                                                                                      <  2:872584/301‑1 (MQ=255)
                                                                                                                                                                           ccTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCAt                                                                                                      <  1:721967/301‑1 (MQ=255)
                                                                                                                                                                                          gtgGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAAt                                                                                       <  1:535946/301‑1 (MQ=255)
                                                                                                                                                                                                  ggCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTTGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATttcttctt                                                                               <  2:207614/301‑1 (MQ=255)
                                                                                                                                                                                                      aCAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGc                                                                           >  2:287669/1‑301 (MQ=255)
                                                                                                                                                                                                        agcagTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCa                                                                         >  2:273395/1‑301 (MQ=255)
                                                                                                                                                                                                                   aGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCa                                                              >  1:788697/1‑301 (MQ=255)
                                                                                                                                                                                                                             gCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACtt                                                    <  1:273485/301‑1 (MQ=255)
                                                                                                                                                                                                                                                  tGAGGAAGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGCAAAGGCGCCTACCTCGTCATCCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACCGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAACGCTTCATCGACGGCACCCGGGTCCAGGTGCTCAACCGCAATTTCTTCTTCACCCAGTACTATCTCAACGCCTACTTTCCACTCCATCCCaatcattc                               >  2:595106/1‑296 (MQ=255)
                                                                                                                                                                                                                                                       aaGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCTGCTCCATGTGCTCAACCGCAATTTCTTCTTCCGCCAGGACTATCGCAACGGCGACTTTAACCTCCATCCCaatactgtcgacg                          >  2:764803/1‑292 (MQ=255)
                                                                                                                                                                                                                                                       aaGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAgg                          >  1:127303/1‑300 (MQ=255)
                                                                                                                                                                                                                                                        aGTTGTTCAATTGCTTGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCg                         <  1:62665/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ttGAATAACAATGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCCTCCTGCTGGCCAGCCCCGCTTTTGCCGCCAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTAACCTCAATCCCACTACTGGCGAGCGTCAGAGcacgcatc           >  2:49110/1‑294 (MQ=255)
                                                                                                                                                                                                                                                                             aCAATGGTTTTTTTTT‑TTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCCCGCGGGTCCCTGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGCTTCCTTTACCCTCAATCCCAATAATGGccatacgcaaacattgcctcgttcatt    >  2:194818/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                               aaTGG‑TTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTGAACCGCAATTTCTTCTTCAGCCAGGACTATGGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCAGCGCGAATggg  >  1:531002/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                                                                        
AGGTGAGGGCAGTGCTGCCGAGCTGGCCGCGCGGGAAGATATTCAGCATTTCTACCTGGGCGGCAAGGGCTGACTGAGCGGTGGCGACGGCTGTGGAGCGGCGCAGTCATACCGCCTTTTACTGAGTAGGTGATGTTCCATTCAACGCGGCCCAGGCCGCGTTTTTTTCGCCCTGACGAAAGCGCTGTGGCGCCGGCAACAGCAGTTGTGCAGATTGCTTCGCAGACGGAGCATCATAGTGTTGAGGAAGTTGTTCAATTGCTTGAATAACAATGGTTTTTTTTTGTTGGCATGGTTGTTGGAACAACGGGATCTGCAGCGCCTGGAAAGGCGCCTACCTCGTCATTCCAAGGATCACTGCCATGCGCTCTACCGTTTGCTACAGAACAGCTCTTTCTCTGGCCATCCTGCTGGCCAGCCCCGCTTTTGCCGACAGCGATGATGAAGGCTTCATCGACGGCACGCGGGTCCATGTGCTCAACCGCAATTTCTTCTTCAGCCAGGACTATCGCAACGGCGACTTTACACTCAATCCCAATAATGGCGAGCGGCAGAGCCTGCATCGTGAATGGG  >  AE015451/3137500‑3138072

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑