Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,444 C→T 47.9% R36R (AGG→AGA)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4440CT47.9% 111.7 / inf 342R36R (AGG→AGA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (84/92);  new base T (84/79);  total (171/171)
Fisher's exact test for biased strand distribution p-value = 5.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

Only 100 of 347 total aligned reads displayed.

TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTA  >  CP000730/36304‑36583
                                                                                                                                             |                                                                                                                                            
   tttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                                                           <  2:407037/140‑1 (MQ=39)
   tttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                           <  2:416068/140‑1 (MQ=11)
   aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                                                           <  2:475625/141‑1 (MQ=255)
    ttttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTa                                                                                                                                          <  1:951470/141‑1 (MQ=18)
    ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTa                                                                                                                                          <  2:750102/141‑1 (MQ=11)
    ttttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTa                                                                                                                                          <  1:126474/141‑1 (MQ=255)
     tttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCGAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTaa                                                                                                                                         <  2:699041/141‑1 (MQ=255)
     tttGATGTACCTCTCATCAATACGCCACTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGCGCATATTCTTGAACCCA‑ACCGTAGACCGTTGAATTATGAACGTTTACACCACGTTCCCTTaa                                                                                                                                         <  2:808891/141‑1 (MQ=14)
     tttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTaa                                                                                                                                         <  2:21935/141‑1 (MQ=255)
      ttGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAat                                                                                                                                        <  2:784533/141‑1 (MQ=11)
       tGATTTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAata                                                                                                                                       <  2:127791/141‑1 (MQ=18)
        ggcaaggatctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCCTAatat                                                                                                                                      >  2:275609/9‑141 (MQ=11)
        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGTTCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                                      >  1:619626/5‑141 (MQ=11)
        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                                                      >  2:334003/5‑141 (MQ=255)
        gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTAGGGCATATTCTTGAACCCA‑ACGTTAGATCGTTGAATGATGAACGTTTACAACACGTTCTCTTAatat                                                                                                                                      >  1:156234/5‑141 (MQ=11)
         atatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACTACGTTCTCTTAATAtt                                                                                                                                     >  1:879967/4‑141 (MQ=11)
         aTGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAtt                                                                                                                                     >  1:559107/1‑141 (MQ=21)
            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCCCGTTCACTTAATATTTCa                                                                                                                                  >  2:622053/1‑141 (MQ=11)
            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAAGAAATACAAT‑CATCTTGGAATGATGAACGATTGCACCACGTTCCCTTAATATTTCa                                                                                                                                  >  1:1037395/1‑141 (MQ=11)
            tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                                                  >  1:976893/1‑141 (MQ=40)
            aaCGTCTCATCAATACGCCATCTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATACTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCa                                                                                                                                  <  1:567834/140‑1 (MQ=11)
              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTTTTTTTTTTTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTCCACCTCTTTCTCTTAATATTTCAGa                                                                                                                                >  2:31529/1‑141 (MQ=11)
              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                                                >  1:232992/1‑141 (MQ=255)
              cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGa                                                                                                                                >  2:972092/1‑141 (MQ=21)
                gctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGata                                                                                                                              <  2:1060084/140‑1 (MQ=40)
                 ttcatcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                             <  2:194425/140‑1 (MQ=255)
                 ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                                                             <  2:548382/141‑1 (MQ=255)
                    ttcaATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  1:795783/140‑1 (MQ=255)
                    cacaATACGCCATGTTTAATAAGCTTTCTTTTGCTTTTTCTTCCAAATTTGATATAAAATTTGGGCATATTCTTGAACCCA‑CCCGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  2:919179/139‑1 (MQ=11)
                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATc                                                                                                                          <  1:382823/141‑1 (MQ=255)
                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                                                          <  2:442966/141‑1 (MQ=18)
                    atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATc                                                                                                                          >  2:843246/1‑141 (MQ=255)
                    atcaATACGCCAATTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAGTGGAGATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGATCCCCTAATATTTCAGatatcta                                                                                                                          >  2:358341/1‑138 (MQ=11)
                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCCTATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCGGATATATCa                                                                                                                         >  2:186123/1‑141 (MQ=11)
                     tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCa                                                                                                                         >  1:984379/1‑141 (MQ=21)
                     acaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                                                         <  2:869601/140‑1 (MQ=255)
                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAACTTGGGGCATCTTCTTGAACCCA‑AAGGTAGCTCGTTGAATGTTGGAGCTTTACACTGCGTTGTTTTAATATTTCTGatttataag                                                                                                                        >  2:466540/1‑138 (MQ=2)
                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                                                        >  2:144167/1‑141 (MQ=255)
                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATAAc                                                                                                                        >  2:943744/1‑141 (MQ=255)
                      caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAc                                                                                                                        <  2:247809/141‑1 (MQ=39)
                       aaTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGGTCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACg                                                                                                                       >  2:937515/1‑141 (MQ=11)
                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                                                      >  1:420394/1‑141 (MQ=255)
                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACATTTACACCACGTTCTCTTAATGTTTCAGATATATCACGa                                                                                                                      >  1:112365/1‑141 (MQ=11)
                        aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGa                                                                                                                      <  2:849154/141‑1 (MQ=18)
                         tACGCCATTTGTAATAAGCTTTTTTTTTTTTTTTTTTCCAAATTTGATATAAAATTGGGGCATATTTTTGAACGCA‑ACGGTAGATCGTTGAATGCTGAACGTGTACACCACGTTCTCTTAATATTTAAGATATGTTACGAt                                                                                                                     >  2:800466/1‑141 (MQ=11)
                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                                                     <  1:446464/141‑1 (MQ=255)
                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGAt                                                                                                                     >  1:259792/1‑141 (MQ=255)
                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGAt                                                                                                                     >  1:652192/1‑141 (MQ=21)
                         tACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGATTACACCACGTTCCCTTAATATTTCAGATATATCACGAt                                                                                                                     >  1:474862/1‑141 (MQ=18)
                          cccactcttgaactaaGCTTTTTCATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATa                                                                                                                    <  2:727628/128‑1 (MQ=11)
                          aCGCCATTTGTAATAAGCTTTTTTATGCTTTTACTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTGCTCTTAATATTTCAGATGTATCAGGATa                                                                                                                    >  1:434471/1‑141 (MQ=11)
                          aCGCAATTTGTACTAAGCTCTTTTATGTTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCATTAATATTTCAGATATATCACGATa                                                                                                                    <  1:677009/141‑1 (MQ=11)
                           cGCCATTTGTAATAAGCTTTTTTTTTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGATATATCCCGATaa                                                                                                                   >  2:654033/1‑141 (MQ=11)
                           cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATaa                                                                                                                   <  1:854560/141‑1 (MQ=18)
                            gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                                                  <  1:751357/141‑3 (MQ=255)
                            gCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAAc                                                                                                                  <  2:278061/141‑1 (MQ=39)
                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATGTTCTTGAACCCA‑ACGGTAGACCGTTGAATGGTGAACGTTTACACCCCGTTCCCTTACTATTTCAGATATATGACGATAACt                                                                                                                 >  2:582954/1‑141 (MQ=11)
                             ccATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCt                                                                                                                 >  1:908149/1‑141 (MQ=255)
                              cATTTGTCATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCAGGTCCCCTTACTGTTTCTGATATATCACGATCACTc                                                                                                                >  2:596130/1‑141 (MQ=11)
                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTc                                                                                                                >  2:563539/1‑141 (MQ=255)
                              cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                                                <  1:54976/141‑1 (MQ=25)
                               aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCa                                                                                                               <  1:846625/141‑1 (MQ=255)
                                tttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                              <  2:738573/141‑1 (MQ=11)
                                tctGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                                                              <  1:281844/139‑1 (MQ=11)
                                 ttttaataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             <  1:247912/138‑1 (MQ=11)
                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             <  1:413385/141‑1 (MQ=255)
                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAAt                                                                                                             <  2:451898/141‑1 (MQ=11)
                                 ttGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGACCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAAt                                                                                                             <  1:703768/141‑1 (MQ=18)
                                  tGTAATAAGGTTTATTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  1:49450/141‑1 (MQ=11)
                                  tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATg                                                                                                            <  1:1008278/141‑1 (MQ=255)
                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                           <  1:1001540/141‑2 (MQ=11)
                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgc                                                                                                           <  1:105569/141‑2 (MQ=11)
                                   gTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATAAATCACGATAACtaccaga                                                                                                           >  1:392533/1‑135 (MQ=255)
                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCATTgca                                                                                                          >  1:303607/1‑139 (MQ=11)
                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAAAGTTTACACAACGATCTCTTAATATTTCAGATATATCACGATAGCTCaaaaca                                                                                                          >  2:897913/1‑137 (MQ=11)
                                    taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTAGGTCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATgca                                                                                                          <  1:344916/141‑3 (MQ=11)
                                    taataaGCTTTTTTATGCTTTTTCCTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTATTATTTTACATATATCACGATAACTCAATgca                                                                                                          <  1:514735/141‑3 (MQ=11)
                                     aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTGTTAATATTACGGATATATCCAGATAGCTCGATgcaa                                                                                                         >  2:1048947/1‑138 (MQ=11)
                                     aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATGCCTCAATgcaa                                                                                                         >  1:20341/1‑138 (MQ=255)
                                            ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAg                                                                                                 >  1:661431/1‑141 (MQ=35)
                                            ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATCACGATAACTCAATG‑TATATCTTaa                                                                                                 >  1:321067/1‑140 (MQ=25)
                                                  aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAACTCAATG‑TATATCTTAGAtagta                                                                                           >  1:191467/1‑141 (MQ=25)
                                                       ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCaa                                                                                      <  2:1004698/141‑1 (MQ=35)
                                                        tttCCTCCAAATTTGATACAAAATTGGTGCTTATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGccgtt                                                                                     <  1:35865/141‑4 (MQ=255)
                                                             tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTTGCCAACGGCTa                                                                                >  1:586296/1‑141 (MQ=25)
                                                             tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTa                                                                                >  1:1025172/1‑141 (MQ=35)
                                                                  aTTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGTTAGACGGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                           <  2:515535/141‑1 (MQ=21)
                                                                    ttGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                         <  2:506800/141‑1 (MQ=35)
                                                                       aTACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                      >  1:716115/1‑141 (MQ=35)
                                                                              aTTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAAGG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                               <  1:728073/141‑1 (MQ=17)
                                                                               ttcgggaagagtCTTTAACCCA‑ACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                              <  1:87006/130‑1 (MQ=11)
                                                                                tGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                             <  1:147010/141‑1 (MQ=32)
                                                                                      atatTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGtacaaacg                                                       >  2:630935/1‑134 (MQ=255)
                                                                                        atTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                     <  1:773464/141‑1 (MQ=32)
                                                                                                              ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                <  1:66892/141‑1 (MQ=32)
                                                                                                                    aaTGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                          <  2:692814/141‑1 (MQ=17)
                                                                                                                          gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                    <  2:686530/141‑1 (MQ=32)
                                                                                                                            aCGTTCACACCACGTTCGCTTAATATTACTGATATATCA‑GCCCAGTTAACGTTAGAGGTTAGTTAATATCGATCTGCTACACTGAGTACCTCCTTGTTAAATTGTTTATATCTGAAATATTTCATACAGAAAACTCCtttt                  >  2:382093/1‑141 (MQ=2)
                                                                                                                                     ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaag         >  2:645991/1‑140 (MQ=17)
                                                                                                                                          ttgccTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc    <  2:296580/138‑1 (MQ=17)
                                                                                                                                             |                                                                                                                                            
TTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCA‑ACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATG‑TATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTA  >  CP000730/36304‑36583

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑