Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,087 A→G 4.5% C155C (TGT→TGC USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,0870AG 4.5% 331.1 / 4.9 134C155C (TGT→TGCUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (63/61);  new base G (4/3);  total (70/64)
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.39e-01

Only 100 of 134 total aligned reads displayed.

TTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTA  >  CP000730/35949‑36225
                                                                                                                                          |                                                                                                                                          
tttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCGTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACa                                                                                                                                          >  1:338880/1‑141 (MQ=255)
  tatataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGt                                                                                                                                        >  1:134283/1‑141 (MQ=40)
  tatataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGt                                                                                                                                        >  2:1094188/1‑141 (MQ=255)
  tatataGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGt                                                                                                                                        <  1:817400/141‑1 (MQ=255)
   atataGAGCGTAAATACAGTCAATACCTTTTAAATTATTTTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACGTTAATATGACGGTTATCTTGCTTAATAGAGTTATTGAGTTATTCCGATGTACATTGACAGTc                                                                                                                                       <  2:258194/141‑1 (MQ=14)
      taGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAgg                                                                                                                                    <  1:834311/141‑1 (MQ=40)
         agCGTAAATACATTTCATTCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttt                                                                                                                                 <  1:716687/141‑1 (MQ=255)
         agCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttt                                                                                                                                 <  2:461103/141‑1 (MQ=40)
          gCGTAAGTACATTCAGTACGTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTGGCTTAATGAGGTTATTCAGATATTTCGATGTACAGTAACAATctgagtta                                                                                                                                >  2:691697/1‑134 (MQ=11)
          gCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttta                                                                                                                                <  1:1128035/141‑1 (MQ=255)
           cGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaa                                                                                                                               <  1:1040350/141‑1 (MQ=255)
            gTAAATACATGCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaag                                                                                                                              >  2:809995/1‑141 (MQ=255)
            gTAAATACATACAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaag                                                                                                                              >  2:810013/1‑141 (MQ=255)
             tAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGTCAGTCAGgtttaagt                                                                                                                             >  1:191892/1‑141 (MQ=255)
              aaaTACATTCAATACCTTTTAAAGTATGCATTGAAGTGTTGAAACTTTGCTACCTGGTATTTCTTACTTTAATATGACGGTGATCTTGTTTCATGAGGGTATTCAGATATTTCGATGTACAATGACCGTCAGgattaaatt                                                                                                                            >  2:520409/1‑141 (MQ=14)
               aaTACATTCAATACCTTTTGAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttt                                                                                                                           >  1:633577/1‑141 (MQ=255)
                aTACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttta                                                                                                                          >  1:413039/1‑141 (MQ=40)
                 tACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGACGTACAATGACAGTCAGgtttaagtttaa                                                                                                                         <  2:507115/141‑1 (MQ=255)
                 tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtttaa                                                                                                                         >  1:535562/1‑141 (MQ=255)
                  aCATTCAATACCTTTTAAAGTATTCTTTGTTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCGGGTTTAAGTTTaaa                                                                                                                        >  1:50142/1‑141 (MQ=255)
                     ttCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGc                                                                                                                     >  2:1190876/1‑141 (MQ=255)
                          tACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTaa                                                                                                                <  1:37881/141‑1 (MQ=255)
                           aCCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAt                                                                                                               >  1:960131/1‑141 (MQ=255)
                            ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                                                                                                              <  2:241196/141‑1 (MQ=255)
                             cTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAttt                                                                                                             >  1:1049198/1‑140 (MQ=255)
                              ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGCCAGTCAGGTTTAAGTTTAAAAGCTTTAATTAc                                                                                                            >  2:428564/1‑141 (MQ=255)
                              ttttAAAGTATTCTTTACTGTATTGATGCTTTGATACCTTGTCTTTCTTACTTTAATATTACGGTGTTCTTGCTCAATGAGGTTATTAAGATATTTTGATGTACAATGACAGTCAGGTGTAAGTTTAAAAACTTTCATTAc                                                                                                            >  1:592175/1‑141 (MQ=14)
                                 tAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACttt                                                                                                         >  2:777472/1‑141 (MQ=40)
                                      tATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCa                                                                                                    >  1:136856/1‑141 (MQ=40)
                                       aTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCAt                                                                                                   <  1:671422/141‑1 (MQ=255)
                                       aTTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCAt                                                                                                   >  1:185288/1‑141 (MQ=255)
                                        ttCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCAtt                                                                                                  <  2:414496/141‑1 (MQ=40)
                                         tCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACGATGACAGTAAATTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTg                                                                                                 <  1:573768/141‑1 (MQ=255)
                                            ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATTACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGCTa                                                                                              >  1:124734/1‑141 (MQ=255)
                                            ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTACGTTTAAAAGCTTTAATTACTTTGGCCATTGCtc                                                                                              >  1:460047/1‑140 (MQ=255)
                                            ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTa                                                                                              >  1:788923/1‑141 (MQ=255)
                                            ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTa                                                                                              >  1:789013/1‑141 (MQ=255)
                                              gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTAcc                                                                                            <  2:963832/141‑1 (MQ=40)
                                               cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGATAGTCAGGTTTAAGTTTAAAAGATTTAATTAAGTTAGctttggttacat                                                                                           >  2:1065998/1‑130 (MQ=14)
                                               cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCt                                                                                           <  1:822931/141‑1 (MQ=255)
                                                   aTTGATACTTTGATACCTTGTCTTTCTTACTTTAATGTGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGt                                                                                       >  1:828903/1‑141 (MQ=255)
                                                    ttGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAGTGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGtt                                                                                      >  1:798093/1‑141 (MQ=255)
                                                       aTACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGGTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGaa                                                                                   >  1:606248/1‑141 (MQ=255)
                                                        tACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAg                                                                                  <  1:700067/141‑1 (MQ=40)
                                                        tACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAg                                                                                  <  1:1033367/141‑1 (MQ=40)
                                                        tACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAg                                                                                  >  1:485711/1‑141 (MQ=255)
                                                         aCTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAgg                                                                                 >  1:1081243/1‑141 (MQ=40)
                                                         aCTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTAGTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCCTTGAAgg                                                                                 >  1:720868/1‑141 (MQ=255)
                                                             tGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGcc                                                                             <  2:1193472/141‑1 (MQ=39)
                                                               aTACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTg                                                                           >  2:1033851/1‑141 (MQ=39)
                                                                tACCTTGTCTTTCTTACTTTAATATGCCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                          >  2:1099455/1‑141 (MQ=38)
                                                                tACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                          >  1:1041267/1‑141 (MQ=38)
                                                                tACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGa                                                                          >  1:1167750/1‑141 (MQ=39)
                                                                    ttGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTg                                                                      >  1:1056338/1‑141 (MQ=38)
                                                                     ttgctttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATCGTGCGATGTACGATGGCAGTCATGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTAACTCCGTTGAACGTTCCTGATCTGt                                                                     <  2:375529/138‑1 (MQ=11)
                                                                      gtatttcttACTTTAATATGACGGTGCGCGTGTTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTa                                                                    <  2:484946/138‑1 (MQ=11)
                                                                      gTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGTTACGTTCGTTGAAGGTGCCTGATCTGTa                                                                    <  1:760936/141‑1 (MQ=255)
                                                                        ctttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAAt                                                                  >  1:736459/1‑141 (MQ=38)
                                                                          ttcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTa                                                                <  2:1183887/141‑1 (MQ=38)
                                                                           tcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAc                                                               <  1:535634/141‑1 (MQ=38)
                                                                           tcttACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAc                                                               <  1:58402/141‑1 (MQ=38)
                                                                            cttACTTTACTATGACGGTGATCTTTCTTAATGATGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAACAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTAcc                                                              <  1:997782/141‑1 (MQ=14)
                                                                              tACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCtt                                                            <  1:626325/141‑1 (MQ=38)
                                                                                 tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                         >  2:263990/1‑141 (MQ=35)
                                                                                 tttAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTg                                                         <  2:338880/141‑1 (MQ=34)
                                                                                     ctccgacatggATCTTGCTCAATGAGGTTATTCAGATAGTTGGATGTACAATCACAGTCAGGTTTAAGTTTAAAAGCTTTAACTACGTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGt                                                     <  1:386000/131‑1 (MQ=1)
                                                                                      tagtaCGGTGATCTTGCTCAATGAGGTTATGCAGATACTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGtt                                                    <  2:837244/137‑1 (MQ=255)
                                                                                          aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCTTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTAcc                                                >  1:908219/1‑141 (MQ=255)
                                                                                          aCGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGGTTAcc                                                >  2:593310/1‑141 (MQ=255)
                                                                                            ggTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaa                                              >  1:416712/1‑141 (MQ=33)
                                                                                             gTGATCTTGCTCAATGAGGTGATTCAGATATTTCGATTTACAATGACAGTCAGGTTTAATTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCaaa                                             <  2:527358/141‑1 (MQ=12)
                                                                                                 tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                         <  1:416175/141‑1 (MQ=32)
                                                                                                 tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGt                                         <  1:306699/141‑1 (MQ=32)
                                                                                                 tCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAACGTGCCTGATCTGTAATTACCTTTTGAGGTTTGCCAAATTGt                                         >  1:1117689/1‑141 (MQ=17)
                                                                                                     ggtCAATGTGGTTATTCAGATAGTTCGATTTACAATGACAGTCAGGTTTACGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                     <  2:784367/139‑1 (MQ=255)
                                                                                                     gCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTTAGGTTTATGTTTAAAAGCTTTAATTACTTTAGCCACTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTaa                                     <  2:784482/141‑1 (MQ=12)
                                                                                                       tCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATg                                   >  1:1115666/1‑141 (MQ=32)
                                                                                                        cAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATga                                  <  2:1049198/141‑1 (MQ=32)
                                                                                                           tGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAc                               >  1:675014/1‑141 (MQ=32)
                                                                                                            ggggTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACg                              <  1:378653/139‑1 (MQ=17)
                                                                                                            gAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGATTTACCAAATTGTTTAATGAGACg                              >  1:766131/1‑141 (MQ=12)
                                                                                                            gAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACg                              <  1:1027867/141‑1 (MQ=32)
                                                                                                             aGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                             <  1:357208/141‑1 (MQ=32)
                                                                                                             aGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGt                             <  2:633577/141‑1 (MQ=32)
                                                                                                              ggTTATTCATATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGtt                            <  2:2454/141‑1 (MQ=17)
                                                                                                                      aGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGATGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACAAAATTGTTTAATGAGAAGTTTGATaatc                    >  2:969207/1‑139 (MQ=12)
                                                                                                                        atatTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCACCTTATTTAATGAGGCGTTTGATAAAccc                  >  1:24400/1‑139 (MQ=255)
                                                                                                                         tatTTCGATGTACAATGACAGTCAGTTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCa                 <  2:183367/141‑1 (MQ=17)
                                                                                                                              cGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGc            <  2:422438/141‑1 (MQ=32)
                                                                                                                               gATGTACAATGCCAGCCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt           >  2:853842/1‑141 (MQ=17)
                                                                                                                               gATGTACAATGACAGTCAGGTTTAAGTTTAACAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATgat           >  2:853832/1‑139 (MQ=17)
                                                                                                                               gATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCt           <  1:790985/141‑1 (MQ=32)
                                                                                                                                aTGTACAATGACAGTCAGGTCTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTg          <  1:886119/141‑1 (MQ=17)
                                                                                                                                    aCAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATg      >  1:1079478/1‑141 (MQ=32)
                                                                                                                                     cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTGATGACGTTAGCCATTGTTACCTTCGGTGATGTTGCCTGCTCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCCTCTTCTGAATgt     >  2:430927/1‑140 (MQ=11)
                                                                                                                                     cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGATTGATAAACGCAGATGGTGAATGa     >  1:62088/1‑141 (MQ=12)
                                                                                                                                     cAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCAATTGCTACCTTCGTTGAAGGAGCCTGATCTATAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGAAGTTTGATAAACGCATATGCTGAATgt     >  2:524495/1‑140 (MQ=255)
                                                                                                                                      aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATTTGTAATTACCTTTTGAGGTTTTCCCAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt    <  2:936601/141‑1 (MQ=12)
                                                                                                                                      aaTGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGAt    <  1:683844/141‑1 (MQ=32)
                                                                                                                                        tGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTCGTCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTa  <  2:191892/141‑1 (MQ=14)
                                                                                                                                          |                                                                                                                                          
TTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCTTCGTTGAAGGTGCCTGATCTGTAATTACCTTTTGAGGTTTACCAAATTGTTTAATGAGACGTTTGATAAACGCATATGCTGAATGATTA  >  CP000730/35949‑36225

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: