Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,501 G→A 5.2% A17A (GCC→GCT USA300HOU_0026 ← IS431mec transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,5010GA 5.2% 410.2 / 9.9 175A17A (GCC→GCTUSA300HOU_0026IS431mec transposase
Reads supporting (aligned to +/- strand):  ref base G (76/88);  new base A (5/4);  total (82/93)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 175 total aligned reads displayed.

CTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATT  >  CP000730/36363‑36637
                                                                                                                                           |                                                                                                                                        
cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                        <  2:378390/141‑1 (MQ=35)
cttcCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                        <  2:536554/141‑1 (MQ=35)
cttcCAAAGTTGATACAAAAGTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGc                                                                                                                                        <  1:277769/141‑1 (MQ=25)
            aTACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATGTCTCACGATAACTCAATGTATATCTTAGATAGTAGCCAGCGGCTACAGTGATaaa                                                                                                                            >  1:26528/1‑140 (MQ=255)
                  aaTTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCtt                                                                                                                      <  1:883063/141‑1 (MQ=32)
                    ttGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGt                                                                                                                    <  2:714914/141‑1 (MQ=32)
                     tGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCCCGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGtt                                                                                                                   <  2:288636/141‑1 (MQ=17)
                      ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATCACTCAATGTATATCTTAGATAGTAGCCAACGGCTCCAGTGATAACATCCTTGTTa                                                                                                                  >  2:48257/1‑141 (MQ=12)
                      ggggCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATGACATCCTTGTTa                                                                                                                  >  2:134463/1‑141 (MQ=17)
                      ggggCATATTCTTGAACCCAACGGTAGACCGCTGAATGATGAAAGTTTACGCAACGCAACTTTAATACTTTTGTTATATGCTGATAAATTGATCTAAAGCTTAAGGGGTAGCAAACTGCTACAGTTTTAACATCCTTGTTa                                                                                                                  >  2:785553/1‑141 (MQ=2)
                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                               >  2:653764/1‑141 (MQ=32)
                         gCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAt                                                                                                               >  2:880897/1‑141 (MQ=32)
                           atatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTg                                                                                                             <  1:286313/141‑1 (MQ=17)
                            tatTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGt                                                                                                            <  2:27748/141‑1 (MQ=32)
                                 ttGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTtata                                                                                                       >  2:571846/1‑141 (MQ=32)
                                  tGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATCTTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTatat                                                                                                      >  1:153761/1‑141 (MQ=17)
                                   gAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTGTATc                                                                                                     >  2:894903/1‑141 (MQ=17)
                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                                    <  1:512041/141‑1 (MQ=32)
                                    aaCCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCt                                                                                                    <  2:879713/141‑1 (MQ=32)
                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            <  1:530140/141‑1 (MQ=32)
                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            <  1:462828/141‑1 (MQ=32)
                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            >  2:291693/1‑141 (MQ=32)
                                            ggTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAAATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAAAGGCTACAGTGATAACATACTTGTTAAATTGTTTATATCTGAAATAGt                                                                                            >  2:224637/1‑141 (MQ=12)
                                              tAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTc                                                                                          >  1:806273/1‑141 (MQ=32)
                                               aGACCGTTGAATGAAGAACGTTTACACCACGTTCCCTTAATATTTCAGATATAACACGGTAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAGTGTTTATATCTGAAATAGTTCa                                                                                         <  1:659888/141‑1 (MQ=255)
                                                gACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCAt                                                                                        >  1:755353/1‑141 (MQ=32)
                                                  ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATAc                                                                                      >  2:65562/1‑141 (MQ=32)
                                                  ccGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGATACAGTCATACCATCCTTGTTAAATTTGTTATATTTTAAATAGTTCATAc                                                                                      >  2:267128/1‑141 (MQ=11)
                                                      tGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACagaa                                                                                  >  2:100273/1‑141 (MQ=32)
                                                               taCGTTTACACCACGTTCTCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                                         <  2:411153/140‑1 (MQ=17)
                                                               aaCTTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                                         <  2:80749/141‑1 (MQ=17)
                                                               aaCGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttt                                                                         <  2:819449/141‑1 (MQ=32)
                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGATGTTGTTCAGACAGAAGACTCCttttt                                                                       <  2:669330/141‑1 (MQ=255)
                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCttttt                                                                       <  2:73262/141‑1 (MQ=32)
                                                                 cGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATAATTTAGATAGTAGCCAACGGCTACAGTGATAACATCATTGTTAAATTGTTTATATCTGAAATAGTTCATACAAAAGCATCCttttt                                                                       >  2:124551/1‑141 (MQ=11)
                                                                  gTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTg                                                                      >  1:18119/1‑141 (MQ=32)
                                                                     tACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACAACCAAGTTAAATTGTTTATATCTGAACGAGTTCGTACAGAAAACTCCTTTTTGTTa                                                                   >  2:150523/1‑141 (MQ=11)
                                                                      acacCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTaa                                                                  >  1:858728/1‑141 (MQ=32)
                                                                            cgttatcTTAATATTACGGATATATCACGAAAGATCAGTTGTAAACTTAGATTGTGGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATtat                                                            >  1:452373/7‑141 (MQ=11)
                                                                                 ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                                       >  1:588038/1‑141 (MQ=32)
                                                                                 ccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACtat                                                       >  1:377630/1‑141 (MQ=32)
                                                                                   ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAGGACTCCTTTTTGTTAACATTATACTCTaa                                                     >  1:74860/1‑141 (MQ=12)
                                                                                      atatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                                                  <  2:367354/141‑1 (MQ=35)
                                                                                      atatTTCAGATATAACACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAtt                                                  <  2:120770/141‑1 (MQ=25)
                                                                                       aatTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAGTATACTATAAATTc                                                 <  2:641276/140‑1 (MQ=25)
                                                                                           actgATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGTAGACTCCTGTTTGTTAAAATTATACTATAAATTCAACt                                             <  1:187929/138‑1 (MQ=255)
                                                                                              gATATATCACGATAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAATTTCAACTTTg                                          <  1:683885/141‑1 (MQ=255)
                                                                                               atatatCACGATAGCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                                         <  1:587874/141‑1 (MQ=255)
                                                                                               atatatCACGATAGCTCAATGGATCTCATGGAGAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGc                                         >  1:456209/1‑141 (MQ=11)
                                                                                                 atatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATATCATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCaa                                       <  1:538986/141‑1 (MQ=25)
                                                                                                  tatCACGATAGCTCAATGTATATATTAGTTAGTAGCCAACGGCTACGGTGATAACATCCTTGTTAAATTGTTGATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCGACTTTGCAAc                                      <  2:184198/141‑1 (MQ=11)
                                                                                                  tatCACGATAACTCAATGTAAATTGTAGCTAGTATCCAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTATACCTGAAATATTTCATACAGAAGACTCCTTTTTGTTAAAATTCTACTATACATTCAACTTTGCaca                                      >  2:435309/1‑139 (MQ=11)
                                                                                                     cACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCATTTTTGTTAAAATTATACTATAAATTCAACTTTCCacccga                                   >  1:116292/1‑136 (MQ=255)
                                                                                                          tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                              <  1:868682/138‑1 (MQ=9)
                                                                                                          tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                              <  2:727138/138‑1 (MQ=12)
                                                                                                          aagcctaaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                              <  2:324286/135‑1 (MQ=255)
                                                                                                           agcTCAATTCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGta                             <  2:598123/139‑1 (MQ=255)
                                                                                                            gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                            <  2:770435/140‑1 (MQ=16)
                                                                                                            aCTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat                            <  1:817193/141‑1 (MQ=38)
                                                                                                             cTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatt                           >  2:578432/1‑141 (MQ=39)
                                                                                                             cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                           >  1:788610/1‑141 (MQ=12)
                                                                                                             aatcaacttataCCTCAGATAGTAGCCAACGCCTACAGCTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTCCACACCGAAGACCCCTGTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat                            <  1:431876/133‑1 (MQ=11)
                                                                                                              tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta                          <  2:502065/141‑1 (MQ=255)
                                                                                                              tCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtatta                          >  1:378210/1‑141 (MQ=40)
                                                                                                              tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTTTTTTATTAATTCAACTTTGCAACAGAACCGttttt                          >  2:611834/1‑140 (MQ=11)
                                                                                                               cAATGTATAGCGTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtattat                         <  1:823997/141‑1 (MQ=255)
                                                                                                                aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATACAACTTTGCAACAGAACCGtatttta                        >  2:583813/1‑138 (MQ=255)
                                                                                                                 aattgatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATg                        <  2:589180/136‑1 (MQ=255)
                                                                                                                 aaaccatatCTTAAATATTAGCCAACGACTACAGTGATAACATCCTTGTTAAATTGTTTATATATGAAATAGTTCATACAGAATACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaact                        <  2:857539/136‑9 (MQ=11)
                                                                                                                  atgcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAACGtattttat                       >  1:372918/5‑137 (MQ=255)
                                                                                                                   tgcatatCTTATATATTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatg                      <  2:230690/138‑6 (MQ=11)
                                                                                                                     catatatTAGATAGTAGACAACGGCTACAGTGATAACATCCTCGTTAAATTGTTTACCTCTGAAATAGTTCATACAGACGACTATTTTTTGTTAAAATTATATTATAAATGCAACTTTGCATCAGAACCGtattttatgag                    <  2:818040/140‑8 (MQ=11)
                                                                                                                     catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                    <  1:483229/140‑8 (MQ=255)
                                                                                                                     catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgag                    <  1:745899/140‑8 (MQ=255)
                                                                                                                      atatCTTAGATAGTAGCCAACGGCTACAGTTATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATa                   <  1:124551/141‑1 (MQ=255)
                                                                                                                      atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTTAAATTATATTATAAATTCAACTTTGCAACAGAAcatcctaactgtatt                   >  2:586003/1‑127 (MQ=11)
                                                                                                                      atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatt                   <  2:309778/141‑14 (MQ=11)
                                                                                                                       tatCTTAGCTGTTATCCAACAGCTACAGTGATTATTCCCGTGTTCAATTGTTTCTTGCACAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  >  2:429497/1‑127 (MQ=2)
                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATTTGAAATAGTGCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                  <  1:220023/141‑10 (MQ=11)
                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAAcctcctaactgtatta                  >  2:828973/1‑127 (MQ=11)
                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatta                  >  2:162247/1‑127 (MQ=11)
                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                  <  1:235621/141‑10 (MQ=255)
                                                                                                                       tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAAAAGAAccacctaactttatta                  >  2:496240/1‑127 (MQ=11)
                                                                                                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  >  1:298913/1‑127 (MQ=11)
                                                                                                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  <  1:774077/141‑15 (MQ=11)
                                                                                                                       tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                  >  1:299376/1‑127 (MQ=11)
                                                                                                                         tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                >  1:929303/1‑130 (MQ=11)
                                                                                                                         tCTTACATAGTTGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACATAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                <  2:771740/141‑12 (MQ=11)
                                                                                                                           ttAGATAGTAGCCAACGGCTAGAGTGATAACATCCTTGTTAATTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact              <  2:4924/141‑14 (MQ=11)
                                                                                                                           ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact              >  1:300190/1‑128 (MQ=11)
                                                                                                                           ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact              >  1:792446/1‑128 (MQ=11)
                                                                                                                            tAGATAGTAGCCAACGGCTAGAGTGATAACATCCTTGTTAAATGGTTTATATCTGAAATAGTTCATACAGAAGGCTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc             <  2:343271/141‑15 (MQ=11)
                                                                                                                            tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc             <  2:924246/141‑15 (MQ=11)
                                                                                                                            tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACTGTATTATGGAATAGAGATg             <  1:519723/141‑1 (MQ=255)
                                                                                                                            tAGATAGTAGCCAACGGCGACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGCACCGtattttatgagtacacgc             <  1:132767/141‑15 (MQ=11)
                                                                                                                              gatatgcgCCAACGGTCACAGTGATAACATCCTGGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGCCTCCTTTTTGTTAACCTTATACACCAAATCCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGtt           <  2:387371/134‑1 (MQ=17)
                                                                                                                                   tagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTaa      >  2:761238/1‑141 (MQ=255)
                                                                                                                                   tagCCAACCGCTACAGTTATAACATCTTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTaa      <  1:417290/141‑1 (MQ=255)
                                                                                                                                    agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc     >  2:753232/1‑141 (MQ=255)
                                                                                                                                    agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc     >  1:331469/1‑141 (MQ=255)
                                                                                                                                           |                                                                                                                                        
CTTCCAAATTTGATACAAAATTGGGGCATATTCTTGAACCCAACGGTAGACCGTTGAATGATGAACGTTTACACCACGTTCCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATT  >  CP000730/36363‑36637

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: