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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC32.4% 131.7 / 83.8 113T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base C (18/18);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 115 total aligned reads displayed.

ATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36342‑36616
                                                                                                                                           |                                                                                                                                         
ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCATTGtata                                                                                                                                        >  2:54104/1‑141 (MQ=255)
ataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtcta                                                                                                                                        >  2:378104/1‑141 (MQ=255)
       tttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACAAAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGa                                                                                                                                 >  1:625013/1‑141 (MQ=25)
        ttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCAGGTTCCCTTAATATTTAAGATATATCACGATAACTCAATGAATATCTTAGaa                                                                                                                                >  2:130693/1‑140 (MQ=12)
          ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                                                                                              <  1:100280/141‑1 (MQ=35)
          ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAACTTAGAtag                                                                                                                              >  1:597979/1‑141 (MQ=21)
                    tcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACgg                                                                                                                    >  2:547771/1‑141 (MQ=35)
                     cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                   <  2:663807/141‑1 (MQ=35)
                     cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTACATCTTAGATAGTAGCCAACGGc                                                                                                                   >  1:460684/1‑141 (MQ=25)
                        ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAAATAGTAGCCAACggatgc                                                                                                                >  2:478196/1‑137 (MQ=25)
                        ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                                                                                >  1:197381/1‑141 (MQ=25)
                        ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAAATATATCACGATAACTCAATGTATATCTTAGATGGTAGCCAACGGCTAc                                                                                                                >  1:286692/1‑141 (MQ=255)
                            atccgATACAAAATTGGTGCTTATTCTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATGTCATATATATCACGATAACTCAATGTATATCTTAGTTAGTCGCCAACGGCTACAGTg                                                                                                            <  1:342458/137‑1 (MQ=11)
                              ttGATACAAAATTNGGGCATATTCTTGAACNCAGCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTCATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  1:244094/141‑1 (MQ=19)
                              ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTATAGTGAt                                                                                                          <  2:368021/141‑1 (MQ=25)
                              ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  1:265276/141‑1 (MQ=35)
                              ttGATACAAAATTGGGGAATTGTTCTTAACCCAACGGTAGCTAGACCGAACGGGAACGTTCACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  1:651257/141‑1 (MQ=11)
                               tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                         >  1:146886/1‑141 (MQ=35)
                               tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCATAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                         >  1:465783/1‑141 (MQ=25)
                                 aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                       <  2:502788/141‑1 (MQ=35)
                                  tACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACa                                                                                                      >  2:563279/1‑141 (MQ=35)
                                  tACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATAATTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCCACGGCTACAGTGCTAACa                                                                                                      >  2:482998/1‑141 (MQ=255)
                                       aaTTGGTGCATATTCTTGAACCCAACGGTAGACCGTTGAATGAAGACCGTTTTCACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                                                                                 <  2:270466/141‑1 (MQ=255)
                                          tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTGACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                              <  1:123703/141‑1 (MQ=17)
                                                atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTg                                                                                        >  1:263384/1‑141 (MQ=32)
                                                  tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGCTAACATCCTTGTTCAATTGt                                                                                       >  1:599675/1‑141 (MQ=17)
                                                     tCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTta                                                                                    <  1:547771/141‑1 (MQ=32)
                                                          aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTACTATTTCAGATATATCACGATAACTCAATGTATCTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTCAATTGTTTATATCt                                                                               >  2:24713/1‑141 (MQ=12)
                                                          aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                               <  1:519093/141‑1 (MQ=32)
                                                            cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                             <  1:37194/141‑1 (MQ=32)
                                                                  ggTAGAGCGTTGAATGATGAACGTTTACAGCCCGTGCCCTAAATATTTCAGATATATCACGATAACACAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATATCTGAAATAGt                                                                       >  1:489431/1‑141 (MQ=11)
                                                                   gTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                      <  1:481845/141‑1 (MQ=32)
                                                                      gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATAAAGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAAATCCTTGTTAAATTGTTTATATCTGAACTAGttact                                                                   >  1:136728/1‑138 (MQ=11)
                                                                        ccGTTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                 <  1:167627/141‑1 (MQ=17)
                                                                          gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                               >  2:591781/1‑141 (MQ=32)
                                                                          gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACag                                                               >  1:667046/1‑141 (MQ=32)
                                                                          gTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTCTTTATATCTCAAATAGTTCATACag                                                               <  1:334069/141‑1 (MQ=17)
                                                                               aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                           <  1:355511/141‑1 (MQ=32)
                                                                                     gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                                                     >  1:29860/1‑141 (MQ=32)
                                                                                     gaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCTTGTTAAATTGTTTAGATTTGAATTATTTAACAAAGATAACTCCtt                                                     >  2:520407/1‑141 (MQ=11)
                                                                                      aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                    <  2:258131/141‑1 (MQ=32)
                                                                                      aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                    <  2:275691/141‑1 (MQ=32)
                                                                                           ttACACCACGTTCCCTTAATATTTCCGCTATATTACGATCACTCAATGTAGATCTTCGATAGTAGCCAACGGCTACAGGGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACGGAAGACTCCTTTTTGtt                                               <  1:137284/141‑1 (MQ=11)
                                                                                           ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                               <  1:368126/141‑1 (MQ=32)
                                                                                             acacCACGTTCCCTTAATATTTCAGATATATGCCGATAACTCAGTGTATGTCTTAGATAGTAGCCAACGGCTACAGTGATGACAGCCTTGTGAAATTGTTTATATCTGCAATAGTTCATACAGAATCAGCATTTTTGTTaa                                             >  2:282174/1‑141 (MQ=11)
                                                                                               accacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                           >  2:530222/1‑141 (MQ=32)
                                                                                                          ttAATATTTCAGATATATCACGCTAACTCAATGTATAACCTAGATAGTAGCCAACGGCTACAGTGTTAACATCCTTGTCAAATTGTTTATATCTAAAATAGTTCATACAGAAGTGTCGTTTTTGTTAAAATAGTACTATaa                                >  2:48732/1‑141 (MQ=11)
                                                                                                             atctttCAGCCATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACTGCTACAGTGATAACATCCTTGTTAAATTGTTCATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAAtt                             <  1:554685/138‑1 (MQ=1)
                                                                                                              tatTTCAGATATATCACGATAGCTCAACGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTc                            <  1:385854/141‑1 (MQ=255)
                                                                                                               atTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                           <  1:201487/141‑1 (MQ=35)
                                                                                                                tttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCaa                          <  2:496702/141‑1 (MQ=35)
                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                        >  1:174276/1‑141 (MQ=35)
                                                                                                                  tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                        >  1:461650/1‑141 (MQ=35)
                                                                                                                    ggATCTATCACGATAACACAAGTCATAACGTGGAGAGTAGCCAATAGCTACAGTGATAACAGCCTTGTTAAATTGTTTACATCTGAAATAGGTCATACAGAAGACTCCTTCTTGTTAAAATTATACTATAAATTCAACttt                      <  2:19385/140‑1 (MQ=255)
                                                                                                                     gATATATCACGATAGCTAAGTGTAGATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTg                     >  1:490584/1‑141 (MQ=11)
                                                                                                                        atacgcggctagcTCAATGAATAGCTTAGATAGTAGCCAACGGCTATAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCaa                  <  2:520474/129‑1 (MQ=11)
                                                                                                                         tagcACGATAACTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAc                 <  2:183814/138‑1 (MQ=2)
                                                                                                                           tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAg               >  1:288255/1‑141 (MQ=35)
                                                                                                                             aCGATAACTCAATGTATATCTTAGATAGTAGACAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTTTTTTTTTCTATAAATTCGACTTTTGATCAGaa             >  2:140991/1‑141 (MQ=11)
                                                                                                                               gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc           >  2:550425/1‑141 (MQ=34)
                                                                                                                                aTAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCg          >  1:573116/1‑141 (MQ=35)
                                                                                                                                 tagcccAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:227555/136‑1 (MQ=255)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAATGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:57855/138‑1 (MQ=255)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTGAACTTTGCAACAGAACCGt         >  1:424805/4‑141 (MQ=255)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         >  1:350448/4‑141 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:135747/138‑1 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:262110/138‑1 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         >  1:163003/4‑141 (MQ=9)
                                                                                                                                 tagcTCAATGCATATCTTAGATAGTAGCCAACGGATACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGATATAACACTATACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:126915/138‑1 (MQ=11)
                                                                                                                                 tagcTCAATGCATATCTAAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATCGTTTATATCTGAAATAGTTCCTACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:447164/138‑1 (MQ=11)
                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt         <  1:638325/141‑1 (MQ=35)
                                                                                                                                 tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTCTGCAACAGAACCGt         <  1:356719/141‑1 (MQ=37)
                                                                                                                                 aagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:421961/138‑1 (MQ=12)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta        >  1:581353/3‑141 (MQ=12)
                                                                                                                                  agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta        >  2:243286/3‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  1:35693/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  1:607485/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  2:350177/2‑141 (MQ=12)
                                                                                                                                   gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       <  2:133010/140‑1 (MQ=16)
                                                                                                                                   gcTCAATGCATATCTCAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       >  1:333509/2‑141 (MQ=12)
                                                                                                                                   aCTCAATGTATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat       <  2:42401/141‑1 (MQ=38)
                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTCAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAACGtat       >  1:486273/1‑141 (MQ=255)
                                                                                                                                   aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat       >  2:349187/1‑141 (MQ=38)
                                                                                                                                    cacAATGCATATCTTAAATAGTCGCCAACGGCTACAGTGATAACACCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAAAGAACCGtatt      >  2:584778/3‑141 (MQ=11)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTATAATTATTTTATAAATTCAACTTTGCAGCAGAACCGtatt      >  1:247940/1‑141 (MQ=11)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  1:509275/1‑141 (MQ=16)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  2:620215/1‑141 (MQ=16)
                                                                                                                                    cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAACTCAACCTTGCAACAGAACCGtagt      >  1:447866/1‑139 (MQ=11)
                                                                                                                                    atCAATGTATATGTTAGATGGCAGCCAACGGCGACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatt      <  1:24713/140‑1 (MQ=255)
                                                                                                                                     tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt     >  1:77760/1‑140 (MQ=16)
                                                                                                                                      cAATGCATATCTTAGTTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:582522/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGCGAGTAGCCAACGGCTACAGTGATAACATCCGTTTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:166017/141‑3 (MQ=11)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  2:247758/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  2:171153/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:199074/141‑3 (MQ=255)
                                                                                                                                      cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    >  1:332808/1‑139 (MQ=16)
                                                                                                                                      cAATGCATATCTTAGATAGTACAAAACGTCTGCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    >  1:659988/1‑139 (MQ=11)
                                                                                                                                       aaTGTATATCTTAGATAGTAGCCAACGTCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg   <  2:507576/141‑1 (MQ=255)
                                                                                                                                       aaTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg   <  2:613547/141‑1 (MQ=255)
                                                                                                                                        aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat  <  1:136447/141‑5 (MQ=255)
                                                                                                                                           |                                                                                                                                         
ATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36342‑36616

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: