Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,479 T→C 12.6% 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.6% 192.2 / 16.6 93T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (38/41);  new base C (5/7);  total (43/50)
Fisher's exact test for biased strand distribution p-value = 7.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36341‑36616
                                                                                                                                          |                                                                                                                                         
aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGaat                                                                                                                                         >  2:377731/1‑141 (MQ=11)
aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGtat                                                                                                                                         >  2:626945/1‑141 (MQ=255)
   aaGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATc                                                                                                                                      >  1:228653/1‑141 (MQ=38)
            tATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagt                                                                                                                             >  2:35258/1‑141 (MQ=34)
             aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                            >  1:492918/1‑141 (MQ=35)
             aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtagta                                                                                                                            >  1:80131/1‑141 (MQ=35)
               gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                          >  2:134848/1‑141 (MQ=35)
               gCTTTTTCTTCCAAATTTGATACAAAATCGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                          >  1:495119/1‑141 (MQ=25)
                cTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGcg                                                                                                                         >  1:647411/1‑140 (MQ=21)
                    ttcttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACg                                                                                                                     <  1:134848/141‑1 (MQ=35)
                      cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                   <  1:384518/141‑1 (MQ=35)
                         ccAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTATCCAACGGCTAc                                                                                                                <  2:168121/141‑1 (MQ=25)
                            aaTTTGATAAAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGt                                                                                                             <  1:185458/141‑1 (MQ=17)
                                tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                         <  1:506892/141‑1 (MQ=35)
                                 gATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACATCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATaa                                                                                                        >  2:727337/1‑141 (MQ=25)
                                  aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATAACGATAACTCAATGTATAACTTAGATAGTAGCCAACGGCTCCAGTGTTAAc                                                                                                       >  1:175958/1‑141 (MQ=255)
                                    aCAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACac                                                                                                     >  1:363742/1‑140 (MQ=18)
                                         aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                                >  1:376256/1‑141 (MQ=32)
                                         aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCAAACGGCTACAGTGATAACATCCttt                                                                                                >  2:71592/1‑140 (MQ=17)
                                          ttGGGGCATATTCTTGACCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                               <  1:451276/141‑1 (MQ=17)
                                          ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                               <  1:686244/141‑1 (MQ=32)
                                          ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                               <  2:187890/141‑1 (MQ=32)
                                          ttGGGGCATATTCTTCAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                               <  2:565699/141‑1 (MQ=17)
                                            ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTa                                                                                             >  2:453354/1‑141 (MQ=32)
                                                cATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATACTGGTTAAAtt                                                                                         >  1:543031/1‑141 (MQ=17)
                                                  tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                       <  1:431976/141‑1 (MQ=32)
                                                  tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                       <  1:431899/141‑1 (MQ=32)
                                                  tatTCTGGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                       <  1:171503/141‑1 (MQ=17)
                                                      cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtat                                                                                   >  1:360401/1‑141 (MQ=32)
                                                      cTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtat                                                                                   >  1:135965/1‑141 (MQ=32)
                                                       ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                  <  1:302797/141‑1 (MQ=32)
                                                       ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                  >  1:248243/1‑141 (MQ=32)
                                                       ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                  <  2:341040/141‑1 (MQ=32)
                                                         gaccccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                <  2:596757/138‑1 (MQ=17)
                                                           aCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTg                                                                              <  1:293254/141‑1 (MQ=32)
                                                              cAACGGTAGACCGTTGAATTATGAACGTTTACACCGCGTTCCCTCATTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                           <  2:747462/141‑1 (MQ=255)
                                                              cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                           <  1:297646/141‑1 (MQ=32)
                                                              cAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGaaa                                                                           >  1:267058/1‑141 (MQ=32)
                                                                  ggtctaCCGTTGGATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                       >  2:46287/6‑141 (MQ=12)
                                                                    tctaCCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtta                                                                     >  2:425417/4‑140 (MQ=17)
                                                                      gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                   >  2:25809/1‑141 (MQ=32)
                                                                         cGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACa                                                                <  2:360401/141‑1 (MQ=32)
                                                                           ttGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACaga                                                              <  1:297053/141‑1 (MQ=32)
                                                                             gtgtgatgcACCTTTAGACCACGTGTCTTTGATTTTTCAGGGATATGATGATAACTTAGTGTATGTCTTCGATAGGGGCCCACGGCTACAGGGATAACATCCTGGTTAAATTGTTTATATCTGAAATAGTTCATACagaag                                                            <  2:597678/138‑1 (MQ=2)
                                                                                tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                         <  2:91867/141‑1 (MQ=32)
                                                                                tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                         <  1:13135/141‑1 (MQ=32)
                                                                                  atgaACGTTTACACCACGTTCAATGGTGTCAGCATATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCGACGGCTACAGTGATAACATCCTTGTTAACTTGTTTATATGTGAAATAGTTCTTCCAGAAGACTcc                                                       >  1:184189/1‑141 (MQ=11)
                                                                                    gaACGTGTACCCCACGTTGCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtt                                                     <  2:336834/141‑1 (MQ=12)
                                                                                      aCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCtttt                                                   >  2:249089/1‑141 (MQ=32)
                                                                                       cGTTTACATCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                  <  1:727337/141‑1 (MQ=17)
                                                                                       cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                  <  2:134192/141‑1 (MQ=32)
                                                                                       cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                  <  2:134093/141‑1 (MQ=32)
                                                                                        gTTTACACCACGTTCCCTTAATATTTCAGTTATATCACGACAACTCAAGGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGATTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                 <  2:248243/141‑1 (MQ=255)
                                                                                        gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                 >  2:510075/1‑141 (MQ=32)
                                                                                         tgtACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                <  2:385198/139‑1 (MQ=17)
                                                                                          ttACACCACGTTCCCTTAATATTTCATATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAAACAGAAGACTCCTTTTTGtt                                               <  1:465239/141‑1 (MQ=17)
                                                                                          tagcaccacGTTCACTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCTAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                               <  2:570136/138‑1 (MQ=12)
                                                                                           tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                              <  2:495119/141‑1 (MQ=32)
                                                                                           gtgtaaacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                              >  2:539844/7‑141 (MQ=17)
                                                                                            tgtaaacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                             >  2:1732/6‑141 (MQ=17)
                                                                                             caccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTCGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaa                                            >  1:483533/1‑141 (MQ=17)
                                                                                              accacGTTCACTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAGGTAATACCACTGAAAGAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                           <  1:369614/141‑1 (MQ=11)
                                                                                                cgcGTTCCCTTTATATTTCAGATATATCTCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAtt                                         <  2:432688/139‑1 (MQ=12)
                                                                                                 acGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATta                                        >  1:652610/1‑141 (MQ=32)
                                                                                                    ttCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGGTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                                     >  1:724255/1‑141 (MQ=17)
                                                                                                           aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                              <  1:650688/141‑1 (MQ=35)
                                                                                                            atgtgtCAGATATATCACGATAACTACATGAATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                             <  2:485960/136‑1 (MQ=11)
                                                                                                                  cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                       >  1:93647/1‑141 (MQ=35)
                                                                                                                   aGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACttt                      <  1:578819/141‑1 (MQ=35)
                                                                                                                    aatatatCACGATAACTCAATGCATATCTTAGATGGTAGCCAGCAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                     <  1:471335/140‑1 (MQ=255)
                                                                                                                        tatCCCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                 <  1:207680/141‑1 (MQ=25)
                                                                                                                        tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaaa                 >  2:402712/1‑140 (MQ=21)
                                                                                                                        tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                 >  1:484238/1‑141 (MQ=35)
                                                                                                                        tatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                 >  1:745558/1‑141 (MQ=35)
                                                                                                                        taaaagcgtctcTCAATGGATACCATCGTTAATAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAc                 <  1:361769/130‑1 (MQ=11)
                                                                                                                           cACGATAACTCAATTTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa              <  1:336527/141‑1 (MQ=25)
                                                                                                                           cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGATGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa              <  2:492918/141‑1 (MQ=25)
                                                                                                                            aCGCTAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGaa             >  1:552389/1‑141 (MQ=25)
                                                                                                                               aTAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg          <  2:445236/141‑1 (MQ=9)
                                                                                                                                tagctgaaTGCATTTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:113210/135‑1 (MQ=11)
                                                                                                                                tagcTCAATGCATATTTTAGATAGTAGCCAACGGCTAAAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCCTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  2:197692/138‑1 (MQ=11)
                                                                                                                                tAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGt         >  2:8100/1‑141 (MQ=37)
                                                                                                                                caaCGCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTACAACAGAACCGt         >  2:702215/2‑141 (MQ=255)
                                                                                                                                 agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCCGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAAATTTGCAACAGAACCGta        >  2:696215/3‑141 (MQ=11)
                                                                                                                                 agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCATGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta        >  2:160272/3‑141 (MQ=11)
                                                                                                                                 aaCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGta        >  1:221879/1‑141 (MQ=37)
                                                                                                                                  gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTATTTATACCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTCTGCAACAGAACCGtat       <  2:654891/140‑1 (MQ=11)
                                                                                                                                  aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat       >  1:635992/1‑141 (MQ=39)
                                                                                                                                  aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat       <  1:17649/141‑1 (MQ=38)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      <  1:686373/141‑1 (MQ=21)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTGTTGGTAAAAATATATTATAAATTAAACTTTGCCACAGAACCGtatt      >  1:630094/1‑141 (MQ=11)
                                                                                                                                    tCAATGTATATCTTATATAGTAGCCAACGGCTACAGTGATAACATCCTTTTTAAATTTTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta     <  2:184189/141‑1 (MQ=255)
                                                                                                                                     cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    >  1:491910/1‑139 (MQ=255)
                                                                                                                                      aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta   >  1:756034/1‑138 (MQ=255)
                                                                                                                                       aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat  <  2:549090/141‑5 (MQ=255)
                                                                                                                                          |                                                                                                                                         
AATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36341‑36616

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