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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC12.0% 257.9 / 25.6 126T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (58/51);  new base C (8/7);  total (67/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 130 total aligned reads displayed.

TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATG  >  CP000730/36340‑36615
                                                                                                                                            |                                                                                                                                        
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                         <  2:226426/141‑1 (MQ=255)
taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGta                                                                                                                                         >  1:429159/1‑141 (MQ=255)
  ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGgata                                                                                                                                       >  2:630090/1‑141 (MQ=38)
     aGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTGACACAACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCt                                                                                                                                    >  2:519723/1‑141 (MQ=255)
        ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAg                                                                                                                                 >  1:649056/1‑141 (MQ=37)
            ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATACCTTAGAtag                                                                                                                             >  1:617848/1‑141 (MQ=25)
              tttttttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGCTAACTCAATGTATATCTTAGAtagta                                                                                                                           >  2:563336/5‑141 (MQ=21)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGCATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                           >  1:868803/1‑141 (MQ=25)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                           >  1:84650/1‑141 (MQ=35)
              aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                           >  1:184727/1‑141 (MQ=35)
                gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATATTAGc                                                                                                                         >  2:921500/1‑141 (MQ=21)
                    tttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGccgtt                                                                                                                     <  1:544293/141‑4 (MQ=25)
                       cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                  >  2:489987/1‑141 (MQ=35)
                       cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                  >  2:379338/1‑141 (MQ=35)
                            aaaTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAg                                                                                                             >  2:774259/1‑141 (MQ=35)
                                 tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                        >  1:795137/1‑141 (MQ=35)
                                   aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTCTCACGATAACTCAATGTATATCTTAGATAGTAGCCAGCGGCTACAGTGATaaa                                                                                                      >  1:26528/1‑140 (MQ=255)
                                   aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAAAATGATGCAAAATGAACTTTTAGACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                      <  1:443933/141‑1 (MQ=11)
                                          aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                               <  1:289927/141‑1 (MQ=32)
                                            tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                             <  2:288636/141‑1 (MQ=17)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTCCAGTGATAACATCCTTGTTa                                                                                            >  2:48257/1‑141 (MQ=12)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATGACATCCTTGTTa                                                                                            >  2:134463/1‑141 (MQ=17)
                                             ggggCATATTCTTGAACCCAACGGTAGACCGCTGAATGATGAAAGTTTACGCAACGCAACTTTAATACTTTTGTTATATGCTGATAAATTGATCTAAAGCTTAAGGGGTAGCAAACTGCTACAGTTTTAACATCCTTGTTa                                                                                            >  2:785553/1‑141 (MQ=2)
                                              gggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                           <  2:84650/141‑1 (MQ=32)
                                               ggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaaa                                                                                          >  2:256505/1‑141 (MQ=32)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGATTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTACATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                         <  1:630090/141‑1 (MQ=12)
                                                gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                         >  2:653764/1‑141 (MQ=32)
                                                  atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTg                                                                                       <  1:286313/141‑1 (MQ=17)
                                                   tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                      <  2:27748/141‑1 (MQ=32)
                                                    atTCTTGAACCCAACGGTAGGCCGTTGAATGATGAACGTGTACACCACGTTCCCTTAACATTTCAGATATATCACGATCACTCAATGTATATCTTAGAGAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                     <  1:37470/141‑1 (MQ=11)
                                                      tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCGTTAATATTTCAGATATATCACGATAACTAAATTTATTTTTTAGTTAGTGGCTAACAGCTACAGTGGTAACATACTAGTTAAATTGTTta                                                                                   >  2:287431/1‑141 (MQ=11)
                                                        ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                 >  2:571846/1‑141 (MQ=32)
                                                         tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                                >  1:153761/1‑141 (MQ=17)
                                                          tgactcAATGGTAGAACCTAGAAATATGACCGTTCACACCACGTTCTCTTAACATTTCAGATATATCACGATAGATAAGTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                               <  1:132173/136‑1 (MQ=11)
                                                          gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTGTATc                                                                               >  2:894903/1‑141 (MQ=17)
                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                              <  2:879713/141‑1 (MQ=32)
                                                           aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                              <  1:512041/141‑1 (MQ=32)
                                                            aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                             >  1:819449/1‑141 (MQ=32)
                                                             cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                            >  2:538986/1‑141 (MQ=32)
                                                              ccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaa                                                                           <  1:691818/141‑1 (MQ=32)
                                                               cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                          <  2:367838/141‑1 (MQ=32)
                                                                 aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                        >  2:152787/1‑141 (MQ=32)
                                                                   ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                      <  1:530140/141‑1 (MQ=32)
                                                                     tAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                                                                    >  1:806273/1‑141 (MQ=32)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAACTCGTTCAt                                                                  >  2:752821/1‑141 (MQ=17)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                  >  2:817193/1‑141 (MQ=32)
                                                                       gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                  >  1:755353/1‑141 (MQ=32)
                                                                         tcGTTGAAGGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                <  2:701723/140‑1 (MQ=12)
                                                                         ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGATACAGTCATACCATCCTTGTTAAATTTGTTATATTTTAAATAGTTCATAc                                                                >  2:267128/1‑141 (MQ=11)
                                                                             tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                            >  2:100273/1‑141 (MQ=32)
                                                                             tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                            <  1:758136/141‑1 (MQ=32)
                                                                             tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                            >  1:224502/1‑141 (MQ=32)
                                                                                atgattAACGTGTACACCAGGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAc                                                         <  2:410920/141‑1 (MQ=12)
                                                                                  gatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTc                                                       <  1:418700/141‑1 (MQ=32)
                                                                                    tgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATGTTAGATAGTAGCCAACGGCTACAGTGATAACATCCGTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGAATCCt                                                     >  2:113421/1‑141 (MQ=12)
                                                                                      taCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                   <  2:411153/140‑1 (MQ=17)
                                                                                      aaCTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                   <  2:80749/141‑1 (MQ=17)
                                                                                        cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATGTTGTTCAGACAGAAGACTCCttttt                                                 <  2:669330/141‑1 (MQ=255)
                                                                                        cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                 <  2:73262/141‑1 (MQ=32)
                                                                                         gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                >  1:18119/1‑141 (MQ=32)
                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCGTACAGAAGGCTCCTTTTTGtt                                              <  1:489987/141‑1 (MQ=17)
                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                              <  1:379338/141‑1 (MQ=32)
                                                                                            tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCCTACAGAAGCCTCCTTTTTGTTa                                             >  2:303268/1‑141 (MQ=17)
                                                                                            tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAACCAAGTTAAATTGTTTATATCTGAACGAGTTCGTACAGAAAACTCCTTTTTGTTa                                             >  2:150523/1‑141 (MQ=11)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                            >  1:858728/1‑141 (MQ=32)
                                                                                                   cgttatcTTAATATTACGGATATATCACGAAAGATCAGTTGTAAACTTAGATTGTGGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATtat                                      >  1:452373/7‑141 (MQ=11)
                                                                                                   cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTTTTTATATTTGACATAGTTCATAAAGACGACTACTTTTTGTTTAAATtat                                      >  1:523321/1‑141 (MQ=11)
                                                                                                       cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtt                                  >  1:810710/1‑140 (MQ=17)
                                                                                                        ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTCTATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                 >  2:617021/1‑141 (MQ=17)
                                                                                                        ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                 >  1:137314/1‑141 (MQ=32)
                                                                                                        ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                 >  1:377630/1‑141 (MQ=32)
                                                                                                          ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAGGACTCCTTTTTGTTAACATTATACTCTaa                               >  1:74860/1‑141 (MQ=12)
                                                                                                             atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                            <  2:367354/141‑1 (MQ=35)
                                                                                                             atatTTCAGATATAACACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                            <  2:120770/141‑1 (MQ=25)
                                                                                                              aatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATACTATAAATTc                           <  2:641276/140‑1 (MQ=25)
                                                                                                                  actgATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGTAGACTCCTGTTTGTTAAAATTATACTATAAATTCAACt                       <  1:187929/138‑1 (MQ=255)
                                                                                                                     gATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAATTTCAACTTTg                    <  1:683885/141‑1 (MQ=255)
                                                                                                                      atatatCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                   <  1:587874/141‑1 (MQ=255)
                                                                                                                      atatatCACGATAGCTCAATGGATCTCATGGAGAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                   >  1:456209/1‑141 (MQ=11)
                                                                                                                        atatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATATCATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                 <  1:538986/141‑1 (MQ=25)
                                                                                                                         tatCACGATAGCTCAATGTATATATTAGTTAGTAGCCAACGGCTACGGTGATAACATCCTTGTTAAATTGTTGATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCGACTTTGCAAc                <  2:184198/141‑1 (MQ=11)
                                                                                                                         tatCACGATAACTCAATGTAAATTGTAGCTAGTATCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATACCTGAAATATTTCATACAGAAGACTCCTTTTTGTTAAAATTCTACTATACATTCAACTTTGCaca                >  2:435309/1‑139 (MQ=11)
                                                                                                                            cACGATAGCTCAATGAATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGa             <  2:66054/141‑1 (MQ=12)
                                                                                                                            cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATACTATAAATTCAACTTTCCacccga             >  1:116292/1‑136 (MQ=255)
                                                                                                                                 ttgcttaaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:582/135‑1 (MQ=255)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:727138/138‑1 (MQ=12)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  1:868682/138‑1 (MQ=9)
                                                                                                                                 aagcctaaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt        <  2:324286/135‑1 (MQ=255)
                                                                                                                                  agcTCAATTCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta       <  2:598123/139‑1 (MQ=255)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta       >  1:244209/3‑141 (MQ=9)
                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat      >  2:828061/1‑141 (MQ=38)
                                                                                                                                   aCTCAATGTATATCTAAAATAGTAGCCAACGGCTACAGTGATCACATCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat      <  1:880897/141‑1 (MQ=255)
                                                                                                                                    cTCAATGCATATGTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt     >  1:665280/1‑141 (MQ=255)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTAATTTTATTATAAATTCAACTTTGCAACAGAACCGtatt     >  2:882925/1‑141 (MQ=11)
                                                                                                                                     tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta    <  2:502065/141‑1 (MQ=255)
                                                                                                                                     tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTATTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta    <  2:502011/141‑1 (MQ=255)
                                                                                                                                     tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt   >  1:77696/1‑140 (MQ=16)
                                                                                                                                      cAATGTATAGCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat   <  1:823997/141‑1 (MQ=255)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttt  >  2:469023/1‑138 (MQ=255)
                                                                                                                                       aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATACAACTTTGCAACAGAACCGtatttta  >  2:583813/1‑138 (MQ=255)
                                                                                                                                            |                                                                                                                                        
TAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATG  >  CP000730/36340‑36615

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