Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 35,997 T→C 28.0% T185T (ACA→ACG)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073035,9970TC28.0% 704.5 / inf 534T185T (ACA→ACG) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (203/178);  new base C (80/68);  total (287/247)
Fisher's exact test for biased strand distribution p-value = 9.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 536 total aligned reads displayed.

GTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35865‑36136
                                                                                                                                    |                                                                                                                                            
   aGTGTTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAc                                                                                                                                    >  1:18135/1‑140 (MQ=11)
   aGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTAGATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATAc                                                                                                                                    <  1:656331/140‑1 (MQ=2)
    gtgtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGAGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACt                                                                                                                                   >  1:284000/1‑140 (MQ=11)
     tgtTCGCTTAACTTGCTAGCATGCTGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACATTTTAAAGTATTCTTTGCTGTATTGATACtt                                                                                                                                  >  1:210348/1‑140 (MQ=12)
      gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGGTAAAATCCGTAGATCTTAAGAGACCTACGGTTCTTTTTATATAGAGCGGAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACttt                                                                                                                                 <  2:649666/140‑1 (MQ=11)
      gtTCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                                                 <  2:548954/140‑1 (MQ=16)
      gtTCGCTTAACTTGCTTGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACttt                                                                                                                                 <  1:80714/140‑1 (MQ=16)
        tCGTTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGa                                                                                                                               <  2:628678/140‑1 (MQ=16)
           ctcgtattGCGAGCATGATGCGCATTTCGTGGCAGGGCGAAGATTCGTAGATCTGATGAGTCCTGCGGTTTTTGTTATATAGAGCGTAAATATGCTCAATAACTTTTAGGATGTGCTTCGCTTTATTGTTGCTTTGATAc                                                                                                                            <  2:276762/134‑1 (MQ=2)
                cTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGt                                                                                                                       <  1:373544/140‑1 (MQ=18)
                 ttGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGCCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtc                                                                                                                      <  2:614156/140‑1 (MQ=18)
                   gCTAGCAGGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctt                                                                                                                    >  1:420186/1‑140 (MQ=18)
                     tAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGtctttc                                                                                                                  <  1:446172/140‑1 (MQ=18)
                        cATGATGCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATAACTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTa                                                                                                               >  1:276510/1‑140 (MQ=11)
                          tgatgCTAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCGGTATTGATACTTTGATAGCtttcataagtgacg                                                                                                             >  2:136003/1‑128 (MQ=11)
                            atgCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACttt                                                                                                           >  1:558252/1‑140 (MQ=35)
                                tAATTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATAATTTAATATCTTGTCTTTCTTACTTTAatt                                                                                                       >  1:600465/1‑139 (MQ=11)
                                  aTTTCGTGGCATGGGGAAAATCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATAACTTGTCTTTCTTACTTTAATATg                                                                                                     >  2:57899/1‑140 (MQ=255)
                                  aTTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATg                                                                                                     >  1:141332/1‑140 (MQ=21)
                                   tttCGTGGCATGTCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTATATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGa                                                                                                    <  2:593538/140‑1 (MQ=255)
                                    ttCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGAc                                                                                                   <  2:43838/140‑1 (MQ=37)
                                       gTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCCGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAATATTCTTTGCTGTATTGATACTTTGATAACTTGTCTTTCTTACTGTAATATGACGGt                                                                                                >  1:114632/1‑140 (MQ=255)
                                        tGGCATGGCGAAAATCCGTAGATCTGAAGAGATTTGCGGATCTTATCATGTAGTTCGGACTTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATACCGGTg                                                                                               >  2:639685/1‑140 (MQ=11)
                                         ggCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCCTTCTTACTTTAATCTAACGGTGa                                                                                              >  2:109051/1‑140 (MQ=11)
                                              ggCGAAAATCCGTAGATCTGAAGATACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTTAATACCTTTTAAAGTATTCTTTGCTGTATTGAGACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTg                                                                                         <  1:102097/140‑1 (MQ=255)
                                              ggCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGCTACCTTGTCTTTCTTACTTTAATATGACGGTCATCTTg                                                                                         >  2:609582/1‑140 (MQ=255)
                                                 gAAAATCCGTAGATTTGATGATTCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATTCCTTTTAAGGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTc                                                                                      <  1:312505/140‑1 (MQ=14)
                                                   aaaTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTCGCTCaa                                                                                    <  1:499595/140‑1 (MQ=255)
                                                    aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAAt                                                                                   <  1:240639/140‑1 (MQ=39)
                                                    aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCact                                                                                   >  2:565920/1‑138 (MQ=255)
                                                    aaTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCCTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATCTGACGGTGATCTTGCTCAAt                                                                                   <  1:468235/140‑1 (MQ=255)
                                                     aTCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTCAATATGACGGTGATCTTGCTCAATg                                                                                  >  1:407238/1‑140 (MQ=255)
                                                     aTCCGTAGATCTGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCCATg                                                                                  >  2:238670/1‑140 (MQ=14)
                                                        cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGATGATCTTGCTCAATAAgg                                                                               >  2:167978/1‑140 (MQ=255)
                                                        cGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTAATTTAATATGAAGGTGGTCTTGCTCAATGAgg                                                                               >  2:196154/1‑140 (MQ=255)
                                                                tGAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                                                       >  2:204934/1‑140 (MQ=255)
                                                                tGAAGACACCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAg                                                                       <  2:381644/140‑1 (MQ=255)
                                                                 gAAGAGACCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGa                                                                      <  2:345428/140‑1 (MQ=255)
                                                                    gagaCCTACGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCGTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGatat                                                                   <  1:588946/140‑1 (MQ=255)
                                                                     agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACGCTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                                                  <  2:122629/140‑1 (MQ=255)
                                                                     agaCCTGCGGTTCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATAtt                                                                  <  1:393980/140‑1 (MQ=38)
                                                                                tCTTTTTATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGGACaa                                                       >  2:492450/1‑140 (MQ=255)
                                                                                   ttttATATAGAGCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATgg                                                    <  1:337020/140‑2 (MQ=255)
                                                                                    tttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGAc                                                   <  2:263556/140‑1 (MQ=255)
                                                                                     ttATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTAGATGTACAATGACa                                                  >  1:106470/1‑140 (MQ=255)
                                                                                        tataGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTAAGATATTTCGATGTAAAATGACAGTc                                               >  2:423473/1‑140 (MQ=255)
                                                                                          taGGGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑g                                            <  2:116986/140‑1 (MQ=255)
                                                                                          taGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑g                                            >  1:613858/1‑140 (MQ=255)
                                                                                           agagCGTAAATACATTCGATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGTTATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATTATGTTATTCAGAAATTTCGATGTACAATGACAGTCA‑gg                                           >  1:120056/1‑140 (MQ=14)
                                                                                            gagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggt                                          >  1:534740/1‑140 (MQ=40)
                                                                                            gagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggt                                          >  2:518658/1‑140 (MQ=255)
                                                                                              tgttAAATACATTCAATACCTTTGAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggttt                                        <  2:350967/137‑1 (MQ=14)
                                                                                               cGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggttta                                       >  1:300230/1‑140 (MQ=255)
                                                                                               cGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGTCGCGTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggttta                                       >  1:129030/1‑140 (MQ=14)
                                                                                                gTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggtttaa                                      >  2:539978/1‑140 (MQ=255)
                                                                                                gTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaa                                      <  1:647663/140‑1 (MQ=255)
                                                                                                 tAAATTCACTCAATACCTTTTAGATTAGTCTTTGTTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaag                                     <  2:360658/140‑1 (MQ=14)
                                                                                                 tAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCCATGTACAATGACAGTCA‑Ggtttaag                                     >  1:403225/1‑140 (MQ=255)
                                                                                                 tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGATGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaag                                     <  1:67175/140‑1 (MQ=255)
                                                                                                  aaaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaagt                                    <  1:218443/140‑1 (MQ=255)
                                                                                                   aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaagtt                                   >  2:622599/1‑140 (MQ=255)
                                                                                                   aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaagtt                                   >  2:393900/1‑140 (MQ=255)
                                                                                                    aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑Ggtttaagttt                                  <  1:525427/140‑1 (MQ=255)
                                                                                                     tACATTCAATACCTTTTAAAGTATTCTTTGCCGTGTTGATACTTTCGTATATTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggtttaagttta                                 >  1:157335/1‑140 (MQ=14)
                                                                                                     tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGTTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggtttaagttta                                 >  1:453491/1‑140 (MQ=255)
                                                                                                     tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTTTTACTTTAATATTACGGTGATCTTGCTCAATGAGGTAAATCAGATATTACGATATACAATGGCATTCA‑Ggtttaaattta                                 >  1:165056/1‑140 (MQ=14)
                                                                                                      aCATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑Ggtttaagtttaa                                >  1:625836/1‑140 (MQ=40)
                                                                                                       cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTaaa                               >  1:26686/4‑140 (MQ=255)
                                                                                                        actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTaaaa                              >  1:493417/3‑140 (MQ=255)
                                                                                                        actCAATACCTTTTAAAGTAGTCTTTGCTGTATTGATACTTTGATACCTTGTCTGTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTaaaa                              >  1:432231/3‑140 (MQ=255)
                                                                                                          tCAATACCTTTTAAAGTATTCTTTGACGTCTTTATACTTTGATATCTTGTCTTTCTTACTTTAATTTGACGGTGATCTTGCTCATAGATGTTATTCAGATATTTCTATGTACAATTAAAGACA‑GAGTTAAGTTTaaaaac                            >  2:335411/1‑138 (MQ=11)
                                                                                                          tCAATACATTTTAAAGTATTCTTTGCCGTATTGATACTTTGAGATCTTGTCTTTCTTACTATACTATGACGGTGACCTTTCTCAATGAGGGTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTTATTTTaaaatc                            >  2:254109/1‑138 (MQ=14)
                                                                                                          ccAATACCTTTTAAAGTATTCTTAGCAGTCTTGATATTATGATAGCCTGTATTTACTCCTTTAATGTGAAGGTGATCTAGCTCATACTGGTTATTCATATATTTCAATTTACAATG‑CAAACACGGTTTAAGTTTAAAAgt                            >  2:296844/2‑139 (MQ=2)
                                                                                                           cAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCCGATATTTCGATGTACAAATGCAGTGA‑GGTCTAAGTTTAAAAGCt                           <  1:365785/140‑1 (MQ=11)
                                                                                                            aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCtt                          >  2:185472/1‑140 (MQ=255)
                                                                                                             aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCttt                         >  1:382613/1‑140 (MQ=255)
                                                                                                             aTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCttt                         >  2:103936/1‑140 (MQ=255)
                                                                                                              tatcTTTTAAAGTATTCTTTGCTGGATTGATACTTTTATTCCTTGTCGTTCTTACTTTAATATGACGGTGATCTTGTTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTa                        <  2:302571/137‑1 (MQ=14)
                                                                                                              tACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTa                        >  2:502923/1‑140 (MQ=255)
                                                                                                                ccTTTTAAATTATTCTTTGCCGTATTGATACTTTGATATATGGTTTTTCTTACTTTCATCTGACGGTGCTCTTGCTCAATCAGGTTATTCATATATTTCGATGTACAATTACAGTCA‑GGTTTAAGTTTAAAAGCTTTAAt                      <  2:229369/140‑1 (MQ=11)
                                                                                                                ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTagt                      >  2:86896/1‑138 (MQ=255)
                                                                                                                ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAAt                      >  1:250437/1‑140 (MQ=255)
                                                                                                                ccTTTTAAAGTATCCTTTGCCGTATCGATACTTTGATCTCTTCTCTTTCTTACTTCAATATGACGGTGATCTTGCTTAATGACGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAAt                      <  2:37029/140‑1 (MQ=14)
                                                                                                                 tttaaaaggTATTCTTTGCCGTATTGCTACTTTGATATCTTGTGTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAAtt                     >  1:631786/9‑140 (MQ=14)
                                                                                                                  ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTa                    >  2:410901/1‑140 (MQ=40)
                                                                                                                  ttttAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATTACGGTGTTCTTTCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTa                    >  1:327959/1‑140 (MQ=14)
                                                                                                                   tttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTAc                   >  2:47392/1‑140 (MQ=255)
                                                                                                                    ttAAAGTATTCTTTGCTGTATTGATACTTTTATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACt                  <  1:584782/140‑1 (MQ=255)
                                                                                                                      aaaGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACttt                <  1:571966/140‑1 (MQ=255)
                                                                                                                        aGTATTCTTTGCTGTATTGATACGTTGATACGTTGTCTTTCTTACTTTAATATGACGGGGATCTTGCTCAATGTGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAg              <  2:213423/140‑1 (MQ=255)
                                                                                                                          tATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTATTTTAGcc            >  2:269435/1‑140 (MQ=255)
                                                                                                                             tCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACAGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAtt         >  1:200072/1‑140 (MQ=255)
                                                                                                                             tCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGAGGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCAtt         <  2:728/140‑1 (MQ=255)
                                                                                                                               tttGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCATTGGGGTTATGCAGATATTTCGAGTTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGc       <  1:319242/140‑1 (MQ=14)
                                                                                                                                ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      >  1:548246/1‑140 (MQ=255)
                                                                                                                                ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      <  2:449088/140‑1 (MQ=255)
                                                                                                                                 tGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTa     >  2:156886/1‑140 (MQ=255)
                                                                                                                                  gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTCACCATCGCTAc    >  2:328770/1‑140 (MQ=255)
                                                                                                                                   cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   <  2:461007/140‑1 (MQ=255)
                                                                                                                                    tGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCt  <  2:282388/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                            
GTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCA‑GGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  CP000730/35865‑36136

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