Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,473,847 G→T 5.7% V69V (GTC→GTA 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,8470GT 5.7% 235.4 / 2.5 71V69V (GTC→GTAbdcAc‑di‑GMP‑binding biofilm dispersal mediator protein
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base T (2/2);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACG  >  NC_000913/4473710‑4473984
                                                                                                                                         |                                                                                                                                         
aCAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAATGGCGATCAATTTCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGAc                                                                                                                                         <  2:1452400/140‑1 (MQ=255)
aCAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGAc                                                                                                                                         <  2:599568/140‑1 (MQ=255)
aCAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGAc                                                                                                                                         <  2:151034/140‑1 (MQ=255)
     ggCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCaa                                                                                                                                    >  2:668567/1‑140 (MQ=255)
     ggCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCaa                                                                                                                                    >  1:713541/1‑140 (MQ=255)
     ggCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCaa                                                                                                                                    >  2:163/1‑140 (MQ=255)
     ggCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCaa                                                                                                                                    >  2:669053/1‑140 (MQ=255)
         tGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGAc                                                                                                                                >  1:711392/1‑140 (MQ=255)
          gATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACa                                                                                                                               <  1:169345/140‑1 (MQ=255)
          gATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACa                                                                                                                               <  1:169393/140‑1 (MQ=255)
               ggATCAAGAATATTCAATTTGAAAAGGCTATCAATATCGTCTGAAATTAATTCCACCGCCTCGCCAAAGACTCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCTCTCTTACGAAAGACATCAATGACAGCGtc                                                                                                                          <  2:1381373/140‑1 (MQ=255)
               ggAGCATGAATATTTATTTTGAAATGGCGAACAATAACGTCGGCATTTAATTCCAGGGCCTCGCGATAGACGCCAATTCCTGCAATTACCTCCACGATATCCAATGCGCCTCTCTTACGAACCACATCAATGACAGCGta                                                                                                                          <  2:331627/140‑2 (MQ=255)
               ggAGCATGAATATTGTTTTAGAAAAGGCGATCCATATCGTCGGCATTTACATCCAGGGCCTCGCCTAAGACGCCACTACCTGCATATACCACCAGGATAGCCAATGCGCCGCTCTTACGAACGACCTCACTGACAGCGtc                                                                                                                          <  2:373999/140‑1 (MQ=255)
               ggAGCATGAATATTGATTATGAAAAGGCGATCAATAACGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGtc                                                                                                                          <  2:73586/140‑1 (MQ=255)
               cgAGCACGACTATGGATTTTGAAAAGGCGAGAAATATCGTCGGCATTTACGACCAGCGCCTCGCGAAAGACCCCACGCCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTCCGAACGACATCAATGACAGCGtc                                                                                                                          <  2:331683/139‑1 (MQ=255)
                gAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGtct                                                                                                                         >  1:790342/1‑140 (MQ=255)
                   cATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTg                                                                                                                      >  1:463067/1‑140 (MQ=255)
                    aTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCATTGACAGCGTctcttt                                                                                                                     >  2:683331/1‑138 (MQ=255)
                    aTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGt                                                                                                                     >  2:1341987/1‑140 (MQ=255)
                    aTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGt                                                                                                                     >  2:1007512/1‑140 (MQ=255)
                    aTGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAAATCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGt                                                                                                                     >  2:722070/1‑140 (MQ=255)
                          atTGATTTTGAAAAGGCGATCACTATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAAAGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCAc                                                                                                               >  2:142393/1‑140 (MQ=255)
                          atTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCAc                                                                                                               <  1:1627352/140‑1 (MQ=255)
                           ttGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACTACATCAATGACAGCGTCTCTGTCAGCACt                                                                                                              >  1:339622/1‑140 (MQ=255)
                           ttGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACTACATCAATGACAGCGTCTCTGTCAGCACt                                                                                                              >  1:339529/1‑140 (MQ=255)
                            tGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTa                                                                                                             <  1:781018/140‑1 (MQ=255)
                            tGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGACTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTa                                                                                                             <  2:206447/140‑1 (MQ=255)
                             gATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTAt                                                                                                            >  1:789962/1‑140 (MQ=255)
                                   gAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGa                                                                                                      >  2:1157204/1‑140 (MQ=255)
                                   gAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGa                                                                                                      >  2:1157419/1‑140 (MQ=255)
                                             tCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACTACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGt                                                                                            <  2:339622/140‑1 (MQ=255)
                                             tCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACTACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGt                                                                                            <  2:339529/140‑1 (MQ=255)
                                             tCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGt                                                                                            <  2:1126075/140‑1 (MQ=255)
                                             tCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGt                                                                                            <  2:777483/140‑1 (MQ=255)
                                              cAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCCCCAGGATATCCAATGCGCCGCTCTTACTACCGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGGTGTc                                                                                           <  2:277293/140‑1 (MQ=255)
                                                    cgtcgGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCa                                                                                     >  2:696426/1‑140 (MQ=255)
                                                        ggCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGtct                                                                                 <  1:1556785/140‑1 (MQ=255)
                                                                aTTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGcc                                                                         <  2:752887/140‑1 (MQ=255)
                                                                    cAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGc                                                                     >  1:390284/1‑140 (MQ=255)
                                                                           tCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGc                                                              >  2:917148/1‑140 (MQ=255)
                                                                           tCGCCAAAGACGCCAAAACCTGCAATTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGACTCTGTCAGCACTATCTGTGAAAACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCCTTTAGc                                                              >  2:914734/1‑140 (MQ=255)
                                                                            cGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCg                                                             >  1:1109516/1‑140 (MQ=255)
                                                                            cGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCg                                                             >  1:1414634/1‑140 (MQ=255)
                                                                                       ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGa                                                  <  1:1611276/140‑1 (MQ=255)
                                                                                       ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGa                                                  <  1:959368/140‑1 (MQ=255)
                                                                                       ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGa                                                  <  1:122522/140‑1 (MQ=255)
                                                                                       ccAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGa                                                  <  2:1346922/140‑1 (MQ=255)
                                                                                         aaTACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGAcc                                                <  2:894724/140‑1 (MQ=255)
                                                                                            cccagCATTTACCACCAGGATATCCACTGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCg                                             <  1:1672963/136‑1 (MQ=255)
                                                                                            acatGCATTTACGAACAGGATATCCAATGCGCCGCTATTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTTACTGTGAATACTGCTGTCGCTCCAGTCTGTTTTGCAAGGCGTTTAGCGGCATCTTTCGACCCCg                                             <  2:357470/137‑1 (MQ=255)
                                                                                            aCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCg                                             <  2:424846/140‑1 (MQ=255)
                                                                                            aCCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTATCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCg                                             <  2:1376902/140‑1 (MQ=255)
                                                                                              cTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATTCTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCa                                           <  1:815443/140‑1 (MQ=255)
                                                                                                 cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAg                                        <  2:850892/140‑1 (MQ=255)
                                                                                                 cATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAg                                        <  2:602149/140‑1 (MQ=255)
                                                                                                   tttACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGt                                      <  1:1166083/140‑1 (MQ=255)
                                                                                                      accaccAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGaa                                   >  2:912738/1‑140 (MQ=255)
                                                                                                      accaccAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGaa                                   >  2:912802/1‑140 (MQ=255)
                                                                                                                  tCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGc                       >  1:211614/1‑140 (MQ=255)
                                                                                                                       tGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                  >  2:950916/1‑140 (MQ=255)
                                                                                                                       tGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                  >  2:21125/1‑140 (MQ=255)
                                                                                                                       tGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAt                  >  2:21114/1‑140 (MQ=255)
                                                                                                                        gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATc                 >  1:1561695/1‑140 (MQ=255)
                                                                                                                        gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATc                 >  1:1562164/1‑140 (MQ=255)
                                                                                                                        gcgcCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATc                 >  1:911347/1‑140 (MQ=255)
                                                                                                                            cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg             >  2:500133/1‑140 (MQ=255)
                                                                                                                            cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg             >  1:440797/1‑140 (MQ=255)
                                                                                                                            cgcTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTg             >  2:642579/1‑140 (MQ=255)
                                                                                                                              ctctTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATAGGCCCCATCGGTGAc           >  1:1574210/1‑140 (MQ=255)
                                                                                                                                  tACGAACGACATCAATGACAGCGACTCTGTCAGCACTATCAGTGAAGACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCAACTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa       <  2:779493/140‑1 (MQ=255)
                                                                                                                                  tACGAACGACAACAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCAACGGTGACGaaa       <  2:779407/140‑1 (MQ=255)
                                                                                                                                  aacgaacgaCATCAATCACAGCGTCTCTGTCAGCACTATCAGTGAATACTGCTGTCGCGCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGaaa       <  2:1225012/139‑1 (MQ=255)
                                                                                                                                       ccgaCATCACTGACAGCGTCTCTGTCAGCACGATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGGTTAGCGGCATCTTTCGACCCCGCCGAGGTGAATCGTACATTGGCCCCATCGGTGACGAAacgacg  <  2:135810/139‑1 (MQ=255)
                                                                                                                                       acgaCATCAAAGACAGCGTCTCTGTCAGCACTATCACTGAATACTGCTGTCGCTCCAGTCTCTTGTTCCAGGCGTTTAGCGGCATCTTTCGACCCAGCATAGGTGTATCTTACATTGGCCCCATCGGTGAAGAAacgacg  <  2:620525/140‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
ACAGAGGCATGATAAGGAGCATGAATATTGATTTTGAAAAGGCGATCAATATCGTCGGCATTTAATTCCAGGGCCTCGCCAAAGACGCCAATACCTGCATTTACCACCAGGATATCCAATGCGCCGCTCTTACGAACGACATCAATGACAGCGTCTCTGTCAGCACTATCTGTGAATACTGCTGTCGCTCCAGTCTCTTGTGCCAGGCGTTTAGCGGCATCTTTCGACCCCGCATAGGTGAATCGTACATTGGCCCCATCGGTGACGAAACGACG  >  NC_000913/4473710‑4473984

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