Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,479 T→C 22.5% T25A (ACA→GCA)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC22.5% 221.0 / 56.7 133T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (35/65);  new base C (14/16);  total (49/84)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

ATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAT  >  CP000730/36342‑36619
                                                                                                                                         |                                                                                                                                            
ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGAGAGAAAATTGGGGCATTTTGTAGAACCCGACGGTAGATCGGTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGcata                                                                                                                                           >  2:232862/1‑141 (MQ=11)
   aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                        <  2:842324/141‑1 (MQ=38)
    gCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCtt                                                                                                                                       <  2:301953/141‑1 (MQ=37)
    gCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCtt                                                                                                                                       <  2:110127/141‑1 (MQ=38)
      ttttttATGCTTTTTCTTCCAAATTTGATACAAACTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCttga                                                                                                                                     >  2:933238/1‑139 (MQ=25)
          tttttttttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTTCCTTATTATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtaa                                                                                                                                 >  2:149320/7‑140 (MQ=255)
          ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATTTTAGAtag                                                                                                                                 >  1:161980/1‑141 (MQ=21)
          ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                                 >  1:347523/1‑141 (MQ=35)
           tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCGTATACTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtggt                                                                                                                                >  1:531312/1‑141 (MQ=255)
           tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATGCCGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  2:527862/1‑141 (MQ=21)
           tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                                >  1:1053774/1‑141 (MQ=35)
            aTGCTTCTTCTTCCAACTTTGATACAAAATTGGGGCACATTCTTGCACCCAAAGGTTGATCGTTGAATGTTGTACGTTTCCACGCCATTCTCTTATTATTCTCGATATATTACGATCGATCAATGTCTCTCCTAGAtagta                                                                                                                               >  2:1029873/1‑141 (MQ=255)
             tGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGACAACTCACTGTATATCTTAGAtagtag                                                                                                                              >  1:380345/1‑141 (MQ=21)
              gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                             <  1:423579/141‑1 (MQ=21)
              gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                             <  1:509463/141‑1 (MQ=35)
                 ttttCTTCCAAATTTGATGCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCaa                                                                                                                          <  2:151450/141‑1 (MQ=25)
                 ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCTTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCaa                                                                                                                          <  1:252949/141‑1 (MQ=25)
                 ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCaa                                                                                                                          <  2:844571/141‑1 (MQ=35)
                 ttttCTTCCAAATTTGATACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCaa                                                                                                                          <  1:355295/141‑1 (MQ=25)
                     cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                      <  1:304546/141‑1 (MQ=35)
                     cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                      <  1:992564/141‑1 (MQ=35)
                       tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                    <  2:767944/141‑1 (MQ=35)
                           aaTTTGATACAAAATTGGGGCATATTCTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGccgttggctactat                                                                                                                <  2:472622/141‑13 (MQ=11)
                           aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                                <  2:900364/141‑1 (MQ=35)
                                   aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAt                                                                                                        >  1:1109001/1‑141 (MQ=33)
                                    caaaaatcaggCATATTCCTGAACCCAACCGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATCTCAAATATATCACGATAACTCAATTTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATc                                                                                                       <  1:233011/132‑1 (MQ=11)
                                     aaaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTGATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                                      <  2:839647/141‑1 (MQ=17)
                                     aaaaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATcc                                                                                                      <  2:78107/141‑1 (MQ=32)
                                        aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATACATCCTTGt                                                                                                  <  1:325712/141‑1 (MQ=17)
                                         ttGGCGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                  <  1:389810/141‑1 (MQ=17)
                                         tgggggTATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                  <  1:838133/139‑1 (MQ=17)
                                          tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                                 >  2:345057/1‑141 (MQ=32)
                                           ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCTCGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                                <  1:639085/141‑1 (MQ=17)
                                              gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                             >  1:328832/1‑141 (MQ=32)
                                              gCATATTCTTGAACCCAACGGTAGACCGGTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATACCTTAGATAGTAGCCAACGGCTCCAGTGATAACATCCTTGTTaaaa                                                                                             >  2:892995/1‑140 (MQ=12)
                                                   ttCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACCGTGATAACATCCTTGTTAAATTGttt                                                                                        >  1:62142/1‑141 (MQ=17)
                                                         gaCCCAACTGTAGACCGTTGAAGGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                  <  1:227666/140‑1 (MQ=17)
                                                         aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                  <  1:730707/141‑1 (MQ=32)
                                                         aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAGCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                  <  2:291439/141‑1 (MQ=17)
                                                          aCCCAACGGTAGATAGTTGTATGATCAACGTTTACACCACGTTCTATTAATATTTCAGATATATCACAATAGATCAATGTATATCTTAGTTTGTAGCTAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTTTATCTg                                                                                 <  2:752/141‑1 (MQ=11)
                                                          aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATTTCTg                                                                                 <  2:669738/141‑1 (MQ=17)
                                                          aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCATTATATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                                 >  1:1088459/1‑141 (MQ=12)
                                                           cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                                >  1:1030341/1‑141 (MQ=32)
                                                             cAACGGTAGACCGTTGAATGATGAATGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGgaa                                                                              <  2:459244/141‑1 (MQ=17)
                                                             cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGCTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                              <  1:736489/141‑1 (MQ=17)
                                                             cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                              <  1:793892/141‑1 (MQ=17)
                                                             cAACGGTAGACCAGTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAACCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                              <  1:811360/141‑1 (MQ=12)
                                                                  gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTACATTGTTGATATCTGAAATAGtt                                                                         >  1:760133/1‑141 (MQ=17)
                                                                       ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                    >  2:1090532/1‑141 (MQ=32)
                                                                        cGTTGAATGATGAACGTTTACCCCTCGTTCCCTTAATATTTCAGATATAACACCACACGGCAGCGTCTATCTCACATAGTAGCCCTGGGTCACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACa                                                                   <  2:501892/141‑1 (MQ=2)
                                                                        cGTTGAAGAATGAACGTTTACACCACGTTTCCTTAATATTTCAGATAATTCACGATAACTCAATGTATAGCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAcg                                                                   <  2:1019382/141‑2 (MQ=11)
                                                                          ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                                 <  2:536629/141‑1 (MQ=32)
                                                                               tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                            <  2:793340/141‑1 (MQ=32)
                                                                                 atgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTcc                                                          >  2:292963/1‑141 (MQ=32)
                                                                                       gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                    <  1:39121/141‑1 (MQ=32)
                                                                                        tttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                   <  1:181807/141‑1 (MQ=32)
                                                                                         ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                  <  1:933500/141‑1 (MQ=32)
                                                                                         taacacCACGTTCCCTTAATAGTTCAGATATATCACGATAACCCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGCCTCCTTTTTGtt                                                  <  2:158942/139‑1 (MQ=255)
                                                                                          tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                                 <  1:70004/141‑1 (MQ=32)
                                                                                                      ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                     <  1:727868/141‑1 (MQ=32)
                                                                                                      ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                     <  1:579120/141‑1 (MQ=32)
                                                                                                       cTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtata                                    <  2:949799/141‑1 (MQ=32)
                                                                                                         tAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                                  <  1:594835/141‑1 (MQ=34)
                                                                                                         tAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                                  <  1:323422/141‑1 (MQ=33)
                                                                                                           atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                                <  2:230664/141‑1 (MQ=35)
                                                                                                            tatTTCATATATATCACGATAACTCAATGTTTATACTAGATAGTACTCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTc                               <  1:573372/141‑1 (MQ=11)
                                                                                                            aatTTCAGATAAATCACGATAACTCAATGAATAACAAAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTc                               <  1:7589/140‑1 (MQ=11)
                                                                                                               ttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTGCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAAc                            >  1:691350/1‑141 (MQ=25)
                                                                                                                 cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                          >  1:1110529/1‑141 (MQ=34)
                                                                                                                 cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                          >  1:285069/1‑141 (MQ=35)
                                                                                                                    gtatatTACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGc                       <  1:467693/140‑1 (MQ=25)
                                                                                                                       tatCACTATAACGCAATGGATATGTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCCTACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                    <  2:428043/141‑1 (MQ=11)
                                                                                                                       tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                    >  2:1025791/1‑141 (MQ=35)
                                                                                                                       tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                    >  1:659772/1‑141 (MQ=35)
                                                                                                                       tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTAAACTTTGCAAc                    >  1:1046244/1‑141 (MQ=25)
                                                                                                                       tagcACTATAACTCAATGCATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGACATAGTTCATACAGAAGACTCCTTTTTGTTAAAATCATACTATAAATTCAACTTTGCAAc                    <  2:937024/138‑1 (MQ=1)
                                                                                                                            cgaaatcTCAATGCATATCTTAGATAGTACGCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc               <  2:338911/135‑1 (MQ=11)
                                                                                                                            cGATTACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc               <  1:905993/141‑1 (MQ=25)
                                                                                                                            cGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAc               <  1:594818/141‑1 (MQ=35)
                                                                                                                              atagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg             <  1:1001994/137‑1 (MQ=9)
                                                                                                                               tagctaaaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTTTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  2:672962/135‑1 (MQ=11)
                                                                                                                               tagctaaaTGCATATCTTAACTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATCAATTCCCCTTTGCAACAGAACCGt            <  2:1097549/135‑1 (MQ=11)
                                                                                                                               tagcTCACTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTAAAAGTAATTTTATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  1:804366/138‑1 (MQ=11)
                                                                                                                               tagcTCAATGCATATCTTAGCTCGTCGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTATAATTATATTATAAATTCAACTTTGCAAAATAACCGt            >  2:237151/4‑141 (MQ=11)
                                                                                                                               tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  1:185295/138‑1 (MQ=9)
                                                                                                                               tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  1:419738/138‑1 (MQ=12)
                                                                                                                               tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt            <  2:213008/138‑1 (MQ=12)
                                                                                                                                agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta           >  1:576448/3‑141 (MQ=12)
                                                                                                                                 gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCTACTTTGCAACAGAACCGtat          >  1:959992/2‑141 (MQ=12)
                                                                                                                                 gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat          >  1:596919/2‑141 (MQ=12)
                                                                                                                                 gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat          <  1:1010545/140‑1 (MQ=16)
                                                                                                                                  cTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTCCAACAGAACCGtatt         >  2:578053/1‑141 (MQ=255)
                                                                                                                                  cTCAATGCATATCTTAGATAGTAGCCAATGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTATATTTTATGATAAATTCACCTTTGCAAAAGAACCGtatt         >  2:688106/1‑141 (MQ=11)
                                                                                                                                  cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt         >  2:858611/1‑141 (MQ=16)
                                                                                                                                  cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt         <  2:188832/141‑1 (MQ=21)
                                                                                                                                   tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt        >  2:514334/1‑140 (MQ=16)
                                                                                                                                    cAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat       >  1:786355/1‑141 (MQ=255)
                                                                                                                                     aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta      >  2:706520/1‑138 (MQ=255)
                                                                                                                                      aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat     <  2:901512/141‑5 (MQ=255)
                                                                                                                                         tataTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAt  <  1:791990/141‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                            
ATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAT  >  CP000730/36342‑36619

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: