Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,205,402 C→A 8.7% R88S (CGT→AGT)  pinE → e14 prophage; site‑specific DNA recombinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,205,4020CA 8.7% 142.3 / 2.5 46R88S (CGT→AGT) pinEe14 prophage; site‑specific DNA recombinase
Reads supporting (aligned to +/- strand):  ref base C (18/24);  new base A (2/2);  total (20/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 = 7.67e-01

GGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGCACAGGGACGAATTGGTGGACGTCGTCCCAAACTTACACCAGAACAATGGGCACAAGCTGGACGATTAATTGCAGCAGGAACTCCTCGCCAGAAGGTGGCGATTATCTATGATGTTGGTGTGTCAACTTTGTATAAGAG  >  CP009273/1205130‑1205673
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
ggggATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGt                                                                                                                                                                                                                                                                               >  1:187787/1‑275 (MQ=255)
            ggattATTGTCTATGTACGCGTAGCAACAAATGACCAAAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGAGTCGAGCTGATTTNTGAAGACAAGATAAGCGGTACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGa                                                                                                                                                                                                                                                                    <  1:539029/271‑1 (MQ=255)
             gCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGAtt                                                                                                                                                                                                                                                                  >  1:89449/1‑275 (MQ=255)
               ttATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGAtt                                                                                                                                                                                                                                                                  <  1:152328/273‑1 (MQ=255)
               ctagtGGCTATTTACGCGTATCAACCAATGACCAGAAAACAGACCTACAACGTAATGCGCTGAAGTGTGCAGGATGCGGGCGGATTTTTGAAGACAAGATAAGCGGAACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGAtt                                                                                                                                                                                                                                                                  <  2:725740/269‑1 (MQ=255)
                                               acccaaagctaTCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCttttt                                                                                                                                                                                                                                <  2:89448/265‑1 (MQ=255)
                                                                  gttctgggaggAACTGGGCAGGGTGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGcc                                                                                                                                                                                                             <  2:187785/265‑1 (MQ=255)
                                                                  gTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGcc                                                                                                                                                                                                             <  1:16396/275‑1 (MQ=255)
                                                                              aCTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGNAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAAt                                                                                                                                                                                                  <  1:678487/274‑1 (MQ=255)
                                                                                       gATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGTGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGCCGGCTTCAATTGATACCAGCACACCAATGGGCCCCTTTTTATTTCATGGGATGGGTGCCCTTGCTGAAATGGAGCGTGa                                                                                                                                                                                         >  2:678473/1‑274 (MQ=255)
                                                                                          gggAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGTCTGAAAAAACTGCTCAGGACAGTATCTGCAGGTGACACTCTGGATGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACt                                                                                                                                                                                      <  1:406929/272‑1 (MQ=255)
                                                                                            gAGCTGACTCTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGCTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGGGAACGAGGCATCAACTTGCGTAGTCTGGAGGATGCAATTGAGACCAGCACAAAAAAGGGACGATTTTGCGTTCATGTGATTGGTGCCTTGGCTGAAAGGGAGCGGGACTTGAt                                                                                                                                                                                   >  2:721081/1‑275 (MQ=255)
                                                                                                                        gCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGAGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTc                                                                                                                                                                                                                                                                           >  1:369940/1‑159 (MQ=255)
                                                                                                                            aaaGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGAGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGAc                                                                                                                                                                                                                                                                       <  2:369933/159‑1 (MQ=255)
                                                                                                                                 ccGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTCACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTATCATGTGATGGGTg                                                                                                                                                                                                               >  1:711284/1‑210 (MQ=255)
                                                                                                                                     aaGGCCGGGACTGATAAAACTGTTCATGTCATTATCGGCAGGTGACACTCTGGGTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTCGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGca                                                                                                                                           <  1:387994/274‑1 (MQ=255)
                                                                                                                                     aaGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCt                                                                                                                                                                                                           <  2:711270/210‑1 (MQ=255)
                                                                                                                                          cGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAg                                                                                                                                                                                              >  1:582618/1‑218 (MQ=255)
                                                                                                                                              aCTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTg                                                                                                                                                                                          <  2:582605/218‑1 (MQ=255)
                                                                                                                                                          gCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATAcc                                                                                                                                                                                                                                                       >  1:23527/1‑145 (MQ=255)
                                                                                                                                                              aGGACATTATCGGCAGGTGACACTCGGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGca                                                                                                                                                                                                                                                   <  2:23526/145‑1 (MQ=255)
                                                                                                                                                                          gCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTAGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTg                                                                                                                                                                                          >  1:304689/1‑190 (MQ=255)
                                                                                                                                                                          gCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTAGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTg                                                                                                                                                                                          >  1:99985/1‑190 (MQ=255)
                                                                                                                                                                             ggTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGCTGGGTGCCCTGACTGAAATTGAGCGTGCACTGAttgt                                                                                                                                                                                >  2:258593/1‑197 (MQ=255)
                                                                                                                                                                              gTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTAGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACt                                                                                                                                                                                      <  2:304684/190‑1 (MQ=255)
                                                                                                                                                                              gTGACAATCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTAGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACt                                                                                                                                                                                      <  2:99984/190‑1 (MQ=255)
                                                                                                                                                                                 acacTCTGGTTGTCGGGAAGCTGGAGCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTgaa                                                                                                                                                                            <  1:258596/197‑1 (MQ=255)
                                                                                                                                                                                        ggTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGATCGTGAACTGCTTGTTGTACGACCACAATCTGGACTGGAAACTGCTCGTGCACAGGGACGAATTGTTCGAcctcg                                                                                                                    >  1:359890/1‑246 (MQ=255)
                                                                                                                                                                                        ggTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCg                                                                                                                                               >  2:649988/1‑219 (MQ=255)
                                                                                                                                                                                         gTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTtttctttc                                                                                                                                                                                                                           >  1:321151/1‑142 (MQ=255)
                                                                                                                                                                                            gTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCTAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGca                                                                                                                                           <  1:650002/219‑1 (MQ=255)
                                                                                                                                                                                             tCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCAtgt                                                                                                                                                                                                                       <  2:321145/142‑1 (MQ=255)
                                                                                                                                                                                                                       agtaTGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGAGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAATTGATTGTTGAACGCACAAAAGCTGGCATGGAAACTGCTCGTGCAAAGGGAAGAATTGGTGGACGGAGTCACAAACTTACACAAGAAAAATGGGCCAAAGCTGGCCGATTAAATTAAGCAAGaaacg                                                         >  2:655460/1‑271 (MQ=255)
                                                                                                                                                                                                                         taTGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGCACAGGGACGAATTGGTGGACGTCGTCCCAAACTTACACCAGAACAATGGGCACAAGCTGGACGATTAATTGCAGCAGGAActcct                                                       <  1:721095/274‑1 (MQ=255)
                                                                                                                                                                                                                                   cTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATa                                                                                                                                                                                                                                                         >  2:494972/1‑70 (MQ=255)
                                                                                                                                                                                                                                   cTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGc                                                                                                                                            >  1:253095/1‑179 (MQ=255)
                                                                                                                                                                                                                                       tCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCGACTTTCGTAGTCTGACGGATTCAATTGATACCAg                                                                                                                                                                                                                                                     <  1:494981/70‑1 (MQ=255)
                                                                                                                                                                                                                                       tCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGCACAg                                                                                                                                        <  2:253092/179‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGaa                                                                                                                                                                                        >  1:590679/1‑126 (MQ=255)
                                                                                                                                                                                                                                            cTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCg                                                                                                                                               >  2:394899/1‑167 (MQ=255)
                                                                                                                                                                                                                                                tggAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGa                                                                                                                                                                                    <  2:590666/126‑1 (MQ=255)
                                                                                                                                                                                                                                                tggAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGca                                                                                                                                           <  1:394906/167‑1 (MQ=255)
                                                                                                                                                                                                                                                      agTTGCGCGAATGAGGCATCAACTTTCGTATTCTGACGTATTCTATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCTTGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGCACAGGGACGAATTGGTGGACGTCGTCCCAAACTTACACCAGAACAATTGGCACAAGCTGGACGATTAATTGCAGCAGGAACTCCTCGCCAGAAGGTGGCGATTATCTATGATGtt                         <  1:44639/275‑1 (MQ=255)
                                                                                                                                                                                                                                                         tGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCg                                                                                                                                               >  1:375935/1‑154 (MQ=255)
                                                                                                                                                                                                                                                             cgAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGca                                                                                                                                           <  2:375928/154‑1 (MQ=255)
                                                                                                                                                                                                                                                                              ttCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTTGAAACTGCTCGTGCACAGGGACGAATTGGTGGACGTCGTCCCAAACTTACACCAGAACAATGGGCACAAGCTGGACGATTAATTGCAGCAGGAACTCCTCGCCAGAAGGTGGCGATTATCTATGATGTTGGTGTGTCAACTTTGTATAagag  <  1:655474/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
GGGGATTTAACATGCTTATTGGCTATGTACGCGTATCAACAAATGACCAGAACACAGATCTACAACGTAATGCGCTGAACTGTGCAGGATGCGAGCTGATTTTTGAAGACAAGATAAGCGGCACAAAGTCCGAAAGGCCGGGACTGAAAAAACTGCTCAGGACATTATCGGCAGGTGACACTCTGGTTGTCTGGAAGCTGGATCGGCTGGGGCGTAGTATGCGGCATCTTGTCGTGCTGGTGGAGGAGTTGCGCGAACGAGGCATCAACTTTCGTAGTCTGACGGATTCAATTGATACCAGCACACCAATGGGACGCTTTTTCTTTCATGTGATGGGTGCCCTGGCTGAAATGGAGCGTGAACTGATTGTTGAACGAACAAAAGCTGGACTGGAAACTGCTCGTGCACAGGGACGAATTGGTGGACGTCGTCCCAAACTTACACCAGAACAATGGGCACAAGCTGGACGATTAATTGCAGCAGGAACTCCTCGCCAGAAGGTGGCGATTATCTATGATGTTGGTGTGTCAACTTTGTATAAGAG  >  CP009273/1205130‑1205673

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