Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,950,616 A→G 7.7% L320P (CTG→CCG) ‡ flhB ← flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,950,6160AG 7.7% 912.5 / 43.3 299L320P (CTG→CCG) ‡flhBflagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base A (135/140);  new base G (14/9);  total (150/149)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 304 total aligned reads displayed.

GGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGG  >  NC_002947/4950467‑4950765
                                                                                                                                                     |                                                                                                                                                     
ggTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCa                                                                                                                                                       <  2:2902810/150‑1 (MQ=255)
ggTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCa                                                                                                                                                       <  1:6354999/150‑1 (MQ=255)
ggTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCa                                                                                                                                                       <  2:5801931/150‑1 (MQ=255)
  tcgtcCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAgc                                                                                                                                                     >  2:3499897/1‑150 (MQ=255)
  tcgtcCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAgc                                                                                                                                                     >  1:6451317/1‑150 (MQ=255)
     tcCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGGGGAGTACTAGATCGCCCGCGCCAGCGcc                                                                                                                                                  >  2:5094292/1‑150 (MQ=255)
       cTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCggg                                                                                                                                                <  2:2785971/150‑1 (MQ=255)
       cTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCgg                                                                                                                                                <  1:6149355/150‑1 (MQ=255)
        ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGc                                                                                                                                               >  1:1117359/1‑150 (MQ=255)
        ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGc                                                                                                                                               >  1:5974626/1‑150 (MQ=255)
        ttCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGc                                                                                                                                               >  1:4932610/1‑150 (MQ=255)
           cgTGGCTCTGGGCGCGTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGAt                                                                                                                                            <  2:5570504/149‑1 (MQ=255)
             tGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTc                                                                                                                                          <  2:5188424/150‑1 (MQ=255)
             tGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTc                                                                                                                                          <  1:4856464/150‑1 (MQ=255)
              ggCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGTTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtcta                                                                                                                                         <  1:5662785/150‑4 (MQ=255)
               gCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                        <  2:717764/150‑1 (MQ=255)
               gCTCTGGGCGCGTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTAGAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGa                                                                                                                                        <  1:2730933/150‑1 (MQ=255)
                ctctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCTAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAg                                                                                                                                       <  1:6243969/150‑1 (MQ=255)
                ctctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAg                                                                                                                                       <  2:2459358/150‑1 (MQ=255)
                ctctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAg                                                                                                                                       <  1:967308/150‑1 (MQ=255)
                 tctGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCGGCGCTGGCGCGGGCGAtctatta                                                                                                                                      <  1:2156848/150‑7 (MQ=255)
                  ctGGGCGCGTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATCGATCGCCCGCGCCAGCGCCGGCGATTCGAGGa                                                                                                                                     <  1:3806903/150‑1 (MQ=255)
                   tGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGAt                                                                                                                                    <  1:2246501/150‑1 (MQ=255)
                    gggCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATc                                                                                                                                   >  1:553218/1‑150 (MQ=255)
                        gcTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGa                                                                                                                               <  1:1714537/150‑1 (MQ=255)
                            gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                           <  2:2581449/150‑1 (MQ=255)
                            gcccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgtt                                                                                                                           <  2:6602401/150‑1 (MQ=255)
                              ccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgttgt                                                                                                                         >  2:4553424/1‑150 (MQ=255)
                              ccttgccGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGAtgttgt                                                                                                                         >  2:1709496/1‑150 (MQ=255)
                                 ggccGGCGCGGTATGGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCt                                                                                                                      <  2:5003762/149‑1 (MQ=255)
                                     ggCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCAc                                                                                                                  >  2:6777759/1‑150 (MQ=255)
                                          ggTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCaa                                                                                                             <  2:4080026/150‑1 (MQ=255)
                                           gTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAt                                                                                                            <  2:6286910/150‑1 (MQ=255)
                                            tATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCCGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAAtt                                                                                                           >  2:6026153/1‑150 (MQ=255)
                                               tGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCa                                                                                                        <  1:4407432/150‑1 (MQ=255)
                                               tGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCa                                                                                                        <  2:3586574/150‑1 (MQ=255)
                                                ggCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCAc                                                                                                       <  2:1212233/150‑1 (MQ=255)
                                                 gCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGAGCTCCACGCCAATATCACg                                                                                                      >  2:6207670/1‑150 (MQ=255)
                                                     aTCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATc                                                                                                  <  1:2521494/150‑1 (MQ=255)
                                                         ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                              >  2:32/1‑150 (MQ=255)
                                                         ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                              >  1:4328326/1‑150 (MQ=255)
                                                         ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                              >  2:2228920/1‑150 (MQ=255)
                                                         ggAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCa                                                                                              >  1:4011607/1‑150 (MQ=255)
                                                          gAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCaa                                                                                             <  2:3050456/150‑1 (MQ=255)
                                                           aaCACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAAc                                                                                            <  1:1255961/150‑1 (MQ=255)
                                                            acacATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATATCACGAATCTTCAACg                                                                                           >  2:1434361/1‑150 (MQ=255)
                                                             cacaGACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGc                                                                                          <  1:4037965/150‑1 (MQ=255)
                                                              acaTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGcc                                                                                         >  1:5571755/1‑150 (MQ=255)
                                                                aTACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCAt                                                                                       >  1:5332802/1‑150 (MQ=255)
                                                                 tACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATGGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATg                                                                                      <  2:3677934/150‑1 (MQ=255)
                                                                  aCGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTACTAGATCGCCCGCGCCAGCGCCGGCGCTTCGAGGAGCTGGATGGGGTGCTCCACGCCAATTTCACGACTCTTCAACGCCCTGa                                                                                     >  2:4600109/1‑150 (MQ=255)
                                                                    gCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGGTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGaac                                                                                   >  2:5726373/1‑149 (MQ=255)
                                                                     ccAGCACCTGAGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAATTTCGGTGGAGTAATAGATCGACCGCGCCAGCGCCGGCGATACGCGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAAt                                                                                  >  2:345334/1‑150 (MQ=255)
                                                                      cAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATc                                                                                 <  1:3852375/150‑1 (MQ=255)
                                                                       aGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCg                                                                                <  1:3003371/150‑1 (MQ=255)
                                                                        gCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATCTCACGAATCTTCAACGCCATGAAATCGc                                                                               <  1:6117797/150‑1 (MQ=255)
                                                                          aCCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTa                                                                             <  1:2037227/150‑1 (MQ=255)
                                                                                gcCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTg                                                                       <  1:3302270/150‑1 (MQ=255)
                                                                                 ccACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTgg                                                                      >  1:2833667/1‑150 (MQ=255)
                                                                                  cACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGc                                                                     >  1:2432712/1‑150 (MQ=255)
                                                                                   aCCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGcc                                                                    >  1:4964016/1‑150 (MQ=255)
                                                                                       ccAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCcagca                                                                <  2:5636559/150‑1 (MQ=255)
                                                                                            tACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGATCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAAc                                                           >  2:1880688/1‑150 (MQ=255)
                                                                                                      ccgGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCAc                                                 >  2:6652661/1‑150 (MQ=255)
                                                                                                         ggATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGc                                              >  1:5402942/1‑150 (MQ=255)
                                                                                                         ggATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGc                                              <  1:1816105/150‑1 (MQ=255)
                                                                                                          gATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGcc                                             >  1:5856568/1‑150 (MQ=255)
                                                                                                               ccTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTtct                                        >  1:891666/1‑150 (MQ=255)
                                                                                                                 tGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTcc                                      <  2:4375349/150‑1 (MQ=255)
                                                                                                                   cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATc                                                                                                  >  1:3597999/1‑88 (MQ=255)
                                                                                                                   cTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCgg                                    <  1:4528450/150‑1 (MQ=255)
                                                                                                                    tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                   >  2:3947285/1‑150 (MQ=255)
                                                                                                                    tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                   <  1:4695554/150‑1 (MQ=255)
                                                                                                                    tCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCggg                                   >  2:2559103/1‑150 (MQ=255)
                                                                                                                         gtCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTa                              <  2:3850719/149‑1 (MQ=255)
                                                                                                                           cGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTAtt                            <  2:4739076/150‑1 (MQ=255)
                                                                                                                             gTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGc                          >  1:1600808/1‑150 (MQ=255)
                                                                                                                             gTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGc                          <  2:6451317/150‑1 (MQ=255)
                                                                                                                                gAGTAATAGATCGCCCGCGCCAGCGCCGGCGTTTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                       <  1:34998/150‑1 (MQ=255)
                                                                                                                                gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                       <  2:931613/150‑1 (MQ=255)
                                                                                                                                gAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCaaa                       <  2:1790338/150‑1 (MQ=255)
                                                                                                                                   tAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGcc                    >  1:3448524/1‑150 (MQ=255)
                                                                                                                                   cctcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGcc                    >  1:1468607/7‑150 (MQ=255)
                                                                                                                                    ctcgaaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                   >  2:3351478/6‑150 (MQ=255)
                                                                                                                                    aaTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAACTCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCa                   <  1:5726373/150‑1 (MQ=255)
                                                                                                                                     aTAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCAc                  >  1:6332836/1‑150 (MQ=255)
                                                                                                                                        gATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGc               >  1:6837702/1‑150 (MQ=255)
                                                                                                                                        aaTCGCCGGCGCTGGCGCGGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGc               >  1:516102/2‑150 (MQ=255)
                                                                                                                                         aTCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTc                                                                                                         <  1:2142455/59‑1 (MQ=255)
                                                                                                                                           cGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATa            <  2:3981923/150‑1 (MQ=255)
                                                                                                                                            gCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATaa           <  1:1405558/150‑1 (MQ=255)
                                                                                                                                              ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATg         >  2:5507847/1‑150 (MQ=255)
                                                                                                                                              ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATg         >  1:2599583/1‑150 (MQ=255)
                                                                                                                                              ccGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATg         >  1:1044948/1‑150 (MQ=255)
                                                                                                                                                gcgcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCg       >  2:2070886/1‑150 (MQ=255)
                                                                                                                                                gcgcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCg       >  1:5055499/1‑150 (MQ=255)
                                                                                                                                                  gcCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTc     >  2:4038704/1‑150 (MQ=255)
                                                                                                                                                   ccAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCg    <  2:2952906/150‑1 (MQ=255)
                                                                                                                                                    cAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCgg   <  2:4932610/150‑1 (MQ=255)
                                                                                                                                                    cAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCgg   >  1:1530005/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTCGTCCTTCAGTGGCTCTGGGCGCTTGCCCTTGCCGGCGCGGTATTGGCGAATCTGGAACACATACGCCAGCACCTGCGCCACCGCCAGGTACAGCCCCGCCGGGATCTCCTGCTCAAGTTCGGTGGAGTAATAGATCGCCCGCGCCAGCGCCGGCGATTCGAGGATCTGGATGTTGTGCTCCACGCCAATTTCACGAATCTTCAACGCCATGAAATCGCTACCCTTGGCCAGCAGCAACGGCGCCGCACCGCCTTTCTCCGGGTCGTATTGCAAAGCCACCGCATAATGCGTCGGG  >  NC_002947/4950467‑4950765

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