Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,975,709 A→T 5.8% L244M (TTG→ATG)  PP_3503 ← sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,975,7090AT 5.8% 993.3 / 32.0 297L244M (TTG→ATG) PP_3503sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (146/133);  new base T (10/8);  total (156/141)
Fisher's exact test for biased strand distribution p-value = 8.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 307 total aligned reads displayed.

TTCGCGGAAACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAG  >  NC_002947/3975560‑3975857
                                                                                                                                                     |                                                                                                                                                    
ttCGCGGAAACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCaa                                                                                                                                                      <  2:8947319/150‑1 (MQ=255)
ttCGCGGAAACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCaa                                                                                                                                                      >  2:7272719/1‑150 (MQ=255)
ttCGCGGAAACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCaa                                                                                                                                                      >  1:1052405/1‑150 (MQ=255)
 tCGCGGAAACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAAc                                                                                                                                                     >  1:1723083/1‑150 (MQ=255)
            cccTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGcc                                                                                                                                          <  1:8169414/150‑1 (MQ=255)
            cccTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGcc                                                                                                                                          >  2:1342319/1‑150 (MQ=255)
               ttGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGAt                                                                                                                                       >  2:6555247/1‑150 (MQ=255)
                  cTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtc                                                                                                                                    >  1:8441491/1‑150 (MQ=255)
                   tCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtca                                                                                                                                   >  2:5264911/1‑150 (MQ=255)
                    cGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcat                                                                                                                                  <  2:5632323/150‑1 (MQ=255)
                     gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                 <  1:8164247/150‑1 (MQ=255)
                     gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                 >  1:3756091/1‑150 (MQ=255)
                     gATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATCTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatc                                                                                                                                 >  2:8427973/1‑150 (MQ=255)
                      aTGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTtcatca                                                                                                                                <  1:3978918/150‑1 (MQ=255)
                       tGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCaa                                                                                                                               >  1:7158074/1‑150 (MQ=255)
                       tGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCaa                                                                                                                               <  1:3453668/150‑1 (MQ=255)
                          cccGCTCAAGGTCCACATGGGTGGCCGCCAACACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTTGCCGATTTCATCAAGaa                                                                                                                            <  2:2308726/150‑1 (MQ=255)
                              cTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACaa                                                                                                                        <  1:251137/150‑1 (MQ=255)
                              cTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACaa                                                                                                                        <  2:7380305/150‑1 (MQ=255)
                              cTCAAGGTCCACAAGGGTGGCCGCCAGCACCCGTACATCCACCGCTATGCGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCACCAACAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTACATCAAGAAACaa                                                                                                                        <  2:7576663/150‑1 (MQ=255)
                               tCAAGGTCCTCATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAAt                                                                                                                       <  1:7401313/150‑1 (MQ=255)
                               tCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAAt                                                                                                                       >  1:3041939/1‑150 (MQ=255)
                                    gTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTgcc                                                                                                                  <  2:3296266/150‑1 (MQ=255)
                                      ccACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTgccgc                                                                                                                <  2:4358437/150‑1 (MQ=255)
                                          aTGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATc                                                                                                            <  2:1125914/150‑1 (MQ=255)
                                           tGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCg                                                                                                           <  2:3054715/150‑1 (MQ=255)
                                            gggTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCgg                                                                                                          <  1:862576/150‑1 (MQ=255)
                                            gggTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCgg                                                                                                          <  1:5118380/150‑1 (MQ=255)
                                                 gccgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCt                                                                                                     >  1:9015668/1‑150 (MQ=255)
                                                   cgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTc                                                                                                   <  2:4202264/150‑1 (MQ=255)
                                                   cgccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTc                                                                                                   <  1:8796518/150‑1 (MQ=255)
                                                     ccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGa                                                                                                 <  2:269730/150‑1 (MQ=255)
                                                     ccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGa                                                                                                 <  2:8360572/150‑1 (MQ=255)
                                                     ccAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGa                                                                                                 <  2:4742867/150‑1 (MQ=255)
                                                          aCCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                            <  2:4507831/150‑1 (MQ=255)
                                                            ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                          >  2:3004359/1‑150 (MQ=255)
                                                            ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                          <  1:6282618/150‑1 (MQ=255)
                                                            ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                          >  1:3127442/1‑150 (MQ=255)
                                                            ccGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGccc                                                                                          >  2:762890/1‑150 (MQ=255)
                                                              gTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCgg                                                                                        >  2:8333341/1‑150 (MQ=255)
                                                               tACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGt                                                                                       <  2:3814256/150‑1 (MQ=255)
                                                               tACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGt                                                                                       <  1:4908920/150‑1 (MQ=255)
                                                                      aCCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTgcgc                                                                                >  1:1423129/1‑150 (MQ=255)
                                                                      aCCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTgcgc                                                                                >  2:1796644/1‑150 (MQ=255)
                                                                      aCCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTGCGATACGCCCGGTCTTgcgc                                                                                >  2:7906097/1‑150 (MQ=255)
                                                                        cGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTg                                                                              <  2:7457992/150‑1 (MQ=255)
                                                                         ggTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTgg                                                                             <  1:6658871/150‑1 (MQ=255)
                                                                         ggTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTgg                                                                             <  2:1253889/150‑1 (MQ=255)
                                                                          gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                            <  1:2067148/150‑1 (MQ=255)
                                                                          gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                            <  2:6960079/150‑1 (MQ=255)
                                                                          gTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGt                                                                            <  1:5257108/150‑1 (MQ=255)
                                                                            aTGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGa                                                                          >  2:1105285/1‑150 (MQ=255)
                                                                            aTGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGa                                                                          <  1:7354805/150‑1 (MQ=255)
                                                                              gggCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGc                                                                        <  1:1339056/150‑1 (MQ=255)
                                                                                  tggctgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCg                                                                    >  2:6434644/1‑150 (MQ=255)
                                                                                      tgCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaa                                                                >  2:5379586/1‑150 (MQ=255)
                                                                                       gCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaaa                                                               <  1:6426503/150‑1 (MQ=255)
                                                                                       gCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTaaa                                                               <  1:994068/150‑1 (MQ=255)
                                                                                        ccacccccCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCACTTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCGCCGGTAAAg                                                              >  2:5063613/1‑150 (MQ=255)
                                                                                        ccacccacCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAg                                                              >  1:4208885/1‑150 (MQ=255)
                                                                                          acccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTCAAGgc                                                            >  1:2641975/1‑150 (MQ=255)
                                                                                          acccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgc                                                            <  1:7053512/150‑1 (MQ=255)
                                                                                          acccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgc                                                            >  1:1573717/1‑150 (MQ=255)
                                                                                           cccacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGgcg                                                           >  2:206213/1‑150 (MQ=255)
                                                                                             cacccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGcc                                                         >  1:2257921/1‑150 (MQ=255)
                                                                                              acccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGccc                                                        >  2:2517335/1‑150 (MQ=255)
                                                                                               cccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCt                                                       >  2:700970/1‑150 (MQ=255)
                                                                                               cccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCt                                                       >  2:1455400/1‑150 (MQ=255)
                                                                                                ccGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCtt                                                      >  1:8875305/1‑150 (MQ=255)
                                                                                                 cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                     <  1:7272719/150‑1 (MQ=255)
                                                                                                 cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                     <  2:3574630/150‑1 (MQ=255)
                                                                                                 cGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtc                                                     <  1:2632244/150‑1 (MQ=255)
                                                                                                   cTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTtctc                                                   <  1:5397973/150‑1 (MQ=255)
                                                                                                     cGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCAt                                                 <  1:6555247/150‑1 (MQ=255)
                                                                                                      gATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATg                                                <  2:461688/150‑1 (MQ=255)
                                                                                                      cctgtgCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCACTTCTCATg                                                >  2:3643012/3‑150 (MQ=255)
                                                                                                              tctcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACa                                        <  2:3715167/150‑1 (MQ=255)
                                                                                                                tcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACatt                                      <  1:7686272/150‑3 (MQ=255)
                                                                                                                tcCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGt                                      <  2:8842734/150‑1 (MQ=255)
                                                                                                                   tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                            >  2:8225438/1‑93 (MQ=255)
                                                                                                                   tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGc                                                                                            <  1:8225438/93‑1 (MQ=255)
                                                                                                                   tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                   >  1:93679/1‑150 (MQ=255)
                                                                                                                   tGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCg                                   >  1:8104982/1‑150 (MQ=255)
                                                                                                                   tGCAACACACGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCg                                   <  2:1899793/150‑1 (MQ=255)
                                                                                                                     caacacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGa                                 <  1:9009803/144‑1 (MQ=255)
                                                                                                                      aaGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGAAAAGGCGCCCTTCTCATGGCCGAACAGATCGGAc                                >  1:5852322/1‑150 (MQ=255)
                                                                                                                       acacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGACt                               >  1:1219504/5‑150 (MQ=255)
                                                                                                                        cacaCGCAGCAAGCGCGTCTGCGCCTCCATTGGCATATCGCCGATTTCATCAAGAAACAGCGTGCCGCCATTGGCTTGTTCAAAACGTCCGCTTGAACGGTTCTGTGCGCCAGTAAAGGCACCCTTTTCATGGCCGAACAATTCGGACtc                              <  2:6055493/147‑2 (MQ=255)
                                                                                                                          aaCGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGa                            >  1:578901/1‑150 (MQ=255)
                                                                                                                              cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                        <  1:4183688/150‑1 (MQ=255)
                                                                                                                              cagcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCaa                        <  2:8604180/150‑1 (MQ=255)
                                                                                                                                gcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGt                      >  2:5699667/1‑150 (MQ=255)
                                                                                                                                gcaAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGt                      >  2:8120796/1‑150 (MQ=255)
                                                                                                                                     ttGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCg                 >  1:4270737/1‑150 (MQ=255)
                                                                                                                                       ggCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCggg               <  2:6665848/150‑1 (MQ=255)
                                                                                                                                          cTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGAt            >  1:5050541/1‑150 (MQ=255)
                                                                                                                                           tGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATAGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATc           <  1:3375072/150‑1 (MQ=255)
                                                                                                                                                ttCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcccc      <  2:672588/150‑1 (MQ=255)
                                                                                                                                                ttCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGcccc      <  2:2429185/150‑1 (MQ=255)
                                                                                                                                                   cAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCa   <  2:3394926/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTCGCGGAAACGCCCTTGCTCGATGGCCCGCTCAAGGTCCACATGGGTGGCCGCCAGCACCCGTACATCCACCGGTATGGGCTGGCTGCCACCCACCCGCTCGATGTGCTTCTCCTGCAAGAAACGCAGCAAGTTGGCCTGCAATTCCAACGGCAGGTCGCCGATTTCATCAAGAAACAATGTGCCGCCATCGGCCGCTTCGATACGCCCGGTCTTGCGCTGGTGAGCACCGGTAAAGGCGCCCTTCTCATGGCCGAACAGTTCGGACTGGATCAAGTGCTCGGGGATCGCCCCGCAG  >  NC_002947/3975560‑3975857

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