Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,479 T→C 22.9% T25A (ACA→GCA)  USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,4790TC22.9% 176.5 / 58.7 110T25A (ACA→GCA) USA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base T (46/38);  new base C (10/15);  total (56/54)
Fisher's exact test for biased strand distribution p-value = 2.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

Only 100 of 112 total aligned reads displayed.

AATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36341‑36616
                                                                                                                                           |                                                                                                                                         
aataaGCTTTTTTATGCTTTTTCTTCCAAATTTGATATAAAATTGGGGCATATTCTTGAACCCAACGGTAGATCGTTGAATGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTgaatgaat                                                                                                                                         >  1:693710/1‑141 (MQ=11)
       ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAg                                                                                                                                  >  1:661431/1‑141 (MQ=35)
       ttttttATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATCACGATAACTCAATGTATATCTTaa                                                                                                                                  >  1:321067/1‑140 (MQ=25)
             aTGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAACTCAATGTATATCTTAGAtagta                                                                                                                            >  1:191467/1‑141 (MQ=25)
               gCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGc                                                                                                                          >  2:1053063/1‑141 (MQ=35)
                  ttttCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCaa                                                                                                                       <  2:1004698/141‑1 (MQ=35)
                   tttCCTCCAAATTTGATACAAAATTGGTGCTTATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGccgtt                                                                                                                      <  1:35865/141‑4 (MQ=255)
                        tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTTGCCAACGGCTa                                                                                                                 >  1:586296/1‑141 (MQ=25)
                        tcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTa                                                                                                                 >  1:1025172/1‑141 (MQ=35)
                             aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGTTAGACGGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                            <  2:515535/141‑1 (MQ=21)
                             aTTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCAAGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTg                                                                                                            <  1:580186/141‑1 (MQ=25)
                               ttGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGAt                                                                                                          <  2:506800/141‑1 (MQ=35)
                                tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGCTGCCCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGAGAACTCAATGTATAGCTTAGATAATAGCCAACGGCTACCGTGATa                                                                                                         >  2:692055/1‑141 (MQ=11)
                                tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCACCGGCTCCAGTGATa                                                                                                         >  1:475285/1‑141 (MQ=21)
                                tGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATa                                                                                                         >  1:228455/1‑141 (MQ=35)
                                  aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAAc                                                                                                       >  1:716115/1‑141 (MQ=35)
                                     cAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACTTc                                                                                                    >  1:656641/1‑141 (MQ=17)
                                         aTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTAAAGGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTg                                                                                                <  1:728073/141‑1 (MQ=17)
                                          ttcgggaagagtCTTTAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                               <  1:87006/130‑1 (MQ=11)
                                          ttGGGGCATATTCTTTAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                               <  1:871889/141‑1 (MQ=17)
                                          ttGGGACATATTCTTTAACCCAACGCTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCCGATATATCACGTTGTCCCATTGTATGTCTGACATAGTATCCAACGGCTAGAGTGAGAACATCCTTGt                                                                                               <  2:40639/141‑1 (MQ=11)
                                           tGGGGCATATTCTTGAACCCAACTGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                              <  2:777507/141‑1 (MQ=17)
                                           tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                              <  1:147010/141‑1 (MQ=32)
                                             gggCATATTCTTGAACCCAACGGCAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTaa                                                                                            <  1:1035448/141‑1 (MQ=17)
                                               gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                          >  1:548503/1‑141 (MQ=32)
                                                 atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTg                                                                                        <  1:188848/141‑1 (MQ=32)
                                                 atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTCAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTCGtacaaacg                                                                                        >  2:630935/1‑134 (MQ=255)
                                                   atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                      <  1:709239/141‑1 (MQ=32)
                                                   atTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGtt                                                                                      <  1:773464/141‑1 (MQ=32)
                                                         gagccgaaCGGTAGACCGTTGAATGATGAACGTTTACACGACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATc                                                                                <  2:894835/135‑1 (MQ=255)
                                                            cccAACTGTAGACCGTTGAATGATTAACGTTTACACCACGTTCCCTTAATATTTCAGATATAGCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                             <  2:923171/141‑1 (MQ=12)
                                                            cccAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGa                                                                             >  2:906165/1‑141 (MQ=32)
                                                                aCGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATa                                                                         >  1:414258/1‑141 (MQ=32)
                                                                   ggtgACCGTTAAATGATTAACGGTTACACCACGTTCCCTTCATATCTCAGATATATCACGACACCTCAATGTATATCCTAGATAGTAGCCAACGGCTAAAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGtt                                                                      <  1:489085/138‑1 (MQ=11)
                                                                      gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                   >  1:1031138/1‑141 (MQ=32)
                                                                        ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGTTAAGGCAATCTCCATTTTATATAGTAGCCAACGGCTACAGTGATAACATCCGTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                 >  1:508680/1‑141 (MQ=11)
                                                                        ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                 <  1:66892/141‑1 (MQ=32)
                                                                         cGTTGAATTATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACa                                                                <  1:387279/141‑1 (MQ=17)
                                                                              aaggatgtACGGTTACACCACGTTCTCTCAGTACTTCAGTGATATCTCGCTAGCTATATGTATCGCTTGGATAGTAGCCAACGTCTACAGTGATGTCAGCTTTGTTAAATTGTTTATATCTGAAATAGTTCATGCagaagt                                                           <  2:467281/138‑2 (MQ=2)
                                                                              aaTGATGAACGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                           <  2:692814/141‑1 (MQ=17)
                                                                              aaTGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaaga                                                           <  1:312597/141‑1 (MQ=32)
                                                                                tgatgaACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACt                                                         >  2:742127/1‑141 (MQ=32)
                                                                                      aCGTTCACACCACGTTCGCTTAATATTACTGATATATCAGCCCAGTTAACGTTAGAGGTTAGTTAATATCGATCTGCTACACTGAGTACCTCCTTGTTAAATTGTTTATATCTGAAATATTTCATACAGAAAACTCCtttt                                                   >  2:382093/1‑141 (MQ=2)
                                                                                       cGTTTACACCACGTACCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGGAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGATGAATCCttttt                                                  >  2:461282/1‑141 (MQ=255)
                                                                                         tttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGt                                                <  2:847029/141‑1 (MQ=32)
                                                                                          ttACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAACCGTTCAGACGGAAGACTCCTTTTTGtt                                               <  1:1053063/141‑1 (MQ=255)
                                                                                          tcacatcgcGTTCCCTTTATATTTCAGATATATCACGATAACTCTATGGATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGtt                                               <  2:901211/133‑1 (MQ=1)
                                                                                           tACAGCACGTTCCCTTAATATTTCAGATATGTCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGttc                                              >  2:646080/1‑140 (MQ=17)
                                                                                           tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTa                                              >  2:449144/1‑141 (MQ=32)
                                                                                              accccGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                           <  2:581429/141‑1 (MQ=17)
                                                                                              accacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaa                                           >  2:284929/1‑141 (MQ=32)
                                                                                               ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaaaag                                          >  2:645991/1‑140 (MQ=17)
                                                                                               ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAt                                          >  2:569891/1‑141 (MQ=32)
                                                                                               ccacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGATACAGTGATAACTTCATTGTTACATTGTTTATATATGGGTTAGTTCAGGCGGAATAGACGTTATTGTTAAAAt                                          >  2:728509/1‑141 (MQ=11)
                                                                                                cacGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAtt                                         <  1:1003147/141‑1 (MQ=32)
                                                                                                  cGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGATTCCtttttttttttttttt                                       >  2:496266/1‑133 (MQ=255)
                                                                                                    ttgccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAc                                     <  2:296580/138‑1 (MQ=17)
                                                                                                      cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACta                                   <  2:985869/141‑1 (MQ=32)
                                                                                                      cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTGTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAATTTTTTACtt                                   >  2:703258/1‑140 (MQ=11)
                                                                                                         ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                                >  2:955833/1‑141 (MQ=33)
                                                                                                                 tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                        >  2:379504/1‑141 (MQ=35)
                                                                                                                  cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                       >  1:344295/1‑141 (MQ=35)
                                                                                                                  cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                       >  2:323225/1‑141 (MQ=35)
                                                                                                                  cAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACtt                       >  2:188940/1‑141 (MQ=35)
                                                                                                                      tatatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                   >  2:880522/1‑141 (MQ=35)
                                                                                                                      aatatTACGATAACTCAATGTATATCTTAGATAGTAGACAACGGCTAGAGTAATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCa                   <  1:248578/140‑1 (MQ=255)
                                                                                                                       atatCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCaa                  >  1:156417/1‑141 (MQ=25)
                                                                                                                          tCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCACCAg               >  2:644793/1‑141 (MQ=25)
                                                                                                                           cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa              >  2:289140/1‑141 (MQ=34)
                                                                                                                              gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATAAGACTCTTTTTTGTTAAAATTATACTATAATTTCAACTTTGCAACAGAAcc           <  1:646080/141‑1 (MQ=255)
                                                                                                                              gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc           <  1:449144/141‑1 (MQ=35)
                                                                                                                              gATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcc           <  1:418981/141‑1 (MQ=34)
                                                                                                                               aTAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGAGAACAGCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGACGCCTCCGTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg          <  1:205571/141‑1 (MQ=11)
                                                                                                                                tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATCGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTCTAACTTCAACTTTGCAACGGAACCGt         <  2:619899/138‑1 (MQ=11)
                                                                                                                                tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:520365/138‑1 (MQ=9)
                                                                                                                                tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt         <  1:427652/138‑1 (MQ=9)
                                                                                                                                  gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGCAATAGGTCATACCGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtct       >  2:519792/2‑139 (MQ=11)
                                                                                                                                  gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAACCGtat       >  2:345531/2‑141 (MQ=12)
                                                                                                                                  gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat       <  1:1021151/140‑1 (MQ=16)
                                                                                                                                  gcTCAATGCATATCTTAGATAGTAACCAAAGGCTACAGTGATAACATCCTTGTTAAATTGTATATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATATTGTAACTTCAACTTTGCGACAGAACCGtat       >  2:569962/2‑141 (MQ=11)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAACCGtatg      <  2:650714/141‑2 (MQ=255)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  2:956179/1‑141 (MQ=16)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  1:537002/1‑141 (MQ=16)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  2:823542/1‑141 (MQ=16)
                                                                                                                                   cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt      >  2:111099/1‑141 (MQ=16)
                                                                                                                                    tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta     >  1:73968/1‑141 (MQ=40)
                                                                                                                                    tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta     <  1:906165/141‑1 (MQ=40)
                                                                                                                                    tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta     >  1:1011714/1‑141 (MQ=40)
                                                                                                                                    tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGAAACAGAACCGtatta     >  1:810811/1‑141 (MQ=255)
                                                                                                                                    tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAACCGtattt     <  2:1052761/141‑2 (MQ=255)
                                                                                                                                    tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt     <  1:223217/141‑2 (MQ=12)
                                                                                                                                    tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt     <  1:106745/141‑2 (MQ=16)
                                                                                                                                    tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt     <  1:722059/141‑2 (MQ=16)
                                                                                                                                    tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatgt     >  1:47710/1‑139 (MQ=12)
                                                                                                                                    tCAATGCATATCTTAAATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt     >  1:787303/1‑140 (MQ=255)
                                                                                                                                     cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:52945/141‑3 (MQ=255)
                                                                                                                                     cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  1:575929/141‑3 (MQ=255)
                                                                                                                                     cAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttt    <  2:192865/141‑3 (MQ=255)
                                                                                                                                      aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGTCTTTTg   >  1:1033368/1‑141 (MQ=255)
                                                                                                                                       aTGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAAAAGAACCGTATTATgg  >  2:463064/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                         
AATAAGCTTTTTTATGCTTTTTCTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGG  >  CP000730/36341‑36616

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