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→ACGUSA300HOU_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.

GTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  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)
                                                                                          taGGGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAg                                            <  2:116986/140‑1 (MQ=255)
                                                                                          taGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAg                                            >  1:613858/1‑140 (MQ=255)
                                                                                           agagCGTAAATACATTCGATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGTTATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATTATGTTATTCAGAAATTTCGATGTACAATGACAGTCAgg                                           >  1:120056/1‑140 (MQ=14)
                                                                                            gagCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgt                                          >  1:534740/1‑140 (MQ=40)
                                                                                            gagCGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgt                                          >  2:518658/1‑140 (MQ=255)
                                                                                              tgttAAATACATTCAATACCTTTGAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttt                                        <  2:350967/137‑1 (MQ=14)
                                                                                               cGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgttta                                       >  1:300230/1‑140 (MQ=255)
                                                                                               cGTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGTCGCGTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgttta                                       >  1:129030/1‑140 (MQ=14)
                                                                                                gTAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaa                                      >  2:539978/1‑140 (MQ=255)
                                                                                                gTAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaa                                      <  1:647663/140‑1 (MQ=255)
                                                                                                 tAAATTCACTCAATACCTTTTAGATTAGTCTTTGTTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaag                                     <  2:360658/140‑1 (MQ=14)
                                                                                                 tAAATACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCCATGTACAATGACAGTCAGgtttaag                                     >  1:403225/1‑140 (MQ=255)
                                                                                                 tAAATACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGATGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaag                                     <  1:67175/140‑1 (MQ=255)
                                                                                                  aaaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagt                                    <  1:218443/140‑1 (MQ=255)
                                                                                                   aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtt                                   >  2:622599/1‑140 (MQ=255)
                                                                                                   aaTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagtt                                   >  2:393900/1‑140 (MQ=255)
                                                                                                    aTACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGgtttaagttt                                  <  1:525427/140‑1 (MQ=255)
                                                                                                     tACATTCAATACCTTTTAAAGTATTCTTTGCCGTGTTGATACTTTCGTATATTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttta                                 >  1:157335/1‑140 (MQ=14)
                                                                                                     tACATTCAATACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGTTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagttta                                 >  1:453491/1‑140 (MQ=255)
                                                                                                     tACACTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTTTTACTTTAATATTACGGTGATCTTGCTCAATGAGGTAAATCAGATATTACGATATACAATGGCATTCAGgtttaaattta                                 >  1:165056/1‑140 (MQ=14)
                                                                                                      aCATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGgtttaagtttaa                                >  1:625836/1‑140 (MQ=40)
                                                                                                       cactCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaa                               >  1:26686/4‑140 (MQ=255)
                                                                                                        actCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                              >  1:493417/3‑140 (MQ=255)
                                                                                                        actCAATACCTTTTAAAGTAGTCTTTGCTGTATTGATACTTTGATACCTTGTCTGTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTaaaa                              >  1:432231/3‑140 (MQ=255)
                                                                                                          tCAATACCTTTTAAAGTATTCTTTGACGTCTTTATACTTTGATATCTTGTCTTTCTTACTTTAATTTGACGGTGATCTTGCTCATAGATGTTATTCAGATATTTCTATGTACAATTAAAGACAGAGTTAAGTTTaaaaac                            >  2:335411/1‑138 (MQ=11)
                                                                                                          tCAATACATTTTAAAGTATTCTTTGCCGTATTGATACTTTGAGATCTTGTCTTTCTTACTATACTATGACGGTGACCTTTCTCAATGAGGGTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTTATTTTaaaatc                            >  2:254109/1‑138 (MQ=14)
                                                                                                          ccAATACCTTTTAAAGTATTCTTAGCAGTCTTGATATTATGATAGCCTGTATTTACTCCTTTAATGTGAAGGTGATCTAGCTCATACTGGTTATTCATATATTTCAATTTACAATGCAAACACGGTTTAAGTTTAAAAgt                            >  2:296844/2‑139 (MQ=2)
                                                                                                           cAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCCGATATTTCGATGTACAAATGCAGTGAGGTCTAAGTTTAAAAGCt                           <  1:365785/140‑1 (MQ=11)
                                                                                                            aaTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCtt                          >  2:185472/1‑140 (MQ=255)
                                                                                                             aTACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCttt                         >  1:382613/1‑140 (MQ=255)
                                                                                                             aTACCTTTTAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCttt                         >  2:103936/1‑140 (MQ=255)
                                                                                                              tatcTTTTAAAGTATTCTTTGCTGGATTGATACTTTTATTCCTTGTCGTTCTTACTTTAATATGACGGTGATCTTGTTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTa                        <  2:302571/137‑1 (MQ=14)
                                                                                                              tACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTa                        >  2:502923/1‑140 (MQ=255)
                                                                                                                ccTTTTAAATTATTCTTTGCCGTATTGATACTTTGATATATGGTTTTTCTTACTTTCATCTGACGGTGCTCTTGCTCAATCAGGTTATTCATATATTTCGATGTACAATTACAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      <  2:229369/140‑1 (MQ=11)
                                                                                                                ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTagt                      >  2:86896/1‑138 (MQ=255)
                                                                                                                ccTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      >  1:250437/1‑140 (MQ=255)
                                                                                                                ccTTTTAAAGTATCCTTTGCCGTATCGATACTTTGATCTCTTCTCTTTCTTACTTCAATATGACGGTGATCTTGCTTAATGACGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAt                      <  2:37029/140‑1 (MQ=14)
                                                                                                                 tttaaaaggTATTCTTTGCCGTATTGCTACTTTGATATCTTGTGTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAAtt                     >  1:631786/9‑140 (MQ=14)
                                                                                                                  ttttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    >  2:410901/1‑140 (MQ=40)
                                                                                                                  ttttAAAGTATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATTACGGTGTTCTTTCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTa                    >  1:327959/1‑140 (MQ=14)
                                                                                                                   tttAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTAc                   >  2:47392/1‑140 (MQ=255)
                                                                                                                    ttAAAGTATTCTTTGCTGTATTGATACTTTTATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACt                  <  1:584782/140‑1 (MQ=255)
                                                                                                                      aaaGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACttt                <  1:571966/140‑1 (MQ=255)
                                                                                                                        aGTATTCTTTGCTGTATTGATACGTTGATACGTTGTCTTTCTTACTTTAATATGACGGGGATCTTGCTCAATGTGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAg              <  2:213423/140‑1 (MQ=255)
                                                                                                                          tATTCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGcc            >  2:269435/1‑140 (MQ=255)
                                                                                                                             tCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACAGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCAtt         >  1:200072/1‑140 (MQ=255)
                                                                                                                             tCTTTGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGAGGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCAtt         <  2:728/140‑1 (MQ=255)
                                                                                                                               tttGCCGTATTGATACTTTGATATCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCATTGGGGTTATGCAGATATTTCGAGTTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTATTTTAGCCATTGc       <  1:319242/140‑1 (MQ=14)
                                                                                                                                ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      >  1:548246/1‑140 (MQ=255)
                                                                                                                                ttGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCt      <  2:449088/140‑1 (MQ=255)
                                                                                                                                 tGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTa     >  2:156886/1‑140 (MQ=255)
                                                                                                                                  gCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTCACCATCGCTAc    >  2:328770/1‑140 (MQ=255)
                                                                                                                                   cTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTAcc   <  2:461007/140‑1 (MQ=255)
                                                                                                                                    tGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGGCAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAACCATTGCTACCt  <  2:282388/140‑1 (MQ=255)
                                                                                                                                    |                                                                                                                                            
GTCAGTGTTCGCTTAACTTGCTAGCATGATGCTAATTTCGTGGCATGGCGAAAATCCGTAGATCTGAAGAGACCTGCGGTTCTTTTTATATAGAGCGTAAATACATTCAATACCTTTTAAAGTATTCTTTGCTGTATTGATACTTTGATACCTTGTCTTTCTTACTTTAATATGACGGTGATCTTGCTCAATGAGGTTATTCAGATATTTCGATGTACAATGACAGTCAGGTTTAAGTTTAAAAGCTTTAATTACTTTAGCCATTGCTACCT  >  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: