Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,814 Δ1 bp 41.1% intergenic (+45/+63) USA300HOU_0025 → / ← USA300HOU_0026 hypothetical protein/IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,8140A.41.1% 113.6 / 118.7 73intergenic (+45/+63)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (22/19);  new base . (10/20);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 9.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

ATGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCT  >  CP000730/35675‑35947
                                                                                                                                           |                                                                                                                                     
atGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACTAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATa                                                                                                                                       >  2:356554/1‑140 (MQ=255)
   cAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGaa                                                                                                                                    >  1:51706/1‑140 (MQ=255)
         attCTGAATTTAAAGATGTAGATGATCTCAGTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGTTATAATAAAAATAGAAGGTTCt                                                                                                                              <  2:486831/140‑1 (MQ=255)
           tCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGt                                                                                                                            >  2:254251/1‑140 (MQ=255)
           tCTGAATTTAAAGATGCAGCAGAGCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATATAAGGTTCTGt                                                                                                                            >  2:307388/1‑140 (MQ=255)
            cTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGtt                                                                                                                           >  1:306120/1‑140 (MQ=255)
            cTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTAATTGAGGAACGGCGAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGtt                                                                                                                           >  1:216147/1‑140 (MQ=255)
             tGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTg                                                                                                                          <  2:561812/140‑1 (MQ=255)
              gAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGc                                                                                                                         >  1:57776/1‑140 (MQ=255)
              gAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGc                                                                                                                         <  1:407308/140‑1 (MQ=255)
                   tAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGt                                                                                                                    >  2:441791/1‑140 (MQ=255)
                   gaaaGATGTAGATGATCTCATAGGCAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGAGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGGTCTGATAAAAATAGAAGGATATGTTGACAAGt                                                                                                                    >  2:407481/2‑140 (MQ=255)
                       gATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTAAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTaaaa                                                                                                                <  1:637113/140‑1 (MQ=255)
                          gTAGATGATCTCATTGACAAAGTCTTTAATGTAGAAGCTATTCAGTATAAAGTTTTAGCCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATATAAGGTTCTGTTGTAAAGTaaaaaaa                                                                                                             >  1:58819/1‑140 (MQ=255)
                                tagctcATTGACAAAGTATTTAATGAAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGAGTACTTTAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTACAAAGTAAAAAAATATAGc                                                                                                       <  2:273041/137‑1 (MQ=255)
                                      aTTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATAGTGATCTAATAAAAATAGAAGGATCTGTTGAAAAGTAAAAAAATATAGCtcccca                                                                                                 >  2:115771/1‑135 (MQ=255)
                                               gTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTAt                                                                                        >  2:251520/1‑140 (MQ=255)
                                                 aTTTAATGTTGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGCTTCAAAGTAAAAAAATATAGCTAACCACTAATTTATCa                                                                                      <  1:627828/140‑1 (MQ=255)
                                                     aaTGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATATc                                                                                  >  2:147723/1‑140 (MQ=255)
                                                             aGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATGAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCg                                                                          >  2:236940/1‑140 (MQ=255)
                                                              gCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAAAAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGc                                                                         >  1:518698/1‑140 (MQ=255)
                                                                 aTTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTa                                                                      >  1:349692/1‑140 (MQ=255)
                                                                   tCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAAc                                                                    >  2:635511/1‑140 (MQ=255)
                                                                    cAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACt                                                                   <  2:307882/140‑1 (MQ=255)
                                                                          aacgTTTTTGGCGAAAAATGGTTACTTGAACAAGAGCTAAGATAAAATTGTGATCTAATAAAAATAGAATGTTCTGTTGCAAAGTAAACACATCTAGCTAACCACCAATTTATCACGGCATTGTTCGCTTAATTGGCTAg                                                             <  2:557235/137‑1 (MQ=255)
                                                                                    gCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCACGATGCTa                                                   >  1:336889/1‑140 (MQ=255)
                                                                                     cGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTaa                                                  >  2:84775/1‑140 (MQ=255)
                                                                                         aaaTGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAAtttg                                              >  1:213423/1‑139 (MQ=255)
                                                                                             gATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTgg                                          >  2:169961/1‑140 (MQ=255)
                                                                                                tACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATTGAAGGTTCTGTTGCTTAGTATATTAATATAGCTAGCCACTGATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCAt                                       <  1:280684/140‑1 (MQ=21)
                                                                                                   ttGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGCCAGTGTTCGCTTAACTTGCTAGCATGACGCTAATTTCGTGGCATGGc                                    <  2:334978/140‑1 (MQ=255)
                                                                                                        cAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGaaaa                               <  2:51706/140‑1 (MQ=255)
                                                                                                              cTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGATATGTTTGTGGCATGGCGAAAATCCGTa                         <  2:306120/140‑1 (MQ=255)
                                                                                                              cTAAAAAAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTa                         <  2:216147/140‑1 (MQ=255)
                                                                                                               taaaataaaaTTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAg                        >  1:663515/1‑140 (MQ=255)
                                                                                                                   acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTACAAAATCTAGCTAATCACTAATATATCATACCAATATTCGTTTAACTGGTGAGCATGATGCTAATTTCGTGGCATGGCGAAAATGCGTAGATCTGa                  >  1:288880/15‑140 (MQ=2)
                                                                                                                   acaaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGa                  >  1:20869/15‑140 (MQ=11)
                                                                                                                    caaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGaa                 <  1:265277/127‑1 (MQ=11)
                                                                                                                     aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCTGTGTTTGTTTATCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                <  1:194225/128‑1 (MQ=11)
                                                                                                                     aaaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAag                <  1:628378/128‑1 (MQ=11)
                                                                                                                      aaagtgtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCTCTAATATATCTTGTCAGTGTTCGTTTATCTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga               <  1:74933/129‑1 (MQ=11)
                                                                                                                      aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga               <  1:517866/129‑1 (MQ=11)
                                                                                                                      aaagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga               >  2:159932/12‑140 (MQ=11)
                                                                                                                       aagagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAagag              <  1:534692/130‑1 (MQ=11)
                                                                                                                        atTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAaga               >  2:584782/1‑140 (MQ=255)
                                                                                                                        agagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAgaga             >  2:296117/10‑140 (MQ=11)
                                                                                                                         gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc            <  1:123561/132‑1 (MQ=11)
                                                                                                                         gagtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAc            <  1:83130/132‑1 (MQ=11)
                                                                                                                          agtattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc           <  1:237830/133‑1 (MQ=11)
                                                                                                                          agtattttaataaAAATGATGGCTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc           <  2:294229/133‑1 (MQ=11)
                                                                                                                            tgactatataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATACCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc           <  2:307935/133‑1 (MQ=255)
                                                                                                                            tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc           <  2:630285/140‑1 (MQ=255)
                                                                                                                            tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc           <  1:504728/140‑1 (MQ=255)
                                                                                                                            tgATCTAATAAAAATAGAAGGTTCTGTTACAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc           <  1:384811/140‑1 (MQ=255)
                                                                                                                            tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGCAgagagttg         <  2:576400/135‑5 (MQ=11)
                                                                                                                            tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg         >  1:545401/6‑140 (MQ=11)
                                                                                                                            tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTg         <  2:357206/135‑1 (MQ=11)
                                                                                                                            tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAGCCGTAGATCTGAAGAGACCTg         <  1:34158/135‑1 (MQ=11)
                                                                                                                            tattttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGATAGCATGCTGCTCAGTTCGTGGCATGGCGATAATCCGTAAATCTGCAGAGAcacg         >  2:329856/6‑137 (MQ=11)
                                                                                                                             gATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTATATCTGAAGAGACCt          <  1:356554/140‑1 (MQ=255)
                                                                                                                             attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGGATGGCGAAAATCCGTAGATCTGAAGAGACCTGc        <  1:158799/136‑1 (MQ=11)
                                                                                                                             attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc        >  1:526189/5‑140 (MQ=11)
                                                                                                                             attttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGACAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAAGCTGTAGATCTGAAGAGAACTGc        >  1:118936/5‑140 (MQ=11)
                                                                                                                             attttaataaAAAGGATGGCTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATAGATCATGTCAGTGTTCGTTTAACTTGCTAGCATGACGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGc        <  2:32631/136‑1 (MQ=11)
                                                                                                                              ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTACTATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg       >  1:46612/4‑140 (MQ=11)
                                                                                                                              ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg       <  2:9117/137‑1 (MQ=11)
                                                                                                                              ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg       <  2:14986/137‑1 (MQ=11)
                                                                                                                              ttttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg       <  1:665108/137‑1 (MQ=11)
                                                                                                                                ttaataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTTATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGt     <  1:595201/139‑1 (MQ=11)
                                                                                                                                cTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCg       <  1:635511/140‑1 (MQ=255)
                                                                                                                                 taataaAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTCGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCgg      <  2:354634/140‑1 (MQ=255)
                                                                                                                                 taataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGtt    <  2:45730/140‑1 (MQ=11)
                                                                                                                                  aataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAATATTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTc   >  1:228945/1‑140 (MQ=11)
                                                                                                                                   ataaAAATGATGGTTCTGTTGCAAAGTAAAAAATATAGCTAATCACTAATATATCATGTCAGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCt  <  1:602248/140‑1 (MQ=11)
                                                                                                                                           |                                                                                                                                     
ATGCAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCT  >  CP000730/35675‑35947

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