Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP000730 36,581 C→T 43.6% intergenic (‑30/‑493) USA300HOU_0026 ← / → USA300HOU_0028 IS431mec transposase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00073036,5810CT43.6% 76.4 / 99.5 91intergenic (‑30/‑493)USA300HOU_0026/USA300HOU_0028IS431mec transposase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (25/25);  new base T (16/23);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.62e-01

CCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTT  >  CP000730/36443‑36718
                                                                                                                                          |                                                                                                                                         
cccTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACta                                                                                                                                         >  2:26746/1‑141 (MQ=32)
  cTTAATATTTCAGATATATCACGCTAGCTCAATGTCTATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAAttatattata                                                                                                                                       <  2:345515/141‑1 (MQ=255)
   ttAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATaa                                                                                                                                      >  1:303476/1‑141 (MQ=33)
     aaTATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAAt                                                                                                                                    <  1:113173/141‑1 (MQ=35)
        atgtCAGATATATCACGATAACTCAATGCATATCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                                                                                                                                 <  1:141912/138‑1 (MQ=2)
        atTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAACTAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCa                                                                                                                                 <  1:201775/141‑1 (MQ=21)
          ttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAAc                                                                                                                               >  1:223302/1‑141 (MQ=35)
          ttCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTAAAc                                                                                                                               >  2:139200/1‑141 (MQ=25)
           tCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACt                                                                                                                              <  2:69889/141‑1 (MQ=35)
              gATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTCTACTATAAATTCAACTTTg                                                                                                                           >  2:345023/1‑141 (MQ=25)
              gATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTg                                                                                                                           >  1:342749/1‑141 (MQ=35)
                 atatCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGccg                                                                                                                        <  2:99578/141‑3 (MQ=25)
                     cACGATGACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGa                                                                                                                    >  2:140014/1‑141 (MQ=21)
                         atagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg                                                                                                                <  1:176685/137‑1 (MQ=0)
                         aTAGCTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCg                                                                                                                <  2:65519/141‑1 (MQ=2)
                          tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                               >  1:127458/4‑141 (MQ=9)
                          tagcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGt                                                                                                               <  1:169282/138‑1 (MQ=12)
                           agcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGGTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCACCAGAACCGta                                                                                                              >  1:2160/3‑141 (MQ=11)
                            gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAGCAGAACCGtat                                                                                                             >  1:305970/2‑141 (MQ=12)
                            gcTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtat                                                                                                             >  1:61650/2‑141 (MQ=12)
                            aCTCAATGAATCTCTTAGACAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCTTACAGAAGACTCCTTTCTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGtat                                                                                                             <  2:192364/141‑1 (MQ=11)
                             cTCAATTCACATCTTAGATAGTATCCAACGGCCACATTCATAACATCTTTTGTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTTCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                            <  2:268448/141‑1 (MQ=11)
                             cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATGCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAAACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAGCAGAACCGtatt                                                                                                            >  2:29853/1‑141 (MQ=11)
                             cTCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatt                                                                                                            >  1:87531/1‑141 (MQ=16)
                              tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGCTAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGGCTCCTTTTTGTTAAAATTATATTATAAATTCAATTTTGCAACAGAACCGtattt                                                                                                           <  1:248228/141‑2 (MQ=11)
                              tCAATGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattt                                                                                                           >  2:76585/1‑140 (MQ=12)
                                aaTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtatttta                                                                                                         <  1:122342/141‑4 (MQ=255)
                                 aTGCATATCTTCGTTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat                                                                                                        <  2:125069/141‑5 (MQ=11)
                                 aTGCATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttat                                                                                                        <  1:17275/141‑5 (MQ=255)
                                  agcatacCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAAcctcatttacta                                                                                                       <  2:286067/138‑11 (MQ=11)
                                   gcatttCTTAGATAGTAGCCAACGGGTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                                                                                                      <  2:77520/136‑7 (MQ=11)
                                   gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTTCAACAGAAcctcctaactgta                                                                                                      <  1:198179/139‑12 (MQ=11)
                                   gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                                                                                                      <  1:340636/139‑7 (MQ=255)
                                   gcatatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAAGTGTTTATATCTGAAATAGATCATACATGAGACTCCTTTTTGTTAAAACTATATGATGAATTCAACGTTGCAACTGAAcctccaaacttga                                                                                                      <  2:157095/139‑12 (MQ=1)
                                   gcatatCTTAGATAGTAGCCAACGGCTACAGCGATAACATCCTTGTAAAATTTTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatga                                                                                                      <  1:314581/139‑7 (MQ=11)
                                    tataTCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAAt                                                                                                     <  1:287619/141‑1 (MQ=255)
                                    catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                                                                                                     >  1:152082/2‑129 (MQ=11)
                                    catatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtat                                                                                                     <  1:133837/140‑13 (MQ=11)
                                     atatCTTAGATAGTATCCAAAGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagt                                                                                                    <  2:80092/141‑9 (MQ=11)
                                     atatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctgctaactgtatt                                                                                                    >  2:59097/1‑128 (MQ=11)
                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtatta                                                                                                   >  2:342005/1‑127 (MQ=11)
                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactgtattc                                                                                                   >  1:223711/1‑127 (MQ=11)
                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                                                                                                   >  1:97026/1‑132 (MQ=255)
                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagta                                                                                                   >  1:359421/1‑132 (MQ=255)
                                      tatCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACACAAcctcctaactgtattc                                                                                                   >  2:223297/1‑127 (MQ=11)
                                      tatCTTAGATAGTAGCCAACAGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAATATTATACTATAAATTCAACTTTGCAACAGAAcctcatttactatatg                                                                                                   >  1:135690/1‑127 (MQ=11)
                                       atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATCATAAATTCAACTTTGCAACAGAACCGtattttatgagtac                                                                                                  >  1:327800/1‑131 (MQ=11)
                                       atCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATACTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtac                                                                                                  <  2:277344/141‑11 (MQ=11)
                                        tCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATCGTTGAAAGAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATagag                                                                                                 <  1:181462/141‑1 (MQ=255)
                                        tCTTAGATACTAGCCAACGGCTACATTTATAACATCCTTGTTAAATTTTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtaca                                                                                                 <  1:76690/141‑12 (MQ=11)
                                          ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                                                                                               <  1:340574/141‑14 (MQ=11)
                                          ttAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacact                                                                                               <  1:382091/141‑14 (MQ=11)
                                           tAGGTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                              <  2:122687/141‑15 (MQ=11)
                                           tAGCTAGTTGCCAACGGCTACAGTGTGAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                              <  2:228648/141‑15 (MQ=11)
                                           tAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAACCGtattttatgagtacactc                                                                                              <  2:279051/141‑15 (MQ=11)
                                            aGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGt                                                                                             <  1:26746/141‑1 (MQ=255)
                                             gATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGtt                                                                                            <  2:223302/141‑1 (MQ=255)
                                              aTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactg                                                                                                        >  2:99566/1‑119 (MQ=11)
                                              aTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATATTATAAATTCAACTTTGCAACAGAAcctcctaactg                                                                                                        <  1:99566/128‑10 (MQ=11)
                                              aTAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTg                                                                                           <  1:58922/141‑1 (MQ=255)
                                                agtagCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGt                                                                                         <  2:144801/141‑1 (MQ=255)
                                                   agCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAAc                                                                                      >  2:256881/1‑141 (MQ=255)
                                                      cAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTTTACTGTAAATTCAACTTTGCAACAGAACCGTCTTATGGCATAGAGATGTGGGGAACAtt                                                                                   >  2:105011/1‑141 (MQ=14)
                                                             tACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAg                                                                            >  2:185933/1‑141 (MQ=255)
                                                              aCAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTTTACAgg                                                                           >  1:248339/1‑141 (MQ=255)
                                                                    aTAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCAtt                                                                     <  2:217070/141‑1 (MQ=255)
                                                                             ttGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAg                                                            >  2:276734/1‑141 (MQ=255)
                                                                             ttGTTAAATTGTTTATATCTGAAAGAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATAATATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAg                                                            <  1:30958/141‑1 (MQ=255)
                                                                               gTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGtt                                                          >  1:155406/1‑141 (MQ=255)
                                                                                    aTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAAtt                                                     >  1:238212/1‑141 (MQ=255)
                                                                                    aTTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAAtt                                                     >  1:53474/1‑141 (MQ=255)
                                                                                        tttATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGt                                                 <  2:186337/141‑1 (MQ=255)
                                                                                            tatCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTAGTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGttatt                                             <  2:221414/141‑1 (MQ=255)
                                                                                                  aaaTAGTTCATACATAAGCCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGaa                                       <  1:129719/141‑1 (MQ=255)
                                                                                                    aTAGTTCATACAGAAGACTCCTTTTTTTTTTTTTTATACAATAAAATCAACTTTACAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTCTACTAAAGTTTAATTTCGTTATTACATAAcc                                     >  2:223390/1‑141 (MQ=37)
                                                                                                    aTAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATGATTATGGTTAAGTTTACTTTCGTTGTTACGGAAcc                                     >  1:151972/1‑141 (MQ=25)
                                                                                                    aTAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTTCTTATTACAGaatc                                     >  2:244500/1‑139 (MQ=255)
                                                                                                        ttCATACAGAAGGCTCCGTTTTGTTAAAATTATACGAGAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATGATACTTAAGTTTAATTTCGTTATTACAGAACcacac                                 <  1:321464/141‑1 (MQ=38)
                                                                                                        ttCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACcacac                                 <  2:199529/141‑1 (MQ=255)
                                                                                                         tCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACATGATCATTATACTTAAGATTAATTTCGTTATTGACGAACcacaca                                >  2:214588/1‑141 (MQ=37)
                                                                                                                    aCTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACAACACATTACAACCaga                     >  2:344959/1‑141 (MQ=255)
                                                                                                                        cttttttttttttttATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAagag                 >  2:377028/14‑141 (MQ=255)
                                                                                                                        ctgtttgggaaaaGTATACTATAAATTCAACTTTGCAACAGATCCGTATTATGGAATAGATATGTTGGTAACATTTATACATGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAagag                 <  2:342102/132‑1 (MQ=255)
                                                                                                                        cTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAagag                 <  2:99461/141‑1 (MQ=255)
                                                                                                                         tttttGTTAAAATTATACTATAAATTCAACTTTGCAACCGAACCGGATTATGGAATAGAGATGTTGGTAACCTTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGCAgaga                <  2:342749/141‑1 (MQ=255)
                                                                                                                         tttttGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAgaga                >  1:302831/1‑141 (MQ=255)
                                                                                                                         tttttGTTAAAATTAGACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACATGATCATTATACTTAAGTTTCGTTTCGTTATTACAGAACCACACATTCCAACCAGAAgaga                <  2:251273/141‑1 (MQ=255)
                                                                                                                          ttctGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTGTACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGaa               <  1:256881/138‑1 (MQ=255)
                                                                                                                               ttAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTAt          <  2:303476/141‑1 (MQ=255)
                                                                                                                                 aaaaTTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGt        >  2:222641/1‑141 (MQ=255)
                                                                                                                                    aTTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTa     >  2:240880/1‑141 (MQ=255)
                                                                                                                                       ataGTATAAATTCGACTTCGCAACAGAACCGTATTATGGAATAGGGAGGTTGGTCAAATGTATACAGGATCATTATACACAAGTTTAATTTCGGTGTTACAGAAACCCACATTCCAACGACAACAGTAAGTATGtgtgcga  >  2:95166/1‑135 (MQ=255)
                                                                                                                                          |                                                                                                                                         
CCCTTAATATTTCAGATATATCACGATAACTCAATGTATATCTTAGATAGTAGCCAACGGCTACAGTGATAACATCCTTGTTAAATTGTTTATATCTGAAATAGTTCATACAGAAGACTCCTTTTTGTTAAAATTATACTATAAATTCAACTTTGCAACAGAACCGTATTATGGAATAGAGATGTTGGTAACATTTATACAGGATCATTATACTTAAGTTTAATTTCGTTATTACAGAACCACACATTCCAACCAGAAGAGAAAGTATGTCTATTT  >  CP000730/36443‑36718

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: