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,9270CA99.4% 531.6 / ‑3.9 158A43S (GCC→TCC) bdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (86/71);  total (86/72)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

Only 100 of 158 total aligned reads displayed.

GGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTA  >  NC_000913/4473781‑4474073
                                                                                                                                                  |                                                                                                                                                  
ggCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCt                                                                                                                                                <  1:102828/151‑1 (MQ=255)
  ccTCGCCAAAGACGCCAATACCTGCATTTACTACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCttt                                                                                                                                              >  1:2205967/1‑151 (MQ=255)
  ccTCGCCAAAGACGCCAAAACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCGTACGAACAACATCAATTACAGCGTCCCAGTCAGCACTAACTGTGAATACTCCTGTCGCCCCAGTCCCTTGTGCCAGGCGTTTAGCGGAATCtct                                                                                                                                              >  2:851059/1‑149 (MQ=255)
   cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTc                                                                                                                                             <  2:2677105/151‑1 (MQ=255)
   cTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTc                                                                                                                                             >  2:345091/1‑151 (MQ=255)
    tCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCCATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                                                                            <  1:2650169/151‑1 (MQ=255)
    tCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                                                                            >  2:1197926/1‑151 (MQ=255)
    tCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCg                                                                                                                                            <  1:311785/151‑1 (MQ=255)
     caccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGa                                                                                                                                           >  1:1860500/3‑151 (MQ=255)
     cGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGa                                                                                                                                           >  2:1174767/1‑151 (MQ=255)
     cGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGa                                                                                                                                           <  1:1818454/151‑1 (MQ=255)
      gCCAAAGACGCCAATACCTGCATTTCCCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAc                                                                                                                                          <  2:2391220/151‑1 (MQ=255)
       ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAcc                                                                                                                                         >  1:1929380/1‑151 (MQ=255)
       ccAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCCCTTGTGCCAGGCGTTTAGCGGAATCTTTCGAcc                                                                                                                                         >  2:659825/1‑151 (MQ=255)
          aaGACCCCAATACCTTCATTTACCACCAGAATATCCACTGCGCCGCTCTTACGCACCACATCAATCACATCGTCTCTGTCACCACTATCTGTCAATACTGCTGTCACTCCAGTCTCTTGTGCCAGGCGTTTCGCGGAATCTTTCGAcacca                                                                                                                                      >  2:213139/1‑150 (MQ=255)
           aGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGc                                                                                                                                     >  1:2415202/1‑151 (MQ=255)
           aGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGc                                                                                                                                     >  2:508282/1‑151 (MQ=255)
              cGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATa                                                                                                                                  >  2:2030562/1‑151 (MQ=255)
                ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAgg                                                                                                                                >  1:1652990/1‑151 (MQ=255)
                ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGATCGAGTCTCTTGTGCCAGGCGTTTAGGGGAATCTTTCGACCCCGCATggg                                                                                                                                >  2:2229963/1‑151 (MQ=255)
                ccAATACCTGCATTTACCACCAGGAAATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATATTTCGACCCCGCATAgg                                                                                                                                >  2:1895698/1‑151 (MQ=255)
                  aaTACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTg                                                                                                                              <  2:885138/151‑1 (MQ=255)
                  aaTACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTg                                                                                                                              <  2:1667073/151‑1 (MQ=255)
                    tACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGaa                                                                                                                            >  1:565808/1‑151 (MQ=255)
                     aCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAAt                                                                                                                           >  2:1068263/1‑151 (MQ=255)
                      ccTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATc                                                                                                                          >  1:421665/1‑151 (MQ=255)
                      ccTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATc                                                                                                                          >  1:882414/1‑151 (MQ=255)
                       cTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCg                                                                                                                         <  2:2613960/151‑1 (MQ=255)
                          cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTAc                                                                                                                      <  1:2344649/151‑1 (MQ=255)
                           aTTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCCCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACa                                                                                                                     <  2:2660692/151‑1 (MQ=255)
                            tttACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACAt                                                                                                                    >  2:2101307/1‑151 (MQ=255)
                               accaccAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTgg                                                                                                                 <  1:882517/151‑1 (MQ=255)
                                 caccaGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGcc                                                                                                               <  1:963012/151‑1 (MQ=255)
                                     aGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                                                                                                           >  2:2344823/1‑151 (MQ=255)
                                         tatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGGCCCCGCATAGGTGAATCGTACATTGGCCCCATCGGt                                                                                                       <  2:933595/151‑1 (MQ=255)
                                         tatCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGt                                                                                                       <  1:1292460/151‑1 (MQ=255)
                                            ccAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGAc                                                                                                    <  1:959437/151‑1 (MQ=255)
                                             cAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACg                                                                                                   >  2:731078/1‑151 (MQ=255)
                                                  cgccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacg                                                                                              >  2:2661758/1‑151 (MQ=255)
                                                   gccgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacga                                                                                             <  1:113495/151‑1 (MQ=255)
                                                     cgcTCTTACGAACGACATCAATGCCAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgacg                                                                                           <  1:1433810/151‑1 (MQ=255)
                                                     cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCATTTGTGGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgacg                                                                                           <  1:1769290/151‑1 (MQ=255)
                                                     cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgacg                                                                                           <  1:1853789/151‑1 (MQ=255)
                                                        tctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTAc                                                                                        <  2:303190/151‑1 (MQ=255)
                                                                acgaCATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCgg                                                                                <  2:1251334/151‑1 (MQ=255)
                                                                    cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAcc                                                                            <  2:419380/151‑1 (MQ=255)
                                                                    cATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCAcc                                                                            >  2:2464317/1‑151 (MQ=255)
                                                                     aTCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGNTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAACGGCACCg                                                                           >  1:722277/1‑151 (MQ=255)
                                                                       cAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGAt                                                                         <  2:1310552/151‑1 (MQ=255)
                                                                       cAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGAt                                                                         >  2:1304062/1‑151 (MQ=255)
                                                                        aaTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATa                                                                        >  1:1821940/1‑151 (MQ=255)
                                                                        aaTGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAAAGCGGCACCGATa                                                                        >  1:2228387/1‑151 (MQ=255)
                                                                            ccAGCGTCTCTGCCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAACCGACGTACGATAGCGGCACCGATACCAc                                                                    <  2:159972/150‑1 (MQ=255)
                                                                             cAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACg                                                                   >  1:755006/1‑151 (MQ=255)
                                                                                   ctctGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGcc                                                             >  2:1673125/1‑151 (MQ=255)
                                                                                       gTCAGCATTATCTGTGAATACTGCTGTCGCTTCAGTCTTTTGTGCCAGGCGTTTAGCGAAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCTGTGACGAAACGACGTACGATATCGTCTTTGATAACACGACTGCCACCg                                                         <  2:1596043/151‑1 (MQ=255)
                                                                                         cAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAg                                                       <  2:445388/151‑1 (MQ=255)
                                                                                          aGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAgg                                                      <  1:1702102/151‑1 (MQ=255)
                                                                                           gCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAggg                                                     >  2:1467877/1‑150 (MQ=255)
                                                                                           gCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGa                                                     >  1:1352695/1‑151 (MQ=255)
                                                                                            cACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTACCACCGAGGAt                                                    <  1:492288/151‑1 (MQ=255)
                                                                                             aCTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATg                                                   <  1:1297459/151‑1 (MQ=255)
                                                                                              cTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATga                                                  <  1:83061/151‑1 (MQ=255)
                                                                                                aTCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATgaga                                                >  1:1907525/1‑151 (MQ=255)
                                                                                                  cTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAAc                                              <  1:2030590/151‑1 (MQ=255)
                                                                                                  cTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAAc                                              <  2:2415503/151‑1 (MQ=255)
                                                                                                  cTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAAc                                              <  2:890518/151‑1 (MQ=255)
                                                                                                  cTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAAc                                              <  1:359940/151‑1 (MQ=255)
                                                                                                   tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACt                                             >  1:2403074/1‑151 (MQ=255)
                                                                                                   tgtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACt                                             >  1:1669801/1‑151 (MQ=255)
                                                                                                    gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTg                                            <  2:1893726/151‑1 (MQ=255)
                                                                                                    gtgAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTg                                            <  1:2079454/151‑1 (MQ=255)
                                                                                                        aTACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCtt                                        <  1:1199625/151‑1 (MQ=255)
                                                                                                         tACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTa                                       <  1:1189753/151‑1 (MQ=255)
                                                                                                             gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                   <  2:2481396/151‑1 (MQ=255)
                                                                                                             gctgTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTg                                   <  2:1268363/151‑1 (MQ=255)
                                                                                                                    cTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                            >  1:1208509/1‑151 (MQ=255)
                                                                                                                    cTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                            >  1:2037867/1‑151 (MQ=255)
                                                                                                                    cTCCAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAgc                            >  1:1628913/1‑151 (MQ=255)
                                                                                                                      ccAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGCTGAGAACTGTCTTACCTGTAAAAgcgc                          >  2:930252/1‑151 (MQ=255)
                                                                                                                       cAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGcc                         >  1:2103376/1‑151 (MQ=255)
                                                                                                                       cAGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGcc                         <  1:1568017/151‑1 (MQ=255)
                                                                                                                        aGTCTCTTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGccc                        <  1:422068/151‑1 (MQ=255)
                                                                                                                            tctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATaa                    <  2:1852140/151‑1 (MQ=255)
                                                                                                                             ctTGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTAAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGAAAGCGGTACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCCTAAg                   >  2:626253/1‑151 (MQ=255)
                                                                                                                              ttGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                  <  1:1019590/151‑1 (MQ=255)
                                                                                                                              ttGTGCCAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAgg                  <  2:2415166/151‑1 (MQ=255)
                                                                                                                                   ccAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTcc             <  1:7519/151‑1 (MQ=255)
                                                                                                                                   ccAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTcc             <  1:2368174/151‑1 (MQ=255)
                                                                                                                                    cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCt            >  1:714418/1‑151 (MQ=255)
                                                                                                                                    cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCt            >  1:75671/1‑151 (MQ=255)
                                                                                                                                    cAGGCGTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCt            >  1:488676/1‑151 (MQ=255)
                                                                                                                                         gTTTAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAtt       <  1:1082375/151‑1 (MQ=255)
                                                                                                                                          tttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt      <  1:1068278/151‑1 (MQ=255)
                                                                                                                                          tttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGAttt      <  1:1150336/151‑1 (MQ=255)
                                                                                                                                           ttAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATtta     >  1:1453046/1‑151 (MQ=255)
                                                                                                                                            tAGCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttat    >  1:2178699/1‑151 (MQ=255)
                                                                                                                                              gtggCACCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttatta  <  1:847347/149‑1 (MQ=255)
                                                                                                                                              gCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATACGGACTCCTTGATttatta  >  2:477839/1‑151 (MQ=255)
                                                                                                                                              gCGGAATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATttatta  <  2:2008115/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACGTACGATAGCGGCACCGATACCACGACTGCCACCGAGGATGAGAACTGTCTTACCTGTAAAAGCGCCCATAAGGACTCCTTGATTTATTA  >  NC_000913/4473781‑4474073

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