Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,473,927 C→A A43S (GCC→TCC)  bdcA ← c‑di‑GMP‑binding biofilm dispersal mediator protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,473,9270CA98.6% 455.9 / NA 143A43S (GCC→TCC) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (73/68);  minor base T (1/1);  total (74/69)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

GGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATT  >  NC_000913/4473781‑4474072
                                                                                                                                                  |                                                                                                                                                  
ggCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATc                                                                                                                                                 <  2:2095611/150‑1 (MQ=255)
 gCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCt                                                                                                                                                <  1:510430/150‑1 (MQ=255)
   cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTCTCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCttt                                                                                                                                              <  1:354740/150‑1 (MQ=255)
    tCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCACTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTc                                                                                                                                             <  1:95090/150‑1 (MQ=255)
    tCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTc                                                                                                                                             <  1:426761/150‑1 (MQ=255)
     cGCCAAAGACGCCAATACCTGCATTTACCACCATGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                                                                            <  1:391859/150‑1 (MQ=255)
     cGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCAGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                                                                            >  1:1387416/1‑150 (MQ=255)
      gCCAAAGCCGCCAATACCTGCATTTACCACCAGGATATCCAATGCGACGCTCATACGAACGACATCAATGACAGCCTCTCTGTCATCACAATCTAAGAATACTGATGTAGCTCCAGTCTCTTGTCGCAGGCGTTTATCGGAATCTTTAGa                                                                                                                                           >  1:391650/1‑150 (MQ=255)
      gCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGa                                                                                                                                           >  1:1780981/1‑150 (MQ=255)
       ccAAAGACTCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCTATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAc                                                                                                                                          <  2:91082/150‑1 (MQ=255)
       ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAc                                                                                                                                          >  1:1345087/1‑150 (MQ=255)
       ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAc                                                                                                                                          >  1:1078572/1‑150 (MQ=255)
              cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCAt                                                                                                                                   >  1:1394693/1‑150 (MQ=255)
              cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCAt                                                                                                                                   >  1:2161328/1‑150 (MQ=255)
              cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCAt                                                                                                                                   >  1:765386/1‑150 (MQ=255)
                   aTACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTg                                                                                                                              <  1:261600/150‑1 (MQ=255)
                      ccTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCCGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAAt                                                                                                                           >  2:554250/1‑150 (MQ=255)
                      ccTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAAt                                                                                                                           >  2:648248/1‑150 (MQ=255)
                        tGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACCACATCAAAGACAGCGTCTCTGTCAGCACTATCTGAGAATACTGCTGTCACTCCAGTCTCTTTTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCAAACGTGAAACg                                                                                                                         >  1:280639/1‑150 (MQ=255)
                           aTTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                                                                      <  1:1794800/150‑1 (MQ=255)
                           aTTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                                                                      <  1:918745/150‑1 (MQ=255)
                              tACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCACTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAtt                                                                                                                   >  2:1181894/1‑150 (MQ=255)
                               accaccAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTg                                                                                                                  <  2:617565/150‑1 (MQ=255)
                                ccaccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                                 >  1:1609214/1‑150 (MQ=255)
                                ccaccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                                 >  2:1928647/1‑150 (MQ=255)
                                ccaccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                                 >  2:1305301/1‑150 (MQ=255)
                                 caccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCACTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTaaccgcaa                                                                                                                >  2:1510264/1‑143 (MQ=255)
                                  gccaGGATCTCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                                                                              <  2:771835/149‑1 (MQ=255)
                                  accaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                                                                              >  2:1652476/1‑150 (MQ=255)
                                  accaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                                                                              >  1:2112954/1‑150 (MQ=255)
                                   ccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGccc                                                                                                             <  2:2120756/150‑1 (MQ=255)
                                   ccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTACCGGATACTTACGACCCCCAGTAGGGGAATGGTAAATTGGccc                                                                                                             >  1:1833507/1‑150 (MQ=255)
                                     aGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCa                                                                                                           >  2:224224/1‑150 (MQ=255)
                                     aGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCa                                                                                                           >  2:932184/1‑150 (MQ=255)
                                       gATATCCAATGCGCCGCTCTGACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATTCTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGTATCTTTCGACCCCGGATAGGTGAATCGGACATAGGCCCCATc                                                                                                         >  2:87034/1‑150 (MQ=255)
                                         tatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCgg                                                                                                       <  2:213429/150‑1 (MQ=255)
                                           tCCAATGCGCCGCTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg                                                                                                     >  1:221101/1‑150 (MQ=255)
                                            ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACCTTGGCCCCCATCGGTg                                                                                                     >  1:1051249/1‑150 (MQ=255)
                                            ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGa                                                                                                    >  2:936454/1‑150 (MQ=255)
                                                  cgccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAac                                                                                              <  2:80201/150‑1 (MQ=255)
                                                        tctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTa                                                                                        <  1:1667717/150‑1 (MQ=255)
                                                          ttACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACg                                                                                      <  1:478803/150‑1 (MQ=255)
                                                            acgaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAt                                                                                    <  1:1918251/150‑1 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                <  1:2155689/150‑1 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                >  2:1301025/1‑150 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                >  2:2245453/1‑150 (MQ=255)
                                                                    cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAc                                                                            >  2:1411028/1‑150 (MQ=255)
                                                                    cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAc                                                                            >  1:1078868/1‑150 (MQ=255)
                                                                      tCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCg                                                                          <  1:807981/150‑1 (MQ=255)
                                                                         aTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATa                                                                       <  2:1844633/150‑1 (MQ=255)
                                                                             cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCAc                                                                   <  2:401892/150‑1 (MQ=255)
                                                                                cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGAc                                                                >  1:1017924/1‑150 (MQ=255)
                                                                                       ttCAGCACTAGCTGTGAATCCTGCTGTCGCTCCAGTCTCTTGATCCATGCGTTTAGCTTAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCTGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCAcc                                                         <  2:118707/149‑1 (MQ=255)
                                                                                        tCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCg                                                        >  1:1218761/1‑150 (MQ=255)
                                                                                         cAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATATCGGCACCGATACCACGACTGCCACCGa                                                       >  2:372382/1‑150 (MQ=255)
                                                                                         cAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGa                                                       >  2:2202464/1‑150 (MQ=255)
                                                                                              cTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                  <  1:1896252/150‑1 (MQ=255)
                                                                                               tATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATga                                                 >  2:1934893/1‑150 (MQ=255)
                                                                                                  cTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATATTTCGACCCCGCATAGGTGAATAGTACATTGGCCCCATCGGTGACGAAACTACGTACTATAGCGTCACAGATACCACGACTGCCACCGAGGATgagac                                              >  2:208304/1‑149 (MQ=255)
                                                                                                    gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCCCGACTGCCACCGAGGATGAGAACt                                            >  2:1917343/1‑150 (MQ=255)
                                                                                                     tgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTg                                           <  1:448756/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:910793/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:1780981/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:1162923/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCACCAGTCTCTTGTGCCCGGCTTTTAGCGGTATCTTTCGACCCCGCATAGGTGACTCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:391650/150‑1 (MQ=255)
                                                                                                           ctgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAc                                     <  2:5668/150‑1 (MQ=255)
                                                                                                            tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                    <  2:1609214/150‑1 (MQ=255)
                                                                                                            tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                    <  2:1666759/150‑1 (MQ=255)
                                                                                                            tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                    <  1:504010/150‑1 (MQ=255)
                                                                                                             gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                   <  1:2093214/150‑1 (MQ=255)
                                                                                                             gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGTGAATCGTACATTGGCCCCATCTGTGACGAAACGACGTACGATAGCGGCCCCGATACCCCGACGGCCCCCGAGGTTGAGAAGTGTCTTACCTg                                  >  2:275689/1‑150 (MQ=255)
                                                                                                              ctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                  <  1:224224/150‑1 (MQ=255)
                                                                                                               tgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGt                                 >  1:1488774/1‑150 (MQ=255)
                                                                                                                gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTa                                <  2:765386/150‑1 (MQ=255)
                                                                                                                gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTa                                <  2:760180/150‑1 (MQ=255)
                                                                                                                  cGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGGACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaa                              <  2:520799/150‑1 (MQ=255)
                                                                                                                    cTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAg                            <  1:1065415/150‑1 (MQ=255)
                                                                                                                     tCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                           <  1:2067325/150‑1 (MQ=255)
                                                                                                                     accAGTCTCTTGTGCCAGGCGTGTAGCGGAAGCTTTCGACCCCGCATAGGGGAATCGTACATTGGCCCCATCGGGGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                           <  2:1394693/149‑1 (MQ=255)
                                                                                                                      ccAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcg                          >  1:743770/1‑150 (MQ=255)
                                                                                                                      ccAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcg                          >  1:251713/1‑150 (MQ=255)
                                                                                                                      ccAGGCTCTTGTGCCAGGCGTTTAGCGGAATCTATCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcg                          <  2:1759703/150‑1 (MQ=255)
                                                                                                                        aGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACAACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGcc                        >  1:1169514/1‑150 (MQ=255)
                                                                                                                         gcctctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGATACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                       <  2:1104842/148‑1 (MQ=255)
                                                                                                                         gTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                       <  2:1527447/150‑1 (MQ=255)
                                                                                                                         gTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                       <  2:1028274/150‑1 (MQ=255)
                                                                                                                           ctctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCAt                     >  1:63834/1‑150 (MQ=255)
                                                                                                                               tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                 <  1:936454/150‑1 (MQ=255)
                                                                                                                               tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                 >  2:1945123/1‑150 (MQ=255)
                                                                                                                                  gCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACt              >  1:209809/1‑150 (MQ=255)
                                                                                                                                    cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAATGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAATGACTcc            >  2:1103462/1‑150 (MQ=255)
                                                                                                                                      ggCGTTTAGCGGAATCTTTCTACCACGCATAGGTTAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTTTAAAATCGCCCATAAGGACTCCtt          <  1:1429123/150‑1 (MQ=255)
                                                                                                                                      ggCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCtt          <  2:502644/150‑1 (MQ=255)
                                                                                                                                       gCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTg         <  2:1078868/150‑1 (MQ=255)
                                                                                                                                         gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt       >  1:34664/1‑150 (MQ=255)
                                                                                                                                         gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt       <  1:2034285/150‑1 (MQ=255)
                                                                                                                                          tttAGCGGAAGCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGTATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAtt      <  1:577185/150‑1 (MQ=255)
                                                                                                                                           ttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt     >  2:506370/1‑150 (MQ=255)
                                                                                                                                           ttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt     >  2:1067447/1‑150 (MQ=255)
                                                                                                                                              gCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttatt  <  1:856316/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATT  >  NC_000913/4473781‑4474072

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