Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,817 A→T 42.5% intergenic (+48/+60) 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,8170AT42.5% 75.1 / 89.5 73intergenic (+48/+60)USA300HOU_0025/USA300HOU_0026hypothetical protein/IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base T (10/21);  total (31/42)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.41e-02

CAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTAT  >  CP000730/35678‑35953
                                                                                                                                           |                                                                                                                                        
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                 <  1:504728/140‑1 (MQ=255)
                                                                                                                         tgATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGAcc                 <  2:630285/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)
                                                                                                                                        tgAAGGTTCTGTTGAAAAGCAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCttttcgt  <  2:504480/139‑4 (MQ=9)
                                                                                                                                           |                                                                                                                                        
CAAATTATTCTGAATTTAAAGATGTAGATGATCTCATTGACAAAGTATTTAATGTAGAAGCTATTCAGTATAAAGTTTTTGGCGAAAAATGATTACTTGAACAAGAGCTAAAATAAAATTGTGATCTAATAAAAATAGAAGGTTCTGTTGCAAAGTAAAAAAATATAGCTAACCACTAATTTATCATGTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTAT  >  CP000730/35678‑35953

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: