Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,413 C→T 46.6% V47I (GTC→ATC)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4130CT46.6% 139.6 / inf 573V47I (GTC→ATC) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base C (169/132);  new base T (133/132);  total (309/264)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

Only 100 of 579 total aligned reads displayed.

AATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAC  >  CP000730/36273‑36553
                                                                                                                                             |                                                                                                                                            
  tGGCACGATATAATTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCg                                                                                                                                            <  1:14631/141‑1 (MQ=18)
    gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGCTATCAAATTGGGGCCT‑CTTCTTGAACCCAAAAGTAGACCGtt                                                                                                                                          >  1:582690/1‑141 (MQ=11)
    gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGtt                                                                                                                                          >  1:238599/1‑141 (MQ=21)
    gCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGtt                                                                                                                                          >  2:631270/1‑141 (MQ=21)
      aCGATTTACATGGCTCCATTTTCCTGTTATTTTGATGTACTTCTTATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCATAGTTGATTTATAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                                                        <  2:33703/141‑1 (MQ=11)
      aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTGTTTTTTTTTCCAAATTTGATATAAAATTGGGGTAT‑ATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                                                        >  2:47462/1‑141 (MQ=11)
      aCGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGa                                                                                                                                        >  1:444989/1‑141 (MQ=25)
       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGaa                                                                                                                                       >  1:17557/1‑141 (MQ=255)
       cGATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCAGTTTATAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTAAGGCAT‑ATTCTTGAACCCAACGGTAGAACGTTGaa                                                                                                                                       >  1:647994/1‑141 (MQ=11)
        gATATAAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAat                                                                                                                                      >  1:423485/1‑141 (MQ=255)
         atatAAATTGTCCCATTTTCCTTTTACTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAatg                                                                                                                                     <  2:259350/141‑1 (MQ=11)
         atatAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAatg                                                                                                                                     <  2:32706/141‑1 (MQ=21)
          tataAATAGCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAatga                                                                                                                                    <  1:360396/141‑1 (MQ=255)
          gataAATTTCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAatga                                                                                                                                    >  2:660845/2‑141 (MQ=11)
           ataAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAatgat                                                                                                                                   <  2:106249/141‑1 (MQ=40)
            taAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCCAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGAACGttatataatg                                                                                                                                  >  2:62041/1‑133 (MQ=11)
             aaaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTTTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTAATATAAAATTGAGGCAT‑AATCTTGAACCCAACGTTAGCCCGTTGAAtgatga                                                                                                                                 >  2:221585/1‑141 (MQ=11)
             aaaTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCCTTTTCTTCCAAATCTGATACAAAATTGGGGCAC‑AGTCCTAACCCCAACGGCAGACCGTTGACtgatga                                                                                                                                 >  2:157664/1‑141 (MQ=14)
               aTAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTACCAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATCTACAATTGGTGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAc                                                                                                                               >  1:156670/1‑141 (MQ=11)
                tAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCTATTCGCCACTTAGAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACg                                                                                                                              <  1:51315/141‑1 (MQ=11)
                 aGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAAAGGTAGGCCGTGGAATGATGAACGt                                                                                                                             <  1:661431/141‑1 (MQ=255)
                  ggtCCATTTTCCTTGTATTTTGATCTACGTACCATCAATACGCCATTTGTAATAAGCTTTTTTATTCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                                                            <  2:13325/139‑1 (MQ=11)
                  gCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAc                                                                                                                               >  1:55613/1‑138 (MQ=21)
                  gCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAAc                                                                                                                               <  2:55613/138‑1 (MQ=21)
                  gCTCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGtt                                                                                                                            <  2:443732/141‑1 (MQ=255)
                    tCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTTTTTTTTTTCTTCCAAATTTGATATAAAATTGTGGCAT‑ATTCTTTAACCCAACGGAAGATCGTTGATTGTTGAACGtttt                                                                                                                          >  2:243465/1‑140 (MQ=11)
                    tCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                                                                                                                          >  2:498164/1‑141 (MQ=255)
                    tCCATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTa                                                                                                                          <  1:26721/141‑1 (MQ=255)
                     ccATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTac                                                                                                                         >  2:535492/1‑141 (MQ=21)
                      cATTTTCTTTTTATTTTGATGTACGTCTCAATAATACGTCTTTTGGAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTaca                                                                                                                        <  1:123353/141‑1 (MQ=11)
                      cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTcca                                                                                                                        >  2:379057/1‑141 (MQ=255)
                      cATTTTCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTaca                                                                                                                        >  1:162140/1‑141 (MQ=255)
                       aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                                                                                                                       <  2:563621/141‑1 (MQ=21)
                       aTTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTacac                                                                                                                       >  2:207029/1‑141 (MQ=40)
                        ttttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacc                                                                                                                      >  2:481845/1‑141 (MQ=40)
                         tttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTAcacca                                                                                                                     >  2:599878/1‑141 (MQ=25)
                         tttCCTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAcacca                                                                                                                     <  1:537691/141‑1 (MQ=255)
                          ttCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGCACCCAACGGTCGACCGTTGATCGATGaaagattcaacctc                                                                                                                    <  2:389651/141‑13 (MQ=11)
                            ccTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCCCGt                                                                                                                  >  1:218275/1‑141 (MQ=25)
                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  1:564664/1‑141 (MQ=21)
                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  1:137707/1‑141 (MQ=25)
                             cTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGGTGAATGATGAACGTTTACACCACgat                                                                                                                 >  2:295042/1‑139 (MQ=255)
                             cTTTTATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGtt                                                                                                                 >  1:529243/1‑141 (MQ=255)
                              ttttatttttatGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTc                                                                                                                <  2:93805/141‑1 (MQ=18)
                                ttATTTTGATATACGTCTCCTCAATACGCCGTTTGTAATAAGCTCTTTTATGCTTTTTATTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGGTAAATGATGAACGATTaaagccacttctc                                                                                                              >  2:54495/1‑129 (MQ=11)
                                ttATTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTctc                                                                                                              <  2:590916/141‑3 (MQ=255)
                                  atttgcaTGTACGTCTCATCAATATGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATGCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCtt                                                                                                            <  2:45936/135‑1 (MQ=11)
                                  aTTTTGATATACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCtt                                                                                                            <  1:91864/141‑1 (MQ=11)
                                     ttGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAat                                                                                                         <  1:224667/141‑1 (MQ=18)
                                       tatatccGTTTCATCAATTCGCCACTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGGTATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       <  1:61956/135‑1 (MQ=11)
                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       >  1:119626/5‑141 (MQ=11)
                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       >  1:473646/5‑141 (MQ=11)
                                       gatatACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAAACCACGttatcataatct                                                                                                       >  2:126915/5‑131 (MQ=11)
                                       gatatACGTCTCATCAATAAGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATGGTTGAATGATAAACGTTTAAACAACGTTCTCTTAatat                                                                                                       >  2:361192/5‑141 (MQ=11)
                                       gacatACGTCTCACCAATACGCCATTTGTAATAAGCTTTTGTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGTCAT‑ATTCTTGATCCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAatat                                                                                                       <  2:19163/137‑1 (MQ=11)
                                       gATGTACGTCTCATCAATACACCATTTGAAATAAGCTTTTTTATGCTTTTTCTTGCAAATTTGATATAAAATTGGGGCAG‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAatat                                                                                                       >  1:430610/1‑141 (MQ=11)
                                           tACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGTCCGTTGAATGAGGAACGTTTACACCACGTTCCCTTAATATTTCa                                                                                                   <  2:112714/141‑1 (MQ=255)
                                            aCGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAg                                                                                                  <  2:638674/141‑1 (MQ=255)
                                             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAgt                                                                                                 >  1:277428/1‑140 (MQ=255)
                                             cGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGa                                                                                                 <  1:249611/141‑1 (MQ=255)
                                              gTCTGATCAATACGCCATTTGAAATAAGCTTTGTTATGCGTTTTATTCCAAATTTGATATAAAATTGCGGCAG‑CTTCTTGAACCCAACGGTAGATCGATGAGTGATGAACGTTCACACCGTGTACTCATAATCTTTCAgta                                                                                                >  2:83555/1‑139 (MQ=2)
                                              gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGat                                                                                                <  1:635427/141‑1 (MQ=255)
                                              gTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGat                                                                                                >  1:393923/1‑141 (MQ=255)
                                                ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatat                                                                                              <  1:54999/141‑1 (MQ=255)
                                                ctcATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatat                                                                                              >  2:614662/1‑141 (MQ=21)
                                                 tcaccaATATGCCATTTGTAATAAGCCTTTTTAGGCTTTTTCTTCCAAATTCGATATAAATATTGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGCTCTCTTAATATTTCAGatata                                                                                             <  1:344190/141‑1 (MQ=11)
                                                  catcgATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACAGTTGAATGATGAACGTTTACACCACGGTCCCTTAATATTTCAGatatat                                                                                            >  2:403309/1‑141 (MQ=255)
                                                  catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGatatat                                                                                            >  2:42801/1‑141 (MQ=255)
                                                  catcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGatatat                                                                                            <  1:213488/141‑1 (MQ=18)
                                                   atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTGGTCTAAATTTGTGACATAATT‑CGACATGATTCTTGAACCCAACGGTATAGCGTTGAATGATGAACGTTTACACCACGTCCCCTTAAATTTTCAGATATATc                                                                                           >  2:273991/1‑141 (MQ=14)
                                                   atcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATc                                                                                           >  2:148620/1‑141 (MQ=25)
                                                    tcaATACGCCATTTGTAATAGGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAACTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCATGTTCTCTTAATATTTCAGATATATCa                                                                                          >  1:402920/1‑141 (MQ=11)
                                                    tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCGCGTTCCCTTATTTTTTCAAatccagcc                                                                                          >  2:30690/1‑135 (MQ=11)
                                                    tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          <  1:32539/141‑1 (MQ=255)
                                                    tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          >  1:46622/1‑141 (MQ=255)
                                                    tcaATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCa                                                                                          >  2:649262/1‑141 (MQ=255)
                                                     caATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCAc                                                                                         <  2:280174/141‑1 (MQ=255)
                                                      aaTACGCCATTTATAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGAGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTAGACCACGTTCTCTTAATATTTCAGATATATCACg                                                                                        >  2:370939/1‑141 (MQ=11)
                                                       aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                       <  1:576129/141‑1 (MQ=255)
                                                       aTACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAAGTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCActa                                                                                       >  1:457574/1‑139 (MQ=255)
                                                       aTACACCATTTGTAATAAGCTTTTTTTGGCTTTTTCTTCCAAATTTGATATAAAATGGGGGCAT‑ATTCTTAAACCCAACGGTAGATCGTTGAATGATGAACGCTTACACCACGTTCTCTTAATATTTCAGATATATCACGa                                                                                       <  2:511633/141‑1 (MQ=11)
                                                         aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATa                                                                                     <  1:596854/141‑1 (MQ=255)
                                                         aCGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATa                                                                                     >  1:190011/1‑141 (MQ=255)
                                                          cGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATag                                                                                    >  1:386545/1‑140 (MQ=255)
                                                           gCCATTTGTAATAAGCTTTTTTAGGCGTTTTCTTCCAAATTTGATATAAGATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                   <  1:299116/141‑3 (MQ=11)
                                                             cATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTc                                                                                 >  2:355511/1‑141 (MQ=25)
                                                              aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATAACGATAGCTCa                                                                                >  1:129490/1‑141 (MQ=11)
                                                              aTTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCa                                                                                >  1:320481/1‑141 (MQ=25)
                                                               tttGTAATAATCTTTTCTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATCCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCaa                                                                               <  1:49771/141‑1 (MQ=11)
                                                                 tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATagc                                                                                   >  2:188025/1‑133 (MQ=255)
                                                                 tGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATg                                                                             <  2:580205/141‑1 (MQ=18)
                                                                   taataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGAGATATCACGATAACTCAATgca                                                                           <  2:170774/141‑3 (MQ=12)
                                                                               ttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGAtag                                                               <  1:100280/141‑1 (MQ=35)
                                                                                           ttcCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCt                                                   <  1:349945/141‑1 (MQ=35)
                                                                                             ccAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACCCAATGTATATCTTAGATAGTAGCCAACGGCTAc                                                 >  1:197381/1‑141 (MQ=25)
                                                                                             ccAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTAAAATATATCACGATAACTCAATGTATATCTTAGATGGTAGCCAACGGCTAc                                                 >  1:286692/1‑141 (MQ=255)
                                                                                              cAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACa                                                >  2:519093/1‑141 (MQ=35)
                                                                                                                                 cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa              <  1:37194/141‑1 (MQ=32)
                                                                                                                                       ggTAGAGCGTTGAATGATGAACGTTTACAGCCCGTGCCCTAAATATTTCAGATATATCACGATAACACAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATATCTGAAATAGt        >  1:489431/1‑141 (MQ=11)
                                                                                                                                             ccGTTGAATGATGAACGTTTACTCCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc  <  1:167627/141‑1 (MQ=17)
                                                                                                                                             |                                                                                                                                            
AATGGCACGATATAAATAGCTCCATTTTCCTTTTATTTTGATGTACGTCTCATCAATACGCCATTTGTAATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCAT‑ATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAC  >  CP000730/36273‑36553

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: ‑