Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 76,295 T→C 53.3% intergenic (+353/‑781) USA300HOU_0072 → / → USA300HOU_0075 universal stress protein/lysophospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073076,2950TC53.3% 64.2 / inf 428intergenic (+353/‑781)USA300HOU_0072/USA300HOU_0075universal stress protein/lysophospholipase
Reads supporting (aligned to +/- strand):  ref base T (103/95);  new base C (115/110);  total (221/207)
Fisher's exact test for biased strand distribution p-value = 9.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 429 total aligned reads displayed.

ATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCTAATTACAAATGGCG  >  CP000730/76156‑76434
                                                                                                                                           |                                                                                                                                            
ataGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGtt                                                                                                                                             >  2:267952/1‑141 (MQ=38)
 taGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTa                                                                                                                                            >  1:608031/1‑141 (MQ=38)
  aGGTTCTGTTGCAAAGTTGAATTTATAATATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTAt                                                                                                                                           >  2:35738/1‑141 (MQ=0)
   ggTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATc                                                                                                                                          <  1:92884/141‑1 (MQ=35)
            gCAAAGTTGAATCTATAGGATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGatatat                                                                                                                                 >  2:61211/1‑141 (MQ=21)
               aaGTTGAATCTATAGCATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTg                                                                                                                              >  2:574360/1‑141 (MQ=25)
                       tttataatataaTTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATGTCTGAAATATTa                                                                                                                      <  1:472470/139‑1 (MQ=9)
                        ttataGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTaa                                                                                                                     <  1:659973/140‑1 (MQ=17)
                           taatataATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAgag                                                                                                                  <  2:323665/138‑3 (MQ=12)
                           taGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAggg                                                                                                                  >  1:256160/1‑141 (MQ=17)
                            aatataATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAgaga                                                                                                                 <  2:621602/139‑1 (MQ=12)
                               agatttttAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACg                                                                                                              <  1:33088/137‑1 (MQ=12)
                                  aTTTTAACAAAAAGGAGTCTTCTGTATGATGTATACCTGTTTTTTACAATTTAGCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAAttttctgagcaagtgg                                                                                                           <  2:302853/141‑13 (MQ=11)
                                    tttAACAAAAAGGATTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACgtggtg                                                                                                         <  2:336640/141‑1 (MQ=12)
                                       aaCAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACTATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTaa                                                                                                      <  1:417074/141‑1 (MQ=255)
                                       aaCAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTaa                                                                                                      >  1:72819/1‑141 (MQ=32)
                                         cAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTaaag                                                                                                    >  2:36917/1‑140 (MQ=17)
                                          aaaaaGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACg                                                                                                   <  2:159617/141‑1 (MQ=17)
                                          aaaaaGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACg                                                                                                   <  1:399358/141‑1 (MQ=12)
                                             aaGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTtc                                                                                                >  2:40549/1‑141 (MQ=12)
                                             aaGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTtc                                                                                                >  2:382048/1‑141 (MQ=17)
                                              aGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTtca                                                                                               <  1:59170/141‑1 (MQ=17)
                                                  gTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATGTATCTGAAATATTAAGGGGAGGTGGTGTAAAGGTTcatcat                                                                                           >  1:246325/1‑141 (MQ=11)
                                                       ctGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAAc                                                                                      >  2:414245/1‑141 (MQ=17)
                                                        tGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACg                                                                                     <  2:210525/141‑1 (MQ=17)
                                                           aTGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACgatc                                                                                  <  1:626411/141‑4 (MQ=255)
                                                            tGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGtag                                                                                 <  2:568279/141‑3 (MQ=17)
                                                               aCTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTAcc                                                                              >  2:516439/1‑141 (MQ=12)
                                                                 tATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGt                                                                            >  1:378813/1‑141 (MQ=12)
                                                                  aTTTCCGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGtt                                                                           <  1:262413/141‑1 (MQ=255)
                                                                  aTTTCAGATGTAAACAATGTGACAAGGGTGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGGGGTGTAAACGTTCATCATTCAACGATCTACCGtt                                                                           <  2:229293/141‑1 (MQ=11)
                                                                  aTTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGtt                                                                           <  2:666831/141‑1 (MQ=32)
                                                                  aTTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGACCGTGGTGTAAACGTTCATCATTCAACGATCTACCGtt                                                                           <  1:191533/141‑1 (MQ=255)
                                                                   tttCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTg                                                                          >  1:152817/1‑141 (MQ=12)
                                                                    ttCATATAGAATCTATTTAACAATGCTTTTACCACTGTAGCTGCTCGCTACTATCTAAGATATGCACTGCGTTATCTTGTTATCACTCAAATATTAAGGGCACTTGGTGTAATCGTTCATCATTCAACGGTAGACCGTTga                                                                         <  2:217469/141‑2 (MQ=255)
                                                                    ttCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTgg                                                                         <  1:145388/141‑1 (MQ=255)
                                                                     tCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTggg                                                                        <  1:447821/141‑1 (MQ=17)
                                                                     tCAGATATAAACAATTTAACAAGGATGTTATCACTGTAACTGTTGGCTACTATCTAAGATATGCAGTGAGTTATCATGATATATCTGAAATATTAAGGAAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGGTggg                                                                        >  2:258979/1‑141 (MQ=11)
                                                                      cAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTGAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGACCGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGt                                                                       >  2:578743/1‑141 (MQ=17)
                                                                      cAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGt                                                                       >  1:113826/1‑141 (MQ=12)
                                                                      cAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGt                                                                       >  1:230715/1‑141 (MQ=255)
                                                                       aGATATAAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGtt                                                                      >  1:286891/1‑141 (MQ=17)
                                                                       aGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGGACGTGGTGTAAACGTTCAACATTCAACGATCTGACGTTGGGtt                                                                      >  2:353202/1‑141 (MQ=11)
                                                                        aatagAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTc                                                                     <  2:50318/140‑1 (MQ=11)
                                                                         atatAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCa                                                                    <  2:521821/141‑1 (MQ=255)
                                                                         atatAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCa                                                                    <  2:43411/141‑1 (MQ=12)
                                                                          tataAACAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCaa                                                                   >  1:526134/1‑141 (MQ=17)
                                                                           ttaAACAATTTAACACGGATGGTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAg                                                                  <  2:455354/140‑1 (MQ=12)
                                                                            taAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATCCAACGATCTACCGTTGGGCTCAAGa                                                                 >  2:74244/1‑141 (MQ=11)
                                                                              aaCAATTTAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGGGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAat                                                               <  1:628098/141‑1 (MQ=17)
                                                                              aaCAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAACCGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAat                                                               >  1:646573/1‑141 (MQ=255)
                                                                                aaaTTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAGCGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAatat                                                             <  2:314782/140‑1 (MQ=255)
                                                                                 aaTTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATg                                                            <  2:653543/141‑1 (MQ=12)
                                                                                  aTTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGc                                                           <  1:27230/141‑1 (MQ=32)
                                                                                  aTTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCCTCATTCACCGATCTACCGCTGGGTTCAAGAATAAGc                                                           >  2:189269/1‑141 (MQ=11)
                                                                                   tttAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGcc                                                          <  1:504643/141‑1 (MQ=32)
                                                                                   tttAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGcc                                                          <  2:543639/141‑1 (MQ=32)
                                                                                     tAACAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTGATCATTCAACGGCCTACCGTTGGGTTCAAGAATATGcccc                                                        >  2:603962/1‑141 (MQ=12)
                                                                                     tAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGcccc                                                        >  1:601193/1‑141 (MQ=255)
                                                                                       aCAAGGATGTTATCACTGTAGCTGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTAGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCaa                                                      <  2:288805/141‑1 (MQ=17)
                                                                                       aCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCaa                                                      >  2:127290/1‑141 (MQ=12)
                                                                                       aCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCaa                                                      >  2:182529/1‑141 (MQ=12)
                                                                                       aCAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCACTCAACGAGCTACCGTTGGGTTCAAGAATATGCCCCaa                                                      >  1:122512/1‑141 (MQ=11)
                                                                                        cAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGCTATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAAt                                                     >  2:13621/1‑141 (MQ=255)
                                                                                        cAAGGATGTTATCACTGTAGATGTTGGCTACTATCTAAGATATGCATTGAGTTAACGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAAt                                                     >  1:433642/1‑141 (MQ=12)
                                                                                         aaGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCACAAtt                                                    >  2:398570/1‑141 (MQ=255)
                                                                                          aGGATGTTACCACGGTACCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAAttt                                                   <  2:285113/141‑1 (MQ=11)
                                                                                           ggATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAAtttt                                                  >  2:235743/1‑141 (MQ=32)
                                                                                           ggATGTCATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTAGCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAAtttt                                                  <  2:275761/141‑1 (MQ=17)
                                                                                                ttATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTCTATc                                             >  1:411176/1‑141 (MQ=255)
                                                                                                  aTCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCGCCAATTTTATATCaa                                           >  1:568915/1‑141 (MQ=21)
                                                                                                  aTCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAGGAATATGCCCCAATTTTATATCaa                                           >  1:98691/1‑141 (MQ=255)
                                                                                                  aTCACTGTAGCCGTTGGCTACTATATAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCaa                                           >  1:488001/1‑141 (MQ=255)
                                                                                                   acACGGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGCTCTCGTCATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATACCaaa                                          <  2:534712/140‑1 (MQ=11)
                                                                                                         gagatgTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACTTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATAAAAATTGga                                    <  2:149062/136‑1 (MQ=255)
                                                                                                         gaGCCGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGga                                    <  1:562707/140‑1 (MQ=255)
                                                                                                            ccGTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAGTATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAAGATGCCCCAATTTTATATCAAATTTGgaaga                                 >  2:491595/1‑141 (MQ=11)
                                                                                                             tgTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGgaagaa                                <  1:135895/140‑1 (MQ=25)
                                                                                                              gTTGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaa                               >  1:163764/1‑141 (MQ=255)
                                                                                                               ttGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaaa                              >  1:553628/1‑141 (MQ=35)
                                                                                                                tGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTCCCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaaaa                             >  2:88528/1‑141 (MQ=25)
                                                                                                                tGGCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGaaaaa                             >  2:478424/1‑141 (MQ=255)
                                                                                                                  gCTACTATCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGc                           <  1:281332/141‑1 (MQ=255)
                                                                                                                   ctactaCTTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCa                          <  2:642830/141‑1 (MQ=11)
                                                                                                                   ctacgaTCTAAGGTATGGATTGAGTTATCGTGATAGATCTGAAAGATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGGCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCa                          <  2:108027/141‑1 (MQ=11)
                                                                                                                   cgactaTCTAAGCTATGCATTGAGTTAGCGTGATATATCTGAAATATGAAGGGAACGTGGTGTACACGTTCATCATGCAACGGTCTACCGTTGGGTTAAAGAATATGCCCCAATTTTATATAATATTTAGAAGAAAAAGCa                          <  2:255164/139‑1 (MQ=11)
                                                                                                                    tactaTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCAt                         <  1:216659/141‑1 (MQ=34)
                                                                                                                    tactaTCTAAGATACGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCAt                         <  1:356569/141‑1 (MQ=255)
                                                                                                                    gactaTCTAAGATATGAATTGAGCTATCGTTATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTCCATCATTCAACGATCTAACTTTGGGTCCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCAt                         <  1:46023/140‑1 (MQ=11)
                                                                                                                     actaTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATa                       >  1:178821/1‑141 (MQ=35)
                                                                                                                     actaTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATa                       <  1:632438/141‑1 (MQ=255)
                                                                                                                      ctaTCTAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaa                      >  2:598065/1‑141 (MQ=255)
                                                                                                                       taTCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaa                     <  1:530352/141‑1 (MQ=35)
                                                                                                                          tgaaGGTATGTATTGAGGTATCGTGATAGATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaa                  <  2:444428/139‑2 (MQ=11)
                                                                                                                           tAAGATATGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaag                 >  2:494333/1‑139 (MQ=255)
                                                                                                                            aaGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCAGAAAAAGCt                >  1:338915/1‑141 (MQ=38)
                                                                                                                             aGATATGCATTGAGCTATCGTGCCTCTGCGTAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaagct               >  1:661998/1‑137 (MQ=11)
                                                                                                                                tatGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAAGCTAAt             >  1:86243/1‑141 (MQ=18)
                                                                                                                                tatGCATTGAGCTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaagcttat            <  1:316892/141‑8 (MQ=1)
                                                                                                                                 atGCATTGAGTTATCGTGATATATCTGAAATATTAAGAGAACGTGGTGTAAACGTTCATCATTCAACGATCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATaaaaaagcttagt           <  2:80052/141‑9 (MQ=1)
                                                                                                                                           |                                                                                                                                            
ATAGGTTCTGTTGCAAAGTTGAATCTATAGTATAATTTTAACAAAAAGGAGTCTTCTGTATGAACTATTTCAGATATAAACAATTTAACAAGGATGTTATCACTGTAGCCGTTGGCTACTATCTAAGATATGCATTGAGTTATCGTGATATATCTGAAATATTAAGGGAACGTGGTGTAAACGTTCATCATTCAACGGTCTACCGTTGGGTTCAAGAATATGCCCCAATTTTATATCAAATTTGGAAGAAAAAGCATAAAAAGCTAATTACAAATGGCG  >  CP000730/76156‑76434

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: