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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC15.9% 164.8 / 24.6 89T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (37/36);  new base C (9/5);  total (46/43)
Fisher's exact test for biased strand distribution p-value = 3.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618
                                                                                                                                           |                                                                                                                                           
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                            >  2:83144/1‑141 (MQ=255)
taataaGCTTTTTGATGCTTTTTCTTCCAAATTTCATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGTACTTTTACACCACGATCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                            <  1:325140/141‑1 (MQ=14)
 aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCAATGaat                                                                                                                                           >  2:150208/1‑141 (MQ=12)
 aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtaa                                                                                                                                           >  1:455946/1‑140 (MQ=255)
  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGaata                                                                                                                                          >  2:5688/1‑141 (MQ=17)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                       <  1:727624/141‑1 (MQ=39)
     aGCTTTTTTAGTCTTTTTCTTCCATATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                       <  2:694149/141‑1 (MQ=255)
         tttttATGCTTTTTCTTCCAAATTTGATTCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                                   <  2:690943/141‑1 (MQ=25)
          ttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAt                                                                                                                                  <  1:42192/141‑1 (MQ=35)
           tttATGCTTTTTCTTCCAAATTTGATACCAAATTGTTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                                 >  1:404323/1‑141 (MQ=255)
           tttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAta                                                                                                                                 >  1:188222/1‑141 (MQ=35)
              ttttttttttttACAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGAACGTTGAATGATGAACGTATCCACTACGTTCCCTTAATTTTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                              >  2:61471/5‑141 (MQ=11)
               tGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagtag                                                                                                                             >  1:699634/1‑141 (MQ=34)
                 cTTTTTCTTCCAAATTTGATACAAAATTTGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGcc                                                                                                                           <  1:837981/141‑1 (MQ=25)
                         tcCAAATTTGATACAAACTTGGAAATGCTATTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                   >  1:830863/1‑141 (MQ=11)
                         tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                   >  1:762679/1‑141 (MQ=35)
                             aaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                               <  2:556495/141‑1 (MQ=35)
                                ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                            <  1:295031/141‑1 (MQ=35)
                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                           <  1:271504/141‑1 (MQ=35)
                                  gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATaa                                                                                                          >  2:112495/1‑141 (MQ=35)
                                     aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGGATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAt                                                                                                       >  2:92717/1‑141 (MQ=21)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTTCAGTGATAACATCCTTGTTa                                                                                               <  2:595522/141‑1 (MQ=17)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCGTAGATAGTAGCCAACGGCTACAGGGATAACATCCTTGTTa                                                                                               >  1:71894/1‑141 (MQ=17)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                              >  1:18646/1‑141 (MQ=32)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                             >  1:2874/1‑141 (MQ=32)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACATTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                            >  2:834053/1‑141 (MQ=17)
                                                    atTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTTCTTCCCGTTGTCTTAATATTTCTGATCTCTTACGATAACAAAAACTACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                        <  1:401720/141‑1 (MQ=11)
                                                      tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                      <  1:136086/141‑1 (MQ=32)
                                                      tCTGGAACCCAACGTGAGGCCCTTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAGTATTTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                      <  1:650565/141‑1 (MQ=11)
                                                       cTTGAACCCAACGGAAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtat                                                                                     >  1:437071/1‑141 (MQ=17)
                                                          gAACCCAATGGTAGACCGTTGAATGATGAACGTTTACACCACGTGCCCTTAATATGTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATc                                                                                  >  1:581478/1‑141 (MQ=255)
                                                          gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                  >  1:255631/1‑141 (MQ=32)
                                                          gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                  >  2:686882/1‑141 (MQ=32)
                                                          gAACCCAACGGGAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACATTGATAACATCCTTGTTAAATGGTTTATATc                                                                                  >  2:284525/1‑141 (MQ=12)
                                                            aCCCAAACGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                                <  1:179927/141‑1 (MQ=17)
                                                             cccAACGGTAGACCGTTGAATTCTGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTAGATCTGa                                                                               <  2:464436/141‑1 (MQ=12)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                               >  1:837695/1‑141 (MQ=32)
                                                              tcAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATTTTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAGAACATCCTTGTTAAATTGTTTATATTTGaa                                                                              <  2:508739/140‑1 (MQ=12)
                                                               cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATGGTTTATATCTGaaa                                                                             <  1:547741/141‑1 (MQ=17)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        <  1:786618/141‑1 (MQ=32)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        >  1:470960/1‑141 (MQ=32)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        >  1:470750/1‑141 (MQ=32)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        >  2:645567/1‑141 (MQ=32)
                                                                    gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                        <  2:762679/141‑1 (MQ=32)
                                                                         ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                   >  2:171911/1‑141 (MQ=32)
                                                                            gtGAATGATGAACGGTTAAACCCAGTTCTCTTAAAATTACCGGTATATCACGATCGTTCAATGTCGATCTGGGATATTAGACGTGGGCCACAGGTAGAATAGCCTTGGTATATTGGTTATATCTGAAATAGTTCATACaga                                                                <  2:198023/140‑1 (MQ=2)
                                                                             tGAATGATTAACGTTTACACCACGTGCCCTGAAGATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                               <  2:54295/141‑1 (MQ=255)
                                                                                      aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                      <  2:45716/141‑1 (MQ=32)
                                                                                      aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                      <  1:614025/141‑1 (MQ=32)
                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACGTAGTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                                                     >  1:478394/1‑141 (MQ=255)
                                                                                       aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                                                     >  2:691288/1‑141 (MQ=32)
                                                                                         gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                   >  2:425974/1‑141 (MQ=32)
                                                                                          tttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                  >  2:56874/1‑141 (MQ=32)
                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 <  2:484022/141‑1 (MQ=32)
                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                                 >  2:149682/1‑141 (MQ=32)
                                                                                             tcaccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                               <  2:349532/140‑1 (MQ=17)
                                                                                              caccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaa                                              >  1:211831/1‑141 (MQ=32)
                                                                                                  acGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                          >  2:704761/1‑141 (MQ=32)
                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATtat                                         >  2:765449/1‑141 (MQ=17)
                                                                                                        ccTTAATAGTGCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                    <  1:266046/141‑1 (MQ=17)
                                                                                                          ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                                  >  1:30189/1‑141 (MQ=33)
                                                                                                           tAATATTTTATATAAATCACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATaaa                                 <  1:412/141‑1 (MQ=11)
                                                                                                           tAATATTTCAGATATAGCACGATAACTCAATGTATATCTTAGATAATAGCCAACTGCTACAGTTATAATATCCTTGTTAAATTGTTTATATCTCAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaaa                                 <  2:247909/141‑1 (MQ=11)
                                                                                                               atgtCAGCTATATCACGACAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                             <  1:108250/138‑1 (MQ=255)
                                                                                                               atTTCCGATATATCACGATAACTCCATGTACATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                             <  1:358286/141‑1 (MQ=255)
                                                                                                               atTTCAGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                             <  1:17475/141‑1 (MQ=25)
                                                                                                                  ttaGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                          <  1:27951/139‑1 (MQ=25)
                                                                                                                    aGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACttt                        >  1:102941/1‑141 (MQ=35)
                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     >  2:333053/1‑141 (MQ=34)
                                                                                                                       tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  1:83144/141‑1 (MQ=35)
                                                                                                                       taaatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                     <  2:699634/138‑1 (MQ=21)
                                                                                                                             aCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGaa               >  1:96305/1‑141 (MQ=35)
                                                                                                                             aCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACACAAGACTCCTTTTTATTAAAATTGTACTATAAATTAACCTTGGGACCAGaa               >  1:195328/1‑141 (MQ=11)
                                                                                                                              cGATAACTCAATGTATATCTTAGATCGTACCTATCGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCACCAGAAc              >  1:404794/1‑141 (MQ=11)
                                                                                                                               gagaaCTCAATGCATATATTAGATAGTAGCCAACAGCGAAAGTGATAACATCCTTGTTAAATTGTTAATAACTGAAATAGTTCATACAGAAGACCCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc             <  1:535071/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCATCGGCTACAGTGATATCGACCTTGTTAAATTGTTTATATCTCAAATAGTTTATACATAATACTCCTTTTTTTTAAAATTATATTATTAATTCAACTTTGCAACAGAACCGt           <  2:31402/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt           >  2:252535/4‑141 (MQ=2)
                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt           <  1:143588/141‑1 (MQ=37)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATATTATAAATACAACTTTGCAACAGAACCGtat         >  2:824909/2‑141 (MQ=255)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         <  2:528203/140‑1 (MQ=16)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTACTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  1:662117/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAGCCTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat         >  1:795425/2‑141 (MQ=255)
                                                                                                                                    cTCAATGTATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatt        <  1:168676/141‑1 (MQ=38)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt        >  1:86602/1‑141 (MQ=12)
                                                                                                                                    atCAAAGCATATCGTAGATGGTCGACAACGGCTACAGTGATAACATCCTAGTTAAATTGTTTATATCTGAAATAGTTCATAAAGAAGACTTGTTTTTGTTAGAATTATATTATAAATTCAACTTTCCAACAGAACGGtatt        >  2:391030/2‑141 (MQ=11)
                                                                                                                                     tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta       >  1:332464/1‑141 (MQ=40)
                                                                                                                                     tCAATGCATATCTTAGATAGTAGCCAACCGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt       <  1:10620/141‑2 (MQ=255)
                                                                                                                                     tCAATGCATATCTTAGATAATAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt       >  1:469163/1‑140 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTAATTTAAAGTATATTATAAATTCACCTTTGCAACAGAACCGtatttt      >  1:567633/1‑139 (MQ=11)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     >  1:762678/1‑138 (MQ=255)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta     <  1:344343/141‑4 (MQ=255)
                                                                                                                                       aaTGAATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaact     <  2:345749/141‑9 (MQ=11)
                                                                                                                                        aTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATgg    <  1:766548/141‑1 (MQ=255)
                                                                                                                                          gTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGaa  <  1:92717/141‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAA  >  CP000730/36340‑36618

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: