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.7% 246.1 / NA 78A43S (GCC→TCC) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (38/39);  total (38/40)
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

GGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTA  >  NC_000913/4473780‑4474070
                                                                                                                                                   |                                                                                                                                               
gggCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAAt                                                                                                                                               >  1:219824/1‑150 (MQ=255)
    cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCttt                                                                                                                                           >  1:969034/1‑150 (MQ=255)
    cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCttt                                                                                                                                           >  1:301501/1‑150 (MQ=255)
    cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCttt                                                                                                                                           <  2:515566/150‑1 (MQ=255)
      cGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAACCTTTCg                                                                                                                                         <  1:416251/150‑1 (MQ=255)
        ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTATCTTGTGCCAGGCGTTTAGCGGAATCTTTCTAc                                                                                                                                       >  1:507305/1‑150 (MQ=255)
          aaaGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAccc                                                                                                                                     <  1:40735/150‑1 (MQ=255)
            aGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCg                                                                                                                                   <  1:354871/150‑1 (MQ=255)
               cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGGCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCAt                                                                                                                                <  1:189210/150‑1 (MQ=255)
                gCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATa                                                                                                                               <  2:881089/150‑1 (MQ=255)
                 ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAg                                                                                                                              >  2:864895/1‑150 (MQ=255)
                  cAATACCTGCATTTACCACCAGGATATCCTATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGGCGCTCCAGTCTCTTGTGCCAGGCGTTAAGCGGAATCTTTCGACCCCGCATAgg                                                                                                                             <  1:131778/150‑1 (MQ=255)
                  cAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGCCCCCGCATAgt                                                                                                                             >  1:179576/1‑149 (MQ=255)
                       ccTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGaaa                                                                                                                        >  1:74718/1‑149 (MQ=255)
                        cTGCATTTACCACCAGGATATCCAATGCGCCGCTATTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCTACACCGCATAGGTCAATc                                                                                                                       >  2:408827/1‑150 (MQ=255)
                        cTGCATTTACCACCAGGATATCCAATGCGCCGCGCTTACGAACCACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATc                                                                                                                       <  2:508341/150‑1 (MQ=255)
                          gCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGt                                                                                                                     >  2:625723/1‑150 (MQ=255)
                          gCATTTACCACCAGGATATCCAAAGCGCCGCTCTTACTAACGACATAAATGACAGCGTCGCTGTCAGCACAATCTGTGAATACTGCTGTCTCTCAAGTCTCTTGTGCAAGGCGTTTATAGGAATCTTTCGACACCGCCTATGTTAATCGt                                                                                                                     >  1:739051/1‑150 (MQ=255)
                           cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCCGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTa                                                                                                                    >  1:700719/1‑150 (MQ=255)
                           cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTa                                                                                                                    <  2:750355/150‑1 (MQ=255)
                                 ccaccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                              >  1:487802/1‑150 (MQ=255)
                                 ccaccaGGATATCCAAAGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                              <  2:403500/150‑1 (MQ=255)
                                   accaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                                                                            >  1:221671/1‑150 (MQ=255)
                                   accaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                                                                            <  2:367055/150‑1 (MQ=255)
                                            tCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg                                                                                                   <  2:998387/150‑1 (MQ=255)
                                             ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGAACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGa                                                                                                  <  2:691612/150‑1 (MQ=255)
                                              cAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCACTGTTAGCACTATCTGTGAATCCTGCTGTCGCTCCAGGCTCTTGTGCCACGCGTTTTGCGGAATCATTCGACCCCGCAGAGGTGAATCGTACATTGGCCCCATCGGTgtc                                                                                                 >  2:757525/1‑148 (MQ=255)
                                               aaTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACg                                                                                                <  1:850606/150‑1 (MQ=255)
                                                  gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGAGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa                                                                                             <  2:94857/150‑1 (MQ=255)
                                                       gcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGCACATTGGCCCCATCGGTGACGAAacgacg                                                                                        <  2:1019610/150‑1 (MQ=255)
                                                         tctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTa                                                                                      >  1:632882/1‑150 (MQ=255)
                                                           ttACGAACGACATCAATGACAGCGTCTCTGTCAGCACTTTCTGTGAATACTGCTGTCGCTCCAGTCTCTTTTTCCAGTCGTTTATAGGATTCTTTCGACCCCGCATAGGTGAATCGTACCTTGGCCGCTTCGGTGACGAAACGACGTACg                                                                                    <  1:647175/150‑1 (MQ=255)
                                                            tACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCAGACGCTCCAGTCTCTTGTGCCAGGCCTTTAGCGGAATCTTTCGACCACCCAAAGGTGAATCGTACAATGGCCCCATCCGTGACGAAACGACGGACGa                                                                                   >  1:927111/1‑150 (MQ=255)
                                                              cgaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATa                                                                                 <  2:939908/150‑1 (MQ=255)
                                                               gaacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAg                                                                                >  1:359420/1‑150 (MQ=255)
                                                                aacgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTTAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGc                                                                               <  1:939141/150‑1 (MQ=255)
                                                                 acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                              <  2:365606/150‑1 (MQ=255)
                                                                 acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCg                                                                              <  2:219824/150‑1 (MQ=255)
                                                                   gaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGGTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGTTGACGAAACGACGTACGATAGCGGc                                                                            <  1:809964/150‑1 (MQ=255)
                                                                      aTCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAcc                                                                         <  2:394871/150‑1 (MQ=255)
                                                                       tCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCg                                                                        >  2:182095/1‑150 (MQ=255)
                                                                          aTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCCGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATa                                                                     <  2:152204/150‑1 (MQ=255)
                                                                            gACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCAAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATAcc                                                                   >  1:705730/1‑150 (MQ=255)
                                                                                gCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTATGATAGCGGCACCGATACCACGa                                                               >  2:681400/1‑150 (MQ=255)
                                                                                   tctcTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGATTCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTg                                                            >  2:296493/1‑150 (MQ=255)
                                                                                    ctctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGc                                                           >  1:366775/1‑150 (MQ=255)
                                                                                    ctctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGCAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGc                                                           >  1:1009116/1‑150 (MQ=255)
                                                                                       tGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCAc                                                        <  1:625723/150‑1 (MQ=255)
                                                                                                  tCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATgaga                                             <  1:15910/150‑1 (MQ=255)
                                                                                                    tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGAGGAGAAc                                           >  2:435301/1‑150 (MQ=255)
                                                                                                          tACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCtt                                     <  2:921932/150‑1 (MQ=255)
                                                                                                           aCTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTa                                    <  2:505417/150‑1 (MQ=255)
                                                                                                           aCTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTa                                    <  2:65375/150‑1 (MQ=255)
                                                                                                              gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCTCATAGGTGAATCGTACATAGGCCCCATCGGTGACGAAACGACGTACGAAAGCCGCACCGATACCACGACTGCCCCCGAGGATGAGAACAGTCTTaaaa                                 >  1:475860/1‑147 (MQ=255)
                                                                                                              gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCTTCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                 >  2:707543/1‑150 (MQ=255)
                                                                                                              gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                 >  1:718970/1‑150 (MQ=255)
                                                                                                              gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCt                                 <  1:581516/150‑1 (MQ=255)
                                                                                                                tgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACAACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTGTTACCTGt                               >  2:182700/1‑150 (MQ=255)
                                                                                                                 ggcGCTCCAGGCTTTTGTGCCAGGCGGTTAGCGGAATCTTTCGACCCCGCATAGTTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTa                              <  2:700719/148‑1 (MQ=255)
                                                                                                                 gTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTa                              >  2:734181/1‑150 (MQ=255)
                                                                                                                    gCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTaaaa                           >  2:331110/1‑150 (MQ=255)
                                                                                                                         tgtgtctTGGGCCAGGCGTTTAGCGTCAGCTTTGGACCCCGCTTAGGTGAATCGTACATTGGCCCGATCGGTGTCGAAACGTCGTACGATAGCGGTACCGATACCACGACTGCCACCGAGGATGAGAACTATCTTACCTGTATTAGCGcc                      <  2:927111/146‑1 (MQ=255)
                                                                                                                          gTCTCTTTTGCCAGGCGTTTAGCGGAATCTTTAGAACACGCATAGGTTAATCGTACATTTGCCACATCGGTGAAGAAACGCCGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAAATGTCTTACCTGTAAAAGCGccc                     <  2:739051/150‑1 (MQ=255)
                                                                                                                          gTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                     <  2:487802/150‑1 (MQ=255)
                                                                                                                           tctcTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGCTACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCa                    >  2:246412/1‑150 (MQ=255)
                                                                                                                             tctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                  <  1:268309/150‑1 (MQ=255)
                                                                                                                             tctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATa                  <  2:921191/150‑1 (MQ=255)
                                                                                                                                tgtgCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg               >  2:490656/1‑150 (MQ=255)
                                                                                                                                   gCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACt            >  1:872752/1‑150 (MQ=255)
                                                                                                                                   gCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACt            <  1:136870/150‑1 (MQ=255)
                                                                                                                                     cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGGACATTGGCCCCATCGGTGACGAAAAGACGTCCGATAGCGGCACCGATACCCCGCCTGCCACCGAGGAGGAGAACTGGCTTACCTGTAAACGCGCCAATAAGGAGTcc          >  2:76342/1‑150 (MQ=255)
                                                                                                                                      aGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCt         >  2:911385/1‑150 (MQ=255)
                                                                                                                                          gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt     >  1:728738/1‑150 (MQ=255)
                                                                                                                                          gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt     <  2:359420/150‑1 (MQ=255)
                                                                                                                                          gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAt     <  2:728738/150‑1 (MQ=255)
                                                                                                                                           tttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAtt    <  1:652890/150‑1 (MQ=255)
                                                                                                                                            ttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCAAAAGGACTCCTTGAttt   >  2:138431/1‑150 (MQ=255)
                                                                                                                                             tAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATtta  >  2:317124/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTA  >  NC_000913/4473780‑4474070

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