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.

GGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATT  >  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)
                                  gccaGGATCTCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑cc                                                                                                              <  2:771835/149‑1 (MQ=255)
                                  accaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑cc                                                                                                              >  2:1652476/1‑150 (MQ=255)
                                  accaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑cc                                                                                                              >  1:2112954/1‑150 (MQ=255)
                                   ccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑ccc                                                                                                             <  2:2120756/150‑1 (MQ=255)
                                   ccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTACCGGATACTTACGACCCCCAGTAGGGGAATGGTAAATTGG‑ccc                                                                                                             >  1:1833507/1‑150 (MQ=255)
                                     aGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCa                                                                                                           >  2:224224/1‑150 (MQ=255)
                                     aGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCa                                                                                                           >  2:932184/1‑150 (MQ=255)
                                       gATATCCAATGCGCCGCTCTGACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATTCTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGTATCTTTCGACCCCGGATAGGTGAATCGGACATAGG‑CCCCATc                                                                                                         >  2:87034/1‑150 (MQ=255)
                                         tatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCgg                                                                                                       <  2:213429/150‑1 (MQ=255)
                                           tCCAATGCGCCGCTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTg                                                                                                     >  1:221101/1‑150 (MQ=255)
                                            ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACCTTGGCCCCCATCGGTg                                                                                                     >  1:1051249/1‑150 (MQ=255)
                                            ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGa                                                                                                    >  2:936454/1‑150 (MQ=255)
                                                  cgccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAac                                                                                              <  2:80201/150‑1 (MQ=255)
                                                        tctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTa                                                                                        <  1:1667717/150‑1 (MQ=255)
                                                          ttACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACg                                                                                      <  1:478803/150‑1 (MQ=255)
                                                            acgaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGAt                                                                                    <  1:1918251/150‑1 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                <  1:2155689/150‑1 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                >  2:1301025/1‑150 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCg                                                                                >  2:2245453/1‑150 (MQ=255)
                                                                    cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCAc                                                                            >  2:1411028/1‑150 (MQ=255)
                                                                    cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCAc                                                                            >  1:1078868/1‑150 (MQ=255)
                                                                      tCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCg                                                                          <  1:807981/150‑1 (MQ=255)
                                                                         aTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATa                                                                       <  2:1844633/150‑1 (MQ=255)
                                                                             cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCAc                                                                   <  2:401892/150‑1 (MQ=255)
                                                                                cGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGAc                                                                >  1:1017924/1‑150 (MQ=255)
                                                                                       ttCAGCACTAGCTGTGAATCCTGCTGTCGCTCCAGTCTCTTGATCCATGCGTTTAGCTTAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCTGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCAcc                                                         <  2:118707/149‑1 (MQ=255)
                                                                                        tCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCg                                                        >  1:1218761/1‑150 (MQ=255)
                                                                                         cAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATATCGGCACCGATACCACGACTGCCACCGa                                                       >  2:372382/1‑150 (MQ=255)
                                                                                         cAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGa                                                       >  2:2202464/1‑150 (MQ=255)
                                                                                              cTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                  <  1:1896252/150‑1 (MQ=255)
                                                                                               tATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATga                                                 >  2:1934893/1‑150 (MQ=255)
                                                                                                  cTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATATTTCGACCCCGCATAGGTGAATAGTACATTGG‑CCCCATCGGTGACGAAACTACGTACTATAGCGTCACAGATACCACGACTGCCACCGAGGATgagac                                              >  2:208304/1‑149 (MQ=255)
                                                                                                    gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCCCGACTGCCACCGAGGATGAGAACt                                            >  2:1917343/1‑150 (MQ=255)
                                                                                                     tgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTg                                           <  1:448756/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:910793/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:1780981/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:1162923/150‑1 (MQ=255)
                                                                                                       aaTACTGCTGTCGCACCAGTCTCTTGTGCCCGGCTTTTAGCGGTATCTTTCGACCCCGCATAGGTGACTCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTc                                         <  2:391650/150‑1 (MQ=255)
                                                                                                           ctgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAc                                     <  2:5668/150‑1 (MQ=255)
                                                                                                            tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                    <  2:1609214/150‑1 (MQ=255)
                                                                                                            tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                    <  2:1666759/150‑1 (MQ=255)
                                                                                                            tgctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTAcc                                    <  1:504010/150‑1 (MQ=255)
                                                                                                             gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                   <  1:2093214/150‑1 (MQ=255)
                                                                                                             gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATA‑GTGAATCGTACATTGG‑CCCCATCTGTGACGAAACGACGTACGATAGCGGCCCCGATACCCCGACGGCCCCCGAGGTTGAGAAGTGTCTTACCTg                                  >  2:275689/1‑150 (MQ=255)
                                                                                                              ctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                  <  1:224224/150‑1 (MQ=255)
                                                                                                               tgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGt                                 >  1:1488774/1‑150 (MQ=255)
                                                                                                                gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTa                                <  2:765386/150‑1 (MQ=255)
                                                                                                                gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTa                                <  2:760180/150‑1 (MQ=255)
                                                                                                                  cGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGGACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaa                              <  2:520799/150‑1 (MQ=255)
                                                                                                                    cTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTTAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAg                            <  1:1065415/150‑1 (MQ=255)
                                                                                                                     tCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                           <  1:2067325/150‑1 (MQ=255)
                                                                                                                     accAGTCTCTTGTGCCAGGCGTGTAGCGGAAGCTTTCGACCCCGCATAGGGGAATCGTACATTGG‑CCCCATCGGGGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                           <  2:1394693/149‑1 (MQ=255)
                                                                                                                      ccAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcg                          >  1:743770/1‑150 (MQ=255)
                                                                                                                      ccAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcg                          >  1:251713/1‑150 (MQ=255)
                                                                                                                      ccAGGCTCTTGTGCCAGGCGTTTAGCGGAATCTATCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgcg                          <  2:1759703/150‑1 (MQ=255)
                                                                                                                        aGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACAACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGcc                        >  1:1169514/1‑150 (MQ=255)
                                                                                                                         gcctctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGATACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                       <  2:1104842/148‑1 (MQ=255)
                                                                                                                         gTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                       <  2:1527447/150‑1 (MQ=255)
                                                                                                                         gTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                       <  2:1028274/150‑1 (MQ=255)
                                                                                                                           ctctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCAt                     >  1:63834/1‑150 (MQ=255)
                                                                                                                               tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                 <  1:936454/150‑1 (MQ=255)
                                                                                                                               tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                 >  2:1945123/1‑150 (MQ=255)
                                                                                                                                  gCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACt              >  1:209809/1‑150 (MQ=255)
                                                                                                                                    cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAATGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAATGACTcc            >  2:1103462/1‑150 (MQ=255)
                                                                                                                                      ggCGTTTAGCGGAATCTTTCTACCACGCATAGGTTAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTTTAAAATCGCCCATAAGGACTCCtt          <  1:1429123/150‑1 (MQ=255)
                                                                                                                                      ggCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCtt          <  2:502644/150‑1 (MQ=255)
                                                                                                                                       gCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTg         <  2:1078868/150‑1 (MQ=255)
                                                                                                                                         gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt       >  1:34664/1‑150 (MQ=255)
                                                                                                                                         gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt       <  1:2034285/150‑1 (MQ=255)
                                                                                                                                          tttAGCGGAAGCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGTATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAtt      <  1:577185/150‑1 (MQ=255)
                                                                                                                                           ttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt     >  2:506370/1‑150 (MQ=255)
                                                                                                                                           ttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt     >  2:1067447/1‑150 (MQ=255)
                                                                                                                                              gCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttatt  <  1:856316/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGG‑CCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATT  >  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: ‑