Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,170,282 C→G 6.9% W31C (TGG→TGC ldtC ← L,D‑transpeptidase linking Lpp to murein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,170,2820CG 6.9% 199.5 / 6.4 58W31C (TGG→TGCldtCL,D‑transpeptidase linking Lpp to murein
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTA  >  NC_000913/1170136‑1170428
                                                                                                                                                  |                                                                                                                                                  
tAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCa                                                                                                                                                  >  2:282328/1‑149 (MQ=255)
tAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGCCGCCAGCGGGCa                                                                                                                                                  <  3:350554/149‑1 (MQ=255)
 aaGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAg                                                                                                                                                 <  4:269239/149‑1 (MQ=255)
    gATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGgtg                                                                                                                                              >  4:165042/1‑149 (MQ=255)
     aTCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACATCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            <  2:287393/149‑1 (MQ=255)
     aTCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACATCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            <  4:809094/149‑1 (MQ=255)
      tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            <  3:165042/149‑1 (MQ=255)
      tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            <  4:472797/149‑1 (MQ=255)
      tCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            >  3:209772/1‑149 (MQ=255)
       cAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTg                                                                                                                                           <  1:98529/149‑1 (MQ=255)
         aCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGATGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGc                                                                                                                                         >  2:27948/1‑149 (MQ=255)
         aCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGc                                                                                                                                         >  2:667569/1‑149 (MQ=255)
         aCGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGc                                                                                                                                         <  3:777458/149‑1 (MQ=255)
                    aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTa                                                                                                                              <  1:789838/149‑1 (MQ=255)
                    aGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTa                                                                                                                              >  2:789838/1‑149 (MQ=255)
                       cTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGTCTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTagcg                                                                                                                           >  2:976268/1‑149 (MQ=255)
                            aCAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTagcgccagc                                                                                                                      >  4:311571/1‑149 (MQ=255)
                                 gCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCAc                                                                                                                 <  3:779780/149‑1 (MQ=255)
                                     gaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            >  2:1128669/1‑118 (MQ=255)
                                     gaAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGtt                                                                                                                                            <  1:1128669/118‑1 (MQ=255)
                                         gCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAgcgg                                                                                                         <  1:976268/149‑1 (MQ=255)
                                                ttGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGt                                                                                                  <  4:73786/149‑1 (MQ=255)
                                                ttGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGt                                                                                                  >  1:895913/1‑149 (MQ=255)
                                                   tATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaa                                                                                               >  3:220347/1‑149 (MQ=255)
                                                     ttttttGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTaaaaa                                                                                             >  2:394630/1‑149 (MQ=255)
                                                                        agagAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGc                                                                                                                                   <  2:377573/92‑1 (MQ=255)
                                                                        agagAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGc                                                                                                                                   >  1:377573/1‑92 (MQ=255)
                                                                         gagaACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGaaaaa                                                                         <  3:473279/149‑1 (MQ=255)
                                                                              ccaccGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGt                                                                    >  4:687022/1‑149 (MQ=255)
                                                                                accGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGttt                                                                  <  4:209772/149‑1 (MQ=255)
                                                                                  cGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTg                                                                >  4:638112/1‑149 (MQ=255)
                                                                                   gTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGa                                                               <  1:282328/149‑1 (MQ=255)
                                                                                     cATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatc                                                             <  1:1051796/149‑1 (MQ=255)
                                                                                      aTTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatca                                                            <  3:810112/149‑1 (MQ=255)
                                                                                       ttttCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcat                                                           >  4:1056523/1‑149 (MQ=255)
                                                                                        tttCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGatcatc                                                          <  2:372183/149‑1 (MQ=255)
                                                                                            caccacATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCa                                                      >  4:412064/1‑149 (MQ=255)
                                                                                              ccacATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAAc                                                                                                    >  2:146412/1‑101 (MQ=255)
                                                                                              ccacATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAAc                                                                                                    <  1:146412/101‑1 (MQ=255)
                                                                                                  aTGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATc                                                >  1:84339/1‑149 (MQ=255)
                                                                                                        tttGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCataat                                          >  2:1087628/1‑149 (MQ=255)
                                                                                                            ttttCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaat                                      <  2:895913/149‑1 (MQ=255)
                                                                                                             tttCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATaataaata                                     <  3:412064/149‑1 (MQ=255)
                                                                                                                cGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATaa                                  <  1:667569/149‑1 (MQ=255)
                                                                                                                cGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATaa                                  <  3:311571/149‑1 (MQ=255)
                                                                                                                 gCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAg                                 <  1:1087628/149‑1 (MQ=255)
                                                                                                                             acTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGACACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTCTTGtaata                     <  3:1056523/149‑1 (MQ=255)
                                                                                                                              cTGGCCGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGtaataa                    >  3:347708/1‑149 (MQ=255)
                                                                                                                               tggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGt                                                                                         >  2:941404/4‑79 (MQ=38)
                                                                                                                               tggccGCTGCCGCCAGCGGGCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGt                                                                                         <  1:941404/76‑1 (MQ=38)
                                                                                                                                 cccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACa                 <  1:27948/149‑1 (MQ=255)
                                                                                                                                  ccGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCGTCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAg                >  3:302302/1‑149 (MQ=255)
                                                                                                                                        gcggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGc          <  4:220347/149‑1 (MQ=255)
                                                                                                                                          ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt        <  3:687022/149‑1 (MQ=255)
                                                                                                                                          ggcAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCt        >  3:143771/1‑149 (MQ=255)
                                                                                                                                              gCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTAtt                           >  1:1029680/1‑126 (MQ=255)
                                                                                                                                              gCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTAtt                           <  2:1029680/126‑1 (MQ=255)
                                                                                                                                                ggCTAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTa  >  4:975793/1‑149 (MQ=255)
                                                                                                                                                  ccAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAAc                                                                                           <  2:791207/58‑1 (MQ=255)
                                                                                                                                                  ccAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAAc                                                                                           >  1:791207/1‑58 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAAGGATCAACGCCGGGGTTAGCCTGTAACAGAGCGAGAAAGCCGACGTTGTATTTTTTGGCGATGGCTTCCAGAGAACCACCGTCATTTTCCACCACATGAAATTTGTTTTCGCCAACCAGACGACTGCCCGCTGGCGGCAGCGGCCAGGTGTTGGCTTTTGCCGGTAGCGCCAGCGCCACGGCAGCGGCGAACGTAAAAAACGTTAGCCAGCGAGAAAAACGCGTTTTGATCATCACCAAAAATCCATAATAAATATAAGGTTATTGTAATAAACAGATAAGCCTGAATTA  >  NC_000913/1170136‑1170428

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